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ABSTRACT 


Based  on  a  modified  output  regulator  problem,  a  design  oriented 
methodology  is  presented  for  the  construction  of  output  feedback  compensa¬ 
tors  retaining  optimal  eigenvectors  from  a  reference  state 

feedback  regulator.  Viewing  &  as  a  design  parameter,  it  is  known  that  in 
the  case  &>r  this  requires  a  dynamic  compensator  of  dimension  &  -  r  whose 
parameters  are  determined  in  function  of  the  solution  of  an  associated 
output  feedback  pole-placement  problem.  Using  an  iterative  dyadic  pole- 
placement  procedure,  an  algorithm  is  given  which  determines  the  solution 
of  this  pole-placement  problem  without  a  priori  assumptions  on  the  compen¬ 
sator  dimension.  The  methodology  is  also  extended  to  the  class  of  stabiliz- 
able  systems  and  the  required  compensator  shown  to  posse®-  a  separation 
property.  Finally  the  design  methodology  is  illustrated  by  three  nontrivial 
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INTRODUCTION 

In  the  design  of  linear  quadratic  regulators  for  time-invariant 
systems,  the  inaccessibility  of  state  variables  precludes  the  implementation 
of  optimal  state  feedback  control  laws.  Two  approaches  to  the  resolution 
of  this  difficulty  are  the  reconstruction  of  unmeasured  states  by  reduced 
order  observers  and  the  reformulation  of  the  optimization  problem  as  an 
output  feedback  quadratic  regulator  problem.  The  introduction  of  reduced 
order  observers  allows  the  retention  of  the  optimal  state  feedback  control 
law,  but  may  often  result  in  the  use  of  a  compensator  of  much  higher  dimen¬ 
sion  than  is  actually  needed  to  satisfactorily  control  the  system.  The 
reformulation  of  the  problem  under  output  feedback  suffers  a  more  severe  dis¬ 
advantage.  Whereas  there  is  a  wealth  of  literature  available  on  the  state 
feedback  regulator  and  associated  state  reconstruction  problems,  or  estima¬ 
tion  and  filtering  problems  in  the  presence  of  noise,  very  little  is  known 
about  the  existence  and  properties  of  solutions  to  the  output  feedback  reg¬ 
ulator  problem.  It  is  therefore  of  interest  to  develop  implementable 
regulators  which  retain  some  measure  of  optimality  provided  by  the  state 
feedback  regulator  problem,  without  requiring  the  use  of  compensators  of 
high  dimension. 

One  solution  to  this  problem  is  to  design  compensators  for  output 
feedback  regulators  which  retain  as  many  optimal  eigenvectors  of  the  corre¬ 
sponding  state  feedback  regulators  as  possible.  Such  regulators  have  the 
properties  of  achieving  the  optimal  cost  in  the  subspace  spanned  by  the 
retained  eigenvectors,  and  of  providing  an  easily  computed  measure  of  the 
cost  degradation  in  the  remaining  state  space.  It  is  known  that  retention 
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of  r  +  p  optimal  eigenvectors,  where  r  is  the  rank  of  the  output  matrix, 
requires  the  construction  of  a  dynamic  compensator  of  dimension  p,  and  that 
the  compensator  parameters  may  be  determined  in  function  of  the  solution  of 
an  associated  output  feedback  pole -placement  problem  [1]. 

Following  a  review  of  a  design  oriented  methodology  for  the  con¬ 
struction  of  output  feedback  regulators  which  retain  r  +  p  optimal  eigen¬ 
vectors  (piO),  an  algorithm  will  be  presented  which  solves  the  associated 
pole-placement  problem  end  determines  the  dimension  p  of  the  required  com¬ 
pensator  without  a  priori  assumptions.  In  the  event  the  system  is  stabili- 
zable  by  static  output  feedback,  rather  than  construct  a  dynamic  compensator, 
it  may  be  preferable  to  relax  the  requirement  of  retention  of  r  +  p  optimal 
eigenvectors.  The  problem  of  retention  of  fewer  than  r  optimal  eigenvectors 
will  therefore  be  considered  and  shown  to  also  give  rise  to  an  output  feed¬ 
back  pole-placement  problem.  Finally  it  will  be  shown  that  the  design 
methodology  may  be  extended  to  the  class  of  stabilizable  systems. 

In  chapter  one  dyadic  solutions  to  the  static  output  feedback 
pole-placement  problem  are  reviewed.  The  second  chapter  presents  the 
methodology  for  the  design  of  suboptimal  linear  quadratic  regulators  which 
retain  X(lsi^n)  optimal  eigenvectors  from  the  state  feedback  regulator. 
Based  on  the  methodology  proposed  in  [2]  for  solving  the  general  output 
feedback  pole-placement  problem,  an  algorithm  is  obtained  in  chapter  three 
which  solves  the  pole -placement  problem  associated  with  the  design  of  sub- 
optimal  regulators.  In  the  last  chapter  three  nontrivial  examples  illustrate 
the  design  methodology. 


CHAPTER  1 


DYADIC  SOLUTIONS  TO  THE  OUTPUT  FEEDBACK  POLE -PLACEMENT  PROBLEM 

Since  the  problem  of  eigenvalue  assignment  by  state  feedback  was 
resolved  [3]  the  related  problem  of  eigenvalue  assignment  by  output  feed¬ 
back  has  been  the  subject  of  extensive  research.  Aside  from  numerical  methods, 
approaches  to  the  problem  may  be  characterized  as  giving  rise  to  either  dyadic 
or  full  rank  feedback  matrices. 

In  the  dyadic  approach  the  nonlinear  equations  which  describe  the 
complete  solution  to  the  problem  are  rendered  linear  by  arbitrarily  fixing 
certain  otherwise  free  parameters,  and  the  feedback  gain  matrix  is  obtained 
as  a  sum  of  dyadic  products.  Dyadic  solutions  have  been  obtained  for  systems 
represented  in  state  space  [4],  [5],  [6],  [7]  and  for  systems  described  by 
transfer  functions  [8],  [9],  [10],  [11]. 

Other  approaches  to  the  problem  attempt  to  utilize  the  freedom 
discarded  in  the  dyadic  approach,  typically  to  assign  eigenvectors  as  well 
as  eigenvalues,  and  while  usually  resulting  in  full  rank  feedback  matrices, 
tend  to  give  rise  to  solutions  which  are  numerically  difficult.  Extensive 
results  have  been  obtained  using  geometric  approaches  [12],  [13],  [14]. 

Other  approaches  include  the  generalized  root  locus  [15],  [16],  [17],  the 
use  of  generalized  inverses  [18],  [19]  and  the  use  of  Kronecker  products  [20], 

Because  of  the  difficulty  of  obtaining  explicit  solutions  to  the 
output  feedback  pole  placement  problem,  many  numerical  procedures  have  been 
proposed.  In  view  of  the  fact  that  in  general  only  m  +  r-1  eigenvalues  may  be 
arbitrarily  assigned  while  nothing  can  be  said  about  the  resulting  locations 
of  the  remaining  eigenvalues,  it  has  been  suggested  that  the  problem  be 
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reformulated  as  an  optimization  problem  to  minimize  the  deviations  of  all  n 
eigenvalues  from  their  desired  locations  [21],  [22],  Other  authors  have 
considered  the  problem  of  requiring  the  closed  loop  eigenvalues  to  lie  in 
prescribed  regions  of  the  complex  plane  [23],  [24].  This  leads  naturally  to 
recasting  the  problem  as  an  output  feedback  linear  quadratic  regulator  prob¬ 
lem  [24] ,  but  theoretical  results  are  lacking.  In  chapter  2  one  approach 
to  designing  suboptimal  linear  quadratic  regulators  will  be  seen  to  lead 
directly  back  to  the  output  feedback  pole-placement  problem. 

In  this  chapter  dyadic  solutions  to  the  output  feedback  pole- 
placement  problem  are  discussed  in  detail  as  they  will  be  used  in  chapter  3 
to  solve  the  pole-placement  problem  associated  with  the  design  of  suboptimal 
regulators. 

Let  the  triple  (A,B,C)  represent  the  linear  time- invariant  system 

x  =  Ax  +  Bu  ,  y  =  Cx  ,  A€Rnxn,  B€RnXn,  C6Rrxn  (1.1) 

where  B  and  C  are  assumed  to  be  of  full  rank.  Then  except  for  certain 
singular  systems  it  is  known  that  max(n,m  +  r-l)  eigenvalues  may  be  assigned 
"almost"  arbitrarily  by  output  feedback  of  the  form  u=Ky;  K£Rmxr.  The 
proof  of  this  is  constructive  and  relies  on  the  fact  that  uncontrollable 
and/or  unobservable  eigenvalues  of  a  system  are  invariant  under  output  feed¬ 
back.  The  idea  is  that,  after  an  initial  feedback  placing  some  eigenvalues 
at  their  desired  locations,  the  system  may  be  collapsed  to  a  single  input  or 
output  system  having  min(m,r)-l  of  these  eigenvalues  uncontrollable  or  un¬ 
observable.  The  remaining  min(max(m,r),  n-min(m, r) + 1)  desired  eigenvalues 
may  then  be  assigned  by  a  further  feedback.  Formulas  for  the  computation  of 


the  required  gains  are  given  in  the  proof  of  Theorem  1.2  below,  and  Theorem 
1.3  gives  the  final  desired  gain  as  a  sum  of  dyadic  products. 


The  following  definitions  make  precise  the  use  of  the  word  "almost". 
Definition  1.1  f 2 61  :  Let  [^(x)}^^  be  a  finite  set  of  real  valued  poly¬ 
nomials  taking  their  arguments  in  a  parameter  space  Rn.  Then  the  set  of 
common  zeros  V  =  [x  €  Rn  :  ^(x)  =0,  l£i£Nj  is  called  a  proper  variety  pro¬ 
vided  V  #  Rn.  D 

Definition  1.2  f 2 61 :  A  property  is  said  to  hold  for  almost  all  points  in  a 
parameter  space  Rn  provided  the  set  of  points  at  which  the  property  fails  to 
hold  is  the  union  of  a  finite  number  of  proper  varieties.  D 

Also  are  needed: 

Definition  1.3  [131:  A  set  T ^  of  SL  complex  numbers  is  said  to  be  a  symmetric 
set  provided  X  €  if  and  only  if  X*  €  T Cl 
Definition  1.4:  The  matrix  A  is  said  to  have  l  assignable  eigenvalues  if  for 
almost  all  symmetric  sets  T ^  there  exists  an  output  feedback  matrix  K  such 
that  rxca(A  +  BKC). 

A  preliminary  result  will  permit  the  avoidance  of  any  discussion 
of  the  complications  which  arise  if  the  matrix  A  has  repeated  eigenvalues. 
Theorem  1.1  (271 :  If  (A,B,C)  is  a  controllable  and  observable  triple  then 
for  almost  all  matrices  K,  the  eigenvalues  of  A  +  BKC  will  be  distinct. 

Proof;  See  [27] .  O 

The  proof  of  the  next  theorem  provides  formulas  for  computing  the 
factors  of  dyadic  output  feedbacks.  A  result  that  will  be  needed  in  the  proof 

is  contained  in  Lemma  1.1.  To  simplify  the  notation,  let  I A I  denote  the  determinant 
of  A. 

d  x  1  t 

Lemma  1,1  f 281  :  If  x,y€R  are  column  vectors  then  |ln+xy  |  =(l  +  y  x). 
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Proof:  Let  V€R  ^  '  have  as  columns  any  basis  for  the  null  space  of  y  . 

T  T 

Then  (I  +  xy  )V*V  and  it  follows  that  1  is  an  eigenvalue  of  I  +xy  with 
n  n 

algebraic  multiplicity  n  -  1.  Letting  X  denote  the  remaining  eigenvalue  of 

T 

I  +  xy  then 
n 

lln  +  xyT|=\  ,  trace  (In +  xyT)  »  X  +  (n-1)  (1.2) 

\ 

T 

and  since  trace(In  +  xy  )  =  n  +  <x,y>  there  follows 

lln  +  xyT|  =  trace(In  +  xyT)  -  (n-1)  =  (l+yTx)  (1.3)  □ 

Theorem  1.2  f41,f6]:  If  (A,B,C)  is  a  controllable  and  observable  triple 
with  C  having  full  rank  r  and  A  having  distinct  eigenvalues,  then  for  any 
symmetric  set  Tr  there  exists  a  feedback  matrix  K  such  that  r  eigenvalues 
of  A  +  BKC  are  arbitrarily  close  to  the  elements  of  the  set  1^. 

Proof :  Let  Tr  *  [X^j^_^  a  symmetric  set  of  complex  numbers  to  be  assigned 
to  the  spectrum  of  A  +  BKC.  By  Theorem  1.1  it  may  be  assumed  without  loss  of 
generality  that  the  eigenvalues  of  A  are  distinct.  The  theorem  is  proved  by 
finding  a  vector  f  such  that  the  pair  (A,Bf)  is  controllable  and  then  solving 
the  pole-placement  problem  for  the  single  input  system  (A,Bf,C).  The  final 
feedback  gain  will  be  a  dyadic  product  K  =  fg. 

Let  T€RnXn  transform  the  matrix  A  to  Jordan  form  so  that 

T_1AT«dg[ai  ...oj  ,  T_1B  =  [T'1b1  \  ...  \  T_lb  ]  (1.4) 

and  note  that  by  the  assumption  of  controllability,  each  row  of  T  has  at 
least  one  nonzero  entry.  Select  f 6  Rm  X ^  such  that  no  entry  of  T  ^Bf  is 
zero.  This  is  always  possible,  as  the  set  of  vectors  f  such  that  some  entry 
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Defining  P-Ip^)  ...  PQttr)]  €  R*  *  l»nd  S -[s^)  l  ...  j  s(Xr)]  €Rnxr 
and  constraining  pc(X^)  *0  for  i“l, ...,r  results  in  the  square  linear 
system  of  equations 

gCQRS  -  P  (1.10) 

which  will  always  have  a  solution  g  except  for  those  choices  of  Fr  for  which 
IcQRsl  =0.  If  this  determinant  is  zero  and  the  equations  are  inconsistent 
then  an  arbitrarily  small  perturbation  of  the  X^  will  result  in  a  consistent 
set  of  equations,  but  it  should  be  noted  that  in  practice  this  will  result 
in  arbitrarily  large  gains  in  some  or  all  of  the  entries  of  g.  The  theorem 
is  proved  on  writing  the  final  feedback  as  K  =  fg.  □ 

By  duality  there  follows  immediately  from  Theorem  1.2: 

Corollary  1.2.1  f41:  If  (A,B,C)  is  a  controllable  and  observable  triple  and 
B  and  C  are  of  full  rank,  then  for  any  symmetric  set  r  of  p=max(m,r)  com¬ 
plex  numbers  there  exists  a  matrix  K  such  that  the  eigenvalues  of  A  +  BKC 
are  arbitrarily  close  to  the  elements  of  r  .  Q 

It  should  be  noted  that  the  selection  of  f  in  the  proof  of  the 
theorem  was  completely  arbitrary  and  represents  precisely  that  loss  of 
freedom  in  the  specification  of  a  dyadic  feedback  K  which  results  in  a 
linear  system  of  equations.  Also  the  requirement  of  controllability  of 
the  pair  (A,B)  may  be  relaxed  to  the  requirement  that  (CQRS)  be  invertible, 
but  conditions  on  the  matrices  A,  B,  f  under  which  this  will  be  true  are  not 
known.  That  the  matrix  may  fail  to  be  invertible  for  some  uncontrollable 
systems  will  lead  to  the  failure  of  Theorem  1.3  to  hold  for  all  triples 
(A,B,C).  Finally  it  is  remarked  that  the  method  of  the  proof  of  Theorem  1.3 


J 
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fails  in  the  case  Tr  contains  repeated  eigenvalues.  A  straightforward 
modification  of  the  equations  to  handle  this  case  is  given  in  [4]  and  may 
be  of  interest,  for  example,  in  solving  minimum  time  problems  for  discrete 
systems. 

In  general, for  any  fixed  f  the  equation  IcQRsI  =0  defines  a  prop¬ 
er  variety  of  eigenvalues  T  which  are  not  assignable  by  any  finite  output 
feedback.  The  following  example  illustrates  this  point. 

Example  1.1  f!21;  Let  (A,B,C)  be  given  by 


(1.11) 


3  2 

The  open  loop  characteristic  equation  is  pQ(\)=X  +  X  -1. 
f»(fpf2)T  the  gain  K  =  f g  is  defined  by  (1.10): 

0  f. 


(8j.»  82) 


0 

1 


fl  ‘  f2 


0 

f. 


3  2  3  2 

<X1  +  X1  "  X2+X2  '  1} 


For  any  given 


(1.12) 


and  the  equation 


0-  ICQRSI  *  (\2  -X1)(fJ(X1X2+X1  +  X2  +  l)  - 

f^(XL  +  X2  +  1)  -fLf2)  (1.13) 

defines  those  symmetric  sets  (Xj^Xj)  which  are  not  assignable,  though  as 
noted  above,  the  equations  may  be  modified  to  handle  the  case  X^*X2*  a 
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The  main  theorem  may  now  be  given. 

Theorem  1.3  f 51 . T 71 !  For  almost  all  controllable  and  observable  triples 

(A,B,C)  and  for  any  symmetric  set  of  p  =min(n,m  +  r-l)  complex  numbers, 

there  exists  a  matrix  K  such  that  the  eigenvalues  of  A  +  BKC  are  arbitrarily 

close  to  the  elements  of  F  . 

P 

This  result  was  obtained  in  [5],  [7]  and  [12]  and  provides  the  best 
known  bound  on  the  number  of  eigenvalues  assignable  by  output  feedback.  For 
a  geometric  proof  which  constructs  a  full  rank  feedback  the  interested 
reader  is  referred  to  [12], 

Proof:  Let  r  —  j.— ^  be  a  symmetric  set  of  p  =min(n,m  + r-1)  complex  numbers 

to  be  assigned  to  the  spectrum  of  A  +  BKC.  By  Theorem  1.1  it  may  be  assumed 
that  the  eigenvalues  of  A  are  distinct.  The  theorem  will  be  proved  by 
constructing  a  rank  2  gain  K*fjg^  +  f2g2  where  f^  will  be  arbitrary,  g^  chosen 
to  place  at  least  r-1  eigenvalues,  g2  to  render  those  r-1  eigenvalues  un¬ 
observable,  and  t0  place  the  remaining  min(n-r+l,m)  eigenvalues. 

With  A,  b,  c,  ta^  ,  {a^}  ^  >  P*(X),  Pq(\)  to  be  defined 

below  let: 
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,  n-1 


>n-l 


min(n,mfr-l) 


Xmin(n,mfr-1)  X 


(1.14) 


Pl’IPoaL) 


P2  "  Ipo(XP  **•  po(Xmin(n,nri-r-l))] 


„mx  1 


As  in  the  proof  of  Theorem  1.2  let  f^€R  be  any  vector  such 
that  the  pair  (A,b),  b^Bf^is  controllable,  and  denote  the  open  loop 
characteristic  polynomial  of  A  by 


1...  y  l.i  1  , 

Po(X)  “  iSoaiX  »  an  =  1 


(1.15) 


1  x  r 


Then  by  Theorem  1.2,  the  solution  g^ 6 R  of  the  equation 

M0  Wi1  ’  pi 


(1.16) 


will  assign  the  eigenvalues  to  t*16  spectrum  of  A=A  +  Bf^g^C,  subject 

to  a  possible  perturbation  of  the  numbers  to  ensure  the  consistency  of  the 
equations. 


Recalling  that  unobservable  eigenvalues  are  invariant  under  output 
feedback  let  V^€RnX^r  ^  and  r+^  have  as  columns  the  eigen- 

A  .  T“1 

vectors  of  A-  corresponding  respectively  to  anc^  c^e  remaining  eigen¬ 

values  of  A.  Then  if  g2CV^*0,  the  single  output  system  (&,B,c),  c4g2C, 
will  have  r-1  unobservable  eigenvalues  .  The  number  of  nonzero 

A 

entries  of  the  vector  82^2  t*ie  nuin^er  of  eigenvalues  of  A  subject 

to  influence  under  further  feedback.  A  solution  to  82^2*^  always  exists 
as  it  requires  finding  an  r-vector  orthogonal  to  r-1  other  r-vectors,  how¬ 
ever  conditions  under  which  at  least  min(m,n-r+l)  entries  of  ^2^2 


A 
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[S2R2Q2B]£  2  =  P2  (1.18) 

will  assign  the  remaining  eigenvalues  l  to  the  spectrum  of 

A 

A  +  Bf2g2C,  provided  the  equations  are  consistent.  The  final  dyadic  feed¬ 
back  will  then  be  K  =  f-jg^  +  f2g2.  If  the  equations  are  inconsistent,  then 
a  perturbation  of  the  eigenvalues  to  be  assigned  may  possibly  fail  to  make 
invertible,  since  the  factor  R.,Q2B  cannot  be  guaranteed  to  be  of 
full  rank,  the  triple  (A,B,C)  not  being  observable.  □ 

Since  rank-deficient  compensators  have  poor  disturbance  rejection 
properties  in  applications,  the  procedure  of  this  proof  may  be  iterated  to 
obtain  full  rank  feedback  matrices  as  sums  of  min(m,r)  dyadic  products  [2]. 
This  modification  provides  the  basis  for  the  algorithm  given  in  section 
3.1  and  will  be  discussed  in  detail  there. 

The  following  example  illustrates  the  possibility  of  a  system 

structurally  failing  to  allow  the  assignment  of  min(n,nrt-r-l)  eigenvalues. 

2x2  2x1 

Example  1,2;  Let  ApA^  €  R  and  bpb2 ’CpCj  6  R  and  consider  the 
system 
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1 


2  2 

Under  output  feedback  u=Ky,  K€R  the  closed  loop  system  is 


Al  +  kllblCl 

k21b2cl 


k12blc2 


A2  +k22b2C2 


(1.20) 


« 

♦ 

; 5  \ 


<* 


The  result  of  Theorem  1.3  would  predict  that  min(n,nrt-r-l)  = 3  eigenvalues 

T  T 

may  be  assigned.  However  if  either  b^c^  =0  or  b^c^  =  0  the  closed  loop 
system  is  block  triangular  and  only  max(m,r)  =2  eigenvalues  may  be  assigned, 
one  to  each  diagonal  block  by  the  choice  of  k^,  ^22*  ° 

An  alternate  characterization  of  the  vectors  f^,g^  in  the  proof  of 
Theorem  1.3  is  available.  Rather  than  requiring  pc(Xi)=0,  i =  1, . . . ,max(m,r) 
in  the  proof  of  Theorem  1.2,  the  desired  closed  loop  characteristic  equation 
may  be  constrained  to  be  of  the  form 

max(m,r) 

pctt)  -  Pla)P2(X),  (x-xt)  (i.2i) 

giving  rise  to  a  system  of  equations  which  determine  not  only  the  feedback 
gain  but  also  the  remaining  spectrum  of  the  closed  loop  system  as  the  roots 
of  P2(X)  [8].  Since  this  will  require  the  computation  of  the  Markov  param¬ 
eters  CA^B,  a  pole-placement  procedure  incorporating  this  approach  will  be 
referred  to  in  this  thesis  as  a  frequency  domain  solution,  whereas  the  pro¬ 
cedure  contained  in  the  proofs  of  Theorems  1.2  and  1.3  will  be  referred  to 
as  a  state  space  solution.  By  duality  it  suffices  to  define  a  frequency 
domain  solution  by  indicating  formulas  for  pole-placement  in  single-input 
systems  and  for  rendering  poles  uncontrollable. 

Let  the  system  (A,B,C)  have  transfer  function 


14 


1  n  i 
C(XI-A)"b  =  N (X)  /  d(X)  ,  d(X)  =  i^0  a.X1  ,  an  =  1  (1.22) 

and  assume  that  the  eigenvalues  of  A  are  distinct.  Using  Leverrier's 
algorithm  the  transfer  function  may  be  written  in  terms  of  the  Markov 
parameters  CAXB: 


n-1  i  .... 

N(X)  =  C  adj  (XI -A)  S  =  C(.2q  Zq  a^A1'^"'1'  )B 
n-1  n-k-1  ,  ,  . .  .  , 

-  Jo  <j5o  Vj“  '  ™ 

n-k-1  n-k-1- ‘ 

Defining  »  ,2q  a  ^CAn  ^B  equation  (1.22)  becomes 


C(XI-A)_1B  =  (i^Q  N.X1)  /  d(X) 


n-1 


(1.23) 


(1.24) 


The  following  lemma  is  needed  in  obtaining  single  input  systems 

with  prescribed  uncontrollable  eigenvalues. 

Lmima  1.2  fill  :  If  Xq  is  a  root  of  d(X)  then  rank  N(Xo)sl. 

Proof:  Let  T  transform  A  to  the  Jordan  form  T  *AT  = dg(X. , . . . ,X  ).  Then 
-  l  n 


-1. 


N(X)  =»  C  adj  (XI -A) B  =  CT  adj  (dg(X-X1, . . .  ,X"^n))T  B 


-1 


(1.25) 


=  CT  dg(.TJ1(X-X.  ),...,  ^n(X-X.))T  B 


If  X  =*X.  for  some  k,  l£ks:n,then  N(X  )  becomes 
ok  o 


-1„ 


N(Xq)=CT  dg(0,...,0,i^k(X-Xi),0,...,0)T  B 


(1.26) 


and  so  rank  N(X  )  SI.  □ 
o 

Recall  that  a  pole  of  a  single  input  system  is  uncontrollable  if 
the  numerator  of  the  transfer  function  is  zero  evaluated  at  that  pole,  and 
let  f  define  the  single  input  system  (A,Bf,C)  with  transfer  function 
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H  (X)  =  N(X)f  /  d(X) 


N(X)  = 


n-1 

E 

i=0 


NA 

L 


>Ni 


n  . 


n 


gRnxi 


(1.27) 


To  render  m-1  poles  {X^}™_^  uncontrollable  it  suffices  to  satisfy  N(X^)f  =  0, 
i  =  1, . . . ,m-l.  If  the  matrix  M  €  Rm  lxm  is  taken  to  have  as  its  ith  row  the 
linearly  independent  row  of  N(X^)  then  f  may  be  obtained  as  the  solution  of 
Mf  =  0  [2],  [11].  (If  N(X^)  =0  then  X^  is  already  uncontrollable  and  the 

ith  row  may  be  taken  all  zeros). 

*  r  X  x  it 

Let  ke  t^ie  r  P°les  t0  be  assigned  by  the  feedback  g€R 

of  Figure  1.1.  By  summing  at  nodes  1  and  2  in  the  figure  the  closed  loop 

transfer  function  may  be  written 


Hc(\)  =  (I  -H(\)g)-1H(X)  =  H (X)  /  (1  -  gH(X)  )  Cl- 28) 


from  which  the  closed  loop  polynomial  is  given  by 


PC(X)  =  d(\)  -  gN (X) 


(1.29) 


Node  1  Node  2 


Figure  1.1.  Block  diagram  of  transfer  function  given  in 
(1.28). 
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These  equations  determine  the  gain  g,  and  the  remaining  spectrum  of  A  +  BfgC 
as  the  roots  of  p2(X).  As  in  Theorem  1.2  it  may  arise  that  the  coefficient 
matrix  in  (1.35)  is  singular.  If  the  system  of  equations  is  inconsistent 


then  an  arbitrarily  small  perturbation  of  the  X^  will  render  the  coefficient 
matrix  invertible  by  altering  the  lower  (n-r)xn  block. 

To  illustrate  frequency  domain  and  state  space  solutions  of  the 
pole-placement  problem  reconsider  Example  1.1. 

Example  1.1  (continued)  f 12 1 :  Let  A,B,C  be  given  by  (1.11),  and  assume  the 
desired  closed  loop  spectrum  is  (-2,-1+  lj).  Since  m  +  r-l  =  3,  arbitrary 
pole-placement  is  possible. 

A  dyadic  feedback  K  =  +  comPute<*  *-n  two  stages.  At 

the  first  stage  two  poles  are  placed  and  at  the  second  stage  one  of  these 


I  poles  is  rendered  invariant  to  further  feedback  and  two  additional  poles 

are  assigned.  There  are  several  possibilities  for  computing  f^,g^,f2»g2- 


At  the  first  stage  f^  may  be  chosen  arbitrarily  and  g^  computed  from  (1.16) 
to  place  r®2  poles,  or  g^  may  be  selected  arbitrarily  and  f^  computed  from 
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(1.18)  to  place  m  =  2  poles.  At  the  second  stage  there  are  also  two  choices. 
Either  f2  may  be  chosen  to  define  a  single-input  system  with  m-1  uncontrol¬ 
lable  poles  and  g2  computed  from  (1.16)  to  place  r  additional  poles,  or 
may  be  chosen  to  define  a  single-output  system  with  r-1  unobservable  poles 
and  f2  computed  from  (1.18)  to  place  m  additional  poles.  In  this  example, 
for  both  the  state  space  and  the  frequency  domain  solutions,  r  poles  will 


be  assigned  at  the  first  stage.  At  the  second  stage  g2  will  be  chosen  to 

render  r-1  poles  unobservable  and  then  f2  will  be  computed  to  place  the 

remaining  m  poles  at  their  desired  locations.  Because  the  complex  pair 

-1+lj  may  not  be  split,  \  =  -2  must  be  assigned  at  the  first  stage. 

The  state  space  solution  is  as  follows.  Arbitrarily  select 
T 

f^ = (1,0)  and  let  the  spectrum  to  be  assigned  by  g^  be  (-2,0).  With 
1  3  2 

PqU)=X  +  X  -1,  ^  =  -2  ,  \2  =  0  equation  (1.16)  becomes 


=  (-5,-1) 


(1.36) 


which  has  solution  g^  =  (-l,-3).  The  resultant  spectrum  is  a(A)  =  (0,-2, -2) 
and 


-10  1 

A=  1-3-1  ,  Pq(\)  =  \3+4X2+4\ 


(1.37) 


0  10 
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T 

The  repeated  eigenvalue  X  = -2  has  only  one  eigenvector  V^*(-l,-2,l)  . 

To  render  one  eigenvalue  at  -2  unobservable  the  equation  g^CV^ *  0  is 

2  - 

solved  for  g2*'(2,l).  Using  pQ(-l  +  lj)  *  -2  +  2 j  equation  (1.18)  becomes 


“2 j  -1+J  1 

L?j  -i-J  1 


10  0 
4  10 

4  4  1 


0  12 
1  -1  -1 

-2  2  2 

“2  -2  j 
-2+2  j 


0  0 
1  0 
0  1 


(1.38) 


Premultipling  by 


1  1 
2  2 

-ij  lj 

2J  2J 


to  obtain  a  set  of  real  equations,  the  solution  is 


f2  =  (4,-2)  .  The  final  feedback  is  then 


7  1 

-4  -2 


(1.39) 


and  the  spectrum  of  A  +  BKC  is  (-2,-1+lj). 

The  frequency  domain  solution  is  as  follows.  Again  selecting 
T 

fL« (1,0)  the  transfer  function  for  the  single  input  system  (A,Bf^,C)  is 


H(X) 


x3  +  x2-~i 


X  +1 

X2  +x 


(1.40) 


Placing  poles  at  -2,0,  (1.31)  becomes  p1(X)=X2  +  2X  and  (1.35)  gives 


tgi  :  h3 


i 

i 


0 

1 


0  2 


(-1,0,1)  -  (0  :  0,2)  (1.41) 


i 
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The  solution  is  g^=*(-l,-3)  and  h  =  2.  The  spectrum  of  A  is  therefore  (0,-2) 
together  with  the  root  -2  of  p2(X)  *X  +  h.  Since  the  spectrum  of  A  is  known, 
its  characteristic  equation  is  also  known  and  the  transfer  function  for  the 
system  (A,B,C)  may  be  computed  from  the  Markov  parameters  (1.23): 

C(Xl-A)B  =  N(X)/d(X)  =  -r - - - 

X  +  4XZ  +  4X 


X+l  X  +  4\  +  3 
-X 


x2  +x 


(1.42) 


To  render  the  pole  at  -2  unobservable,  g 2  is  computed  as  the  solution 
g2  *  (2,1)  of  g2N(-2)  *0: 


g2N(-2)  =  g2 


-1  -1 

2  2 


(1.43) 


The  transfer  function  of  the  single-output  system  (A,B,g2C)  is  then 


H(X)  = 


X3+4X2+4X 


(X2  +3X  +  2,2X2  +  7X  +  6) 


(1.44) 


To  place  two  poles  at  -1  +  lj,  (1.31)  becomes  p2  (X)  =  X  +  2\  +  2 ,  and  the 
dual  of  (1.35)  gives 


2  2 
3  7 

1  6 


2 

2 

1 


0 

4 

4 


0 

2 

2 


(1.45) 


The  solution  is  f2  *  (4,-2)  ,  h=2.  The  final  feedback  gain  K  is  given  by 
(1.39),  and  the  spectrum  of  A  +  BKC  is  -1  +  lj  together  with  the  root  -2  of 
P2(X)  =X  +  h.  q 


i 
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CHAPTER  2 

DESIGN  OF  SUBOPTIMAL  LINEAR  QUADRATIC  REGULATORS 

This  chapter  considers  the  problem  of  designing  suboptimal  linear 
quadratic  regulators  having  the  property  of  retaining  any  number  X 
(lsXsn)  of  eigenvectors  at  their  optimal  locations  as  defined  by  a 
reference  state  feedback  regulator  problem.  The  system  under  consideration 
will  be  taken  to  be 

x  =  Ax  +  Bu  ,  y  =  Cx  ,  A  6  Rnxn,  B  €  R11*111,  C  €  Rrxn  (2 . 1) 

and  will  be  assumed  to  be  controllable  and  observable. 

In  the  case  X^r  the  solution  will  be  seen  to  be  given  by  a  static 
output  feedback  compensator,  while  in  the  case  X> r  a  dynamic  compensator 
will  be  required.  In  the  complete  methodology  X  is  viewed  as  a  free  para¬ 
meter  to  be  specified  during  the  design  procedure  and  not  fixed  a  priori. 
When  X  r  there  exists  freedom  in  the  design,  which  is  translated  into  the 
choice  of  feedback  gains  if  X  <  r  and  into  the  choice  of  parameters  of  the 
dynamic  compensator  if  X>r.  In  both  cases  this  freedom  is  used  to  shape 
the  complementary  spectral  characteristics  of  the  closed  loop  system  by 
solving  an  associated  output  feedback  pole-placement  problem,  a  solution  to 
which  will  be  given  in  section  3.1. 

2.1.  Review  of  Necessary  Conditions 

In  this  section  solutions  to  three  related  linear  quadratic 


regulator  problems  are  presented  and  their  properties  reviewed. 
Let  Q  €  Rnxn,  R  €  r11™11  where  Q  *Q^  a  0  and  R  ■  RT  >  0. 
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For  reference  purposes  the  state  feedback  linear  quadratic  regu¬ 
lator  problem  is  defined  as: 

r  i 


minimize 


ui4.iiiouit.fc  1  1  00  m  rp 

u  =  -Kx,K€R  \  f  f  (x  Qx  +  u  Ru)dt 
x  =  Ax  +  Bu  o  J 

Defining  the  state  feedback  Ricatti  equation  as 


(2.1.1) 


A^  +  MA.  -MBR"1BTM  +  Q  =  0 


(2.1.2) 


it  is  well  known  that  if  (a.Vq)  is  a  detectable  pair  then  the  minimizing 
control  law  is  given  by  u  = -Kx  where  K = R  "LBTM  and  M  is  the  unique 
symmetric  positive  definite  solution  of  (2.1.2).  If  (A,VQ)  is  not  a  detect¬ 
able  pair  and  none  of  the  eigenvalues  of  the  matrix 


E  = 


A 

-Q 


-1  T 
-BR  B 


-A 


(2.1.3) 


lie  on  the  imaginary  axis,  then  (2.1.2)  has  at  least  two  positive  semi-' 
definite  solutions  M.  More  generally,  if  E  has  p  unobservable  eigenvalues 
with  positive  real  parts  then  there  are  at  least  2p  such  solutions,  in¬ 
cluding  the  optimal  solution  Mo  and  the  unique  stabilizing  solution  Ms,  and 
they  satisfy  [29] 


OSM  SMSMo  (2.1.4) 

o  s 

In  particular,  for  the  minimum  energy  problem  (Q  =  0,  R  =  I)  the  minimizing 
solution  is  Mq *  0, and  if  A  is  unstable  (2.1.2)  has  a  unique  stabilizing 
solution  Mg  with  the  property  of  "reflecting"  unstable  eigenvalues  of  A 
about  the  imaginary  axis. 


The  output  feedback  regulator  problem  is  formulated  under  the 
assumption  that  the  initial  state  Xq  is  a  zero  mean  random  variable  with 
covariance  matrix  Qq  in  order  to  eliminate  the  dependence  of  the  cost 
functional  upon  the  initial  state.  The  problem  is  to 


minimize 

u  “  -Ky,K  €  R1”*1  E 
x  = Ax+Bu 
y  *  Cx 


J  (xTQx  +  uTRu)dt 


(2.1.5) 


Introducing  the  symmetric  positive  definite  matrices  L€Rnxn,  M€RnXn 
defined  by 


L  *  E  (xxT)dt  ,  -|x^x0=  j  Jo(xTQx  +  uTRu)dt 


(2.1.6) 


the  necessary  conditions  for  a  solution  to  the  problem  are  the  coupled 
Ricatti  equations: 


f^+mf+q  +  ctktrkc  =  0 

(2.1.7a) 

FL  +  LFT  +  Q0  =  0 

(2.1.7b) 

K-r”1BTMLCT(CLCT)“1 

(2.1.7c) 

F  -  A  -  BKC 

(2.1.7d) 

Little  is  known  regarding  the  existence  and  properties  of 
solutions  of  these  equations  beyond  the  following  sufficiency  condition. 
Theorem  2.1.1  f 301 :  If  there  exists  an  output  feedback  matrix  K € R™^  such 
that  A  -  BKC  is  a  stable  matrix,  then  there  exists  a  solution  to  equations 
(2.1.7)  for  all  Qq>0»  R>0,  and  QiO  provided  the  pair  (A^VQ)  is 

D 


observable. 
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Thus  if  the  triple  (A,B,C)  may  be  stabilized  by  output  feedback, 
then  the  output  feedback  regulator  problem  has  a  solution.  Furthermore, 
using  the  stabilizing  matrix  KQ  as  an  initial  guess  for  the  solution  of 
(2.1.7),  the  following  numerical  iteration  scheme  has  been  proposed  [25]  . 


FiMi+l +  Mi+lFi  +  Q  +  C^i^i0  =  0  ’  Fi  =  A  -  BKic 


K1+i  -  !l'1Ai+iLi+icI(ni+icI)'1 


Fi+lLi+l  +  Li+lFi+l  +  Q0  0 


(2.1.8a) 

(2.1.8b) 

(2.1.8c) 


Given  K^,  (2.1.8a)  is  solved  for  which  determines  ■  Ki+l^i+l) 

by  (2.1.8b).  Solving  equation  (2.1.8c)  for  Li+^  gives  numerically 

and  completes  the  ith  iteration.  In  practice  this  scheme  frequently 
converges  but  there  is  no  general  convergence  proof. 

For  later  purposes  it  is  noted  here  that  if  an  arbitrary  output 
feedback  K  is  applied  to  the  triple  (A,B,C)  the  associated  cost  is  jtrace(MQQ) 
where  M=M(K)  is  the  solution  of  (2.1.7a).  If  M,  L  is  a  solution  of 
(2.1.7a-d)  then  of  course  this  is  the  optimal  cost. 

In  order  to  obtain  further  insight  into  the  properties  of  the 
output  feedback  regulator  problem,  the  problem  has  recently  been  reformu¬ 
lated  so  as  to  eliminate  the  dependence  of  M  and  L  on  the  covariance  matrix 
Q0  [31].  Noting  that  equation  (2.1.7a)  may  be  rewritten  as 

ATM  +  MA  +  MBR"1BTM  +  Q+W(M,K)  =  0  (2.1.9a) 

W(M,K)  * (r‘1BTM-KC)TR(r’1BTM-KC)  (2.1.9b) 

ie  *  i  * 

it  may  be  shown  that  if  K  minimizes  W(M,K)  then  the  solution  M  =  M  (K  ) 
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of  (2.1.9a)  is  Che  optimal  solution  of  the  output  regulator  problem, 

independent  of  the  distribution  of  initial  states.  However,  a  K*  that 

makes  W(M,K)  minima1,  in  the  positive  semi-definite  sense  exists  only  if 

either  C  is  invertible,  or  if  the  pair  (A,C)  is  completely  aggregable 

(CA  ■  AqC  for  some  Aq)  and  the  weighting  matrix  Q  may  be  decomposed  as 
T 

Q«C  QqC.  Both  of  these  cases  are  equivalent  to  the  state  feedback  reg¬ 
ulator  problem  ,  under  a  transformation  of  basis  in  the  first  case,  and 
under  a  reduction  of  state  in  the  later.  It  was  therefore  suggested  in 
[31]  that  K  be  chosen  to  minimize  the.  term  R  '‘IV^M-KC  in  W(M,K)  with  respect 
to  the  matrix  norm  induced  by  the  inner  product  (x,Ly)  in  order  to  make  the  j 

contribution  of  W(M,K)  in  (2.1.9a)  small.  For  a  given  L  the  minimizing  K  i 

is  given  by  j 

K  =  R"1B1m.CT(CLCT)":L  (2.1.10)  j 

> 

i 

•* 

and  substitution  in  equation  (2.1.9a)  yields  the  equations:  j 

t 

A^+MA  -MBR'1BTM+Q  +  (I-P)TMR_1BTM(I-P)  =  0  (2.1.11a)  < 

i 

P  -  LCT(CLCT)_1C  (2.1.11b)  | 


Thus  for  a  given  L  the  necessary  conditions  for  this  modified  output  feed¬ 
back  regulator  problem  are  the  existence  of  a  positive  definite  matrix  M 
satisfying  equations  (2.1.11).  The  corresponding  feedback  gain  is  given 
by  (2.1.10). 


Relating  this  modified  problem  to  the  output  feedback  regulator 
problem  defined  above  it  is  known  that  if  for  some  (Q,R)  for  which  (A^Vq) 
is  observable  there  exists  an  L>0  such  that  (2.1.11)  has  a  positive 


1 
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definite  solution  M,  then  for  any  Q  i  0,  R>0,  Qq>0,  (A^VQ)  observable,  the 
necessary  conditions  (2.1.7)  for  the  output  feedback  problem  have  a  solution 
L  >0,  M  >0  '[31]  . 

The  properties  of  solutions  of  the  modified  regulator  problem  will 
be  discussed  in  the  next  section. 

2.2.  Retention  of  Optimal  Invariant  Subspaces  by  Static  Output  Feedback 
Compensation 

Since  the  solution  of  the  output  feedback  regulator  problem  does 
not  have  an  analytic  characterization  it  is  of  interest  to  obtain  sub- 
optimal  output  regulators  associated  with  the  state  regulator.  This  section 
gives  a  solution  to  the  problem  of  determining  output  feedback  gains  which 
assign  X  dimensional  invariant  subspaces  (l£X^r)  of  the  optimal  state 
regulator.  Consideration  is  restricted  to  those  output  feedbacks  which 
may  be  obtained  as  "generalized  projections"  of  state  feedbacks: 

K0  =  KsP1  +  P2  (2.2.1) 

-1  T 

Let  Kg  = R  B  Mc  be  the  solution  to  the  state  feedback  regulator 
problem  for  a  given  Q,R,where  Mc  is  the  solution  of  the  Ricatti  equation, 
and  define  the  two  problems: 

a)  For  1  s  X  <  r  determine  an  output  feedback  gain  Kq  such  that  a  pre¬ 
scribed  X  dimensional  invariant  subspace  of  (A-BKg)  is  also  an  invariant 
subspace  of  (A-BKqC). 

b)  For  X  -  r  determine  an  output  feedback  gain  Kq  such  that  an  r  dimensional 
invariant  subspace  of  A-BKg  is  also  an  invariant  subspace  of  A-BKqC  and  Kq  is 
optimal  with  respect  to  the  modified  output  feedback  regulator  problem 
define  above. 


L 
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The  first  problem  has  a  non-unique  solution  and  may  be  solved 
using  generalized  inverses.  This  gives  rise  to  an  output  feedback  pole- 
placement  problem  which  may  be  used  to  shape  the  complementary  spectrum 
of  A-BKqC.  Thus  Kq  may  be  chosen  to  retain  an  JL  dimensional  subspace 
which  is  optimal  with  respect  to  a  state  feedback  regulator  and  the  re¬ 
maining  freedom  in  Kq  used  to  shape  the  spectrum  of  the  closed  loop  system. 

Let  {u^J/^,  uk  6  RnX'1'  and  ^  be  the  eigenvectors 

and  eigenvalues  of  the  optimal  closed  closed  loop  system  F=A-BKs,  where 
the  first  JL  eigenvectors  span  the  JL  dimensional  invariant  subspace  of  F  to 
be  assigned  to  A-BKQC.  It  is  assumed  that  (XX,X is  not  a  complex  pair. 
In  order  to  work  over  the  reals  define  a  transformation  Tr  by 


[row  (T  )]  =  (0,.. .,0,1,0, ...,0)  if  is  real 

t  ith  position 


(2.2.2) 


r°w1 (Tr) 
rowi+l<Tr> 


0>  •  •  •  > 0>2 > “2". 3 >0) 

J^O, . . .  ,0,— ,jj ,  0,...,0J 

f  ith  position 


if  are  a 

complex  pair 


Then  under  Tf  the  complex  pairs 


mapped  to 


Defining 


C  +  ju>  0 
0  a  -  juj 


and  [u  +  jv  u  -  jv]  are 


and  [u  v]  respectively. 


Ux  »  <uL  j  ...  I  u^T.  ,  hl  =  T^1  dg(X1,...,Xi)Tr  (2.2.3) 

the  problem  may  be  stated  as  that  of  finding  a  gain  Kq  ant’  a  matrix  P 
such  that 


(A-BKoC)Ux  =  UXAX 
(A-BKgP)Ux  -  U^A JL 
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or  equivalently,  assuming  B  is  of  full  rank,  such  that 


haJt  ■  KA 


K  P  3  K  C 
s  o 


(2.2.5) 


The  general  solution  of  equations  (2.2.5)  is 


K  =  K  U  [CUi]1+X[I  -CUX(CU  J1],  X  6  RmXr 


P  =  K  K  C 
s  o 


(2.2.6) 


1  ^ 

where  A  denotes  any  1-inverse  of  A  [32].  By  assumption  (A,C)  is  observable 

and  if  CU^  has  rank  i<  r  one  choice  for  (CU^)L  is  [(CU,)1^^)]”1^!^)1. 

It  should  be  noted  that  (CU jj)^ (CV ^)  *1^.  If  also  (A,Vq)  is  observ- 
-1  T 

able  then  K  =  R  B  M  is  of  full  rank  since  M  >0,  and  the  1-inverse  of  K 
s  c  c  s 

1  T  t  -1  1 

may  be  taken  to  be  K  =K  (K  PC )  ,  which  satisfies  K  K  = I  .  However,  if 

J  ss  ss  ssm 

(A^Q)  is  not  observable,  and  in  particular  in  the  case  of  the  minimum 


energy  problem,  M  and  K  need  not  be  of  full  rank  and  the  formula  for  K  is 
c  s  s 


replaced  by 


K_1 

1  KH 

Ks  =  \ 
s  0 


ZK  Y  = 
s 


hi  hi 


(2.2.7) 


where  K^^6R^X^,  Y  €  R™™,  Z€RnXn,  and  rank  (Kg)=p. 

Since  rank  (Ir*CU^(CU^)  ^)-r-X  there  are  only  m(r-i)  degrees  of 
freedom  in  the  matrix  X  in  (2.2.6).  To  eliminate  the  redundancy  let 


X  is  a  1-inverse  of  A  if  AXA=A. 


by  X  =  YS.  Also  note 


S€R^r”'^xr  be  of  full  rank  and  define  Y€Rmx^r  ^ 

that  since  P  is  to  be  eventually  multiplied  on  the  left  by  Ks,  the  factor 

K^K  in  (2.2.6)  may  be  omitted.  The  expression  for  P  then  becomes 
s  s 

P  =  U^CU^) LC  +  KgYS (Ir  -  CUX(CU/)C  (2.2.8) 

Defining  the  matrices : 

*0  -A-mt,«i(c0Jt)1c 

Bn  =  -BK  K1  (=  -Bif  K  is  of  full  rank)  (2.2.9) 

u  s  s  s 

cQ  =  s(ir  -cu  ^cu^Sc 

and  substituting  (2.2.8)  into  (A-BK  P)  gives  the  output  feedback  pole- 

s 

placement  problem  of  finding  Y  gR™*^  ^such  that  the  spectrum  of  Aq  +  BqYCq 
is  satisfactory. 

Since  by  construction  Aq  +  BqYCq  contains  an  1  dimensional 
invariant  subspace  it  is  possible  to  exploit  this  fact  to  reduce  the  di¬ 
mensionality  of  the  problem.  Let  T=[T^  T^]  be  any  invertible  matrix 

with  and  T2  ^  R',n  J  •  Applying  the  transformation  T  to  the  triple 

(A0,B0,C0)  gives: 

,  CqT*[0  I  C1]  (2.2.10) 

,  .  _  _  (n-^)x(n-X)  _  -_(n-X)xm  „  „  „rx(n-X)  .  , 

where  A^  6  Rv  ,  B^6RV  ,  6  R  .  The  pole-placement 

problem  to  be  solved  is  then  that  of  satisfactorily  shaping  the  spectrum  of 

A^  +  B^YC^  by  an  appropriate  Y  €  R^^  ^  .  The  final  output  feedback  gain  will 


t'\t  = 


Ax  i  x 


1  A 

1  Ai 


X 

B, 


be  given  by 


K0  =  KsU,((iU<e)L  +  YSlIr-CUi(CU,)1] 


(2.2.11) 


and  the  corresponding  projection  matrix  P  by  (2.2.8). 

In  summary  equation  (2.2.11)  characterizes  the  m(r-X)  degrees  of 
freedom  available  in  solving  the  problem  of  retaining  a  prescribed  i 
dimension  invariant  subspace  of  the  optimal  state  regulator.  The  associated 
pole -placement  problem  may  be  solved  by  any  pole-placement  procedure,  but 


in  particular  the  dyadic  solution  for  which  software  support  has  been 
developed  is  appropriate. 

If  it  is  desired  to  retain  an  r  dimensional  invariant  subspace  of 
the  state  feedback  regulator,  it  is  possible  to  choose  KQ,P  so  as  to 
additionally  solve  the  modified  output  feedback  regulator  problem  defined 
in  the  previous  section.  Whereas  in  the  case  X<r  it  is  only  guaranteed 
that  using  output  feedback  achieves  the  optimal  state  regulator  cost  in  an 
l  dimensional  subspace,  in  the  case  X=r  the  output  feedback  control  law 
will  also  be  optimal  in  the  sense  of  the  modified  regulator  problem. 

In  order  to  simplify  the  calculations  it  is  assumed  that  a  trans¬ 


formation  has  been  applied  to  the  triple  (A,B,C)  such  that  C  =  [Ij.  0]. 

Denote  by  Mc  the  solution  of  the  state  feedback  Ricatti  equation  and  by 
M(L)  the  solution  of  the  modified  regulator  problem.  Let  anc* 

be  r  eigenvectors  and  eigenvalues  of  (A-BR  S'S^)  and  introduce  the 
partitions : 


lu»  11’  nr 


Y  6  R 


rxr 


Kl  €R 


mxr 


defined  in  (2.2.2) 


Then  the  following  holds  [31] . 

Theorem  2.2.1:  If  the  matrix  Ar  *  (A.^-ZY  ^A^)  has  all  its  eigenvalues  in 
the  left  half  plane  then: 

1)  The  modified  output  feedback  regulator  problem  has  a  solution  M(L)=M(N) 
for  each  L  satisfying  I^l^ll  =  ZY  ^  =  N>  and  the  associated  feedback 
gain  is  given  by 


k*r"1btm(n)p  =  [Kj^+^zy'1  :  0] 


(2.2.13) 


2)  The  spectrum  of  the  closed  loop  output  feedback  system  A-BK  is  given 

by  tX.^}  U  a(Ar) . 

3)  The  cost  matrix  M(N)  may  be  decomposed  as  M(N)  =Mc+D(N),  D(N)s0j 
where  D(N)  represents  the  cost  increase  over  the  optimal  state  feedback 
solution,  and  the  null  space  of  D  is  spanned  by  {u.}  f ^  so  that  in 
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this  r  dimensional  retained  subspace  there  is  no  increase  in  cost 
associated  with  using  the  output  feedback  control.  D  is  given  by 


D 


022  6E(n-r)lI<n-t> 


(2.2.14) 


T 

where  D22=D22>^  so^-ut*-on  t*ie  LyaPunov  equation 

(A22  -NA12)TD22+D22(A22  -N412)  +  r22  -  0  (2.2. 15) 2 

Furthermore  the  matrix  D(N)  provides  a  bound  on  the  optimal  cost 

for  the  output  feedback  regulator  problem.  Defining  J  as  the  cost  associated 

s 

with  the  optimal  state  feedback  regulator  and  Jq  as  the  cost  associated  with 
the  optimal  output  feedback  regulator  then  [31] 

J  £  J  £  J  +  |trace(DQ  )  (2.2.16) 

s  o  s  2  s  o 

In  summary  the  problem  of  determining  a  gain  to  retain  r  optimal 
eigenvectors  has  the  unique  solution  K^K^  +  I^N.  It  is  noted  that  this  is 
precisely  equation  (2.2.11)  under  the  constraints  X  =  r  and  C  ~ [1^  0], 


2.3.  Retention  of  Optimal  Invariant  Subspaces  by  Dynamic  Output  Feedback 
Compensation 

In  the  event  that  it  is  not  possible  to  stabilize  the  system  by 
static  output  feedback  while  retaining  a  desired  X  dimensional  subspace  of  the 
state  feedback  regulator, or  that  the  spectrum  of  the  resultant  system  is  not 
acceptable,  a  dynamic  compensator  of  dimension  p  may  be  designed  which  will 


retain  an  r  +  p  dimensional  subspace  of  the  optimal  state  feedback  regulator. 


In  this  section  the  derivation  given  in  [1]  of  a  design  oriented  approach 
to  the  construction  of  such  a  compensator  of  prespecified  dimension  p  will 
be  presented.  This  approach  reduces  the  specification  of  the  parameters 
of  the  compensator  to  the  solution  of  an  output  feedback  pole-placement 
problem  similiar  to  that  encountered  in  the  previous  section  in  the  design 
of  static  output  feedback  compensators.  In  the  next  chapter  an  algorithm 
will  be  given  which  simultaneously  solves  this  pole-placement  problem  and 
determines  without  apriori  assumptions  the  dimension  of  the  desired  compen¬ 
sator. 

Again  for  simplicity  the  output  matrix  is  assumed  to  be  in  the 

form  C*[I  0].  Introducing  the  compensator  z=Hz+Dy,  H  6  R^x^,  D  6  Rpxr 
r 

and  momentarily  assuming  that  the  matrices  H  and  D  and  the  dimension  p  are 
known,  the  compensator  design  problem  may  be  treated  as  a  pole  placement 
problem  to  which  Theorem  2.2.1  may  be  applied.  Defining  the  matrices: 


back  u  = -K  C  x  which  retains  an  r  +  p  dimensional  subspace  of  the  optimal 
state  feedback  regulator  associated  with  (A,B,Q,R)  and  which  is  optimal  in 
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the  sense  of  the  modified  regulator  problem. 

Since  the  compensator  state  is  not  observable  throughV^T  the 
solution  of  the  state  feedback  regulator  problem  is 


M  = 
c 


M 


,  K  =  r"1^  =  [0  :  R_1BTM  ] 


(2.3.2) 


where  is  the  solution  of  the  Ricatti  equation  associated  with  (A,B,Q,R). 
Note  however  that  if  H  is  not  stable,  then  (A^Y(T )  is  not  a  detectable  pair 
and  though  optimal,  the  solution  (2.3.2)  will  not  be  stabilizing.  (See 
remarks  preceding  equation  (2.1.4).)  This  poses  no  difficulty  in  the  design 
methodology  as  the  optimal  closed  loop  system  will  then  consist  of  the  system 
x  =  (A-BR  '*'BTMc)x  driving  an  unstable  open  loop  compensator.  Even  though 
the  compensator  states  will  diverge,  the  plant  response  will  be  stable.  The 
implemented  output  feedback  control  law  will  stabilize  the  total  closed  loop 
system. 

Let  and  ei-8envectors  an<*  eigenvalues  of 

the  optimal  state  feedback  regulator,  where  the  first  r  +  p  eigenvectors  span 
the  r  +  p  dimensional  subspace  to  be  retained  and  (^r>^r+i)  and  ^r+p’^r+p+P 
are  not  complex  pairs.  Denote: 


(u  :...: 

■  111 

(ur+l 

.:  u  ^  )t  - 
.  r+p'  r  [ 

-1 

A  -  T 
r  r 

dg(Xl  ...  Xf)T 

W  €  RPXP 

and  W  €RpXr 

defined  in  (2.3.2) 


(2.3.3) 


-1 
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(2.3.4) 


By  Theorem  2.2.1,  stabilization  of  the  system  by  output  feedback  of  the 
measured  variables  and  the  compensator  states  requires  the  selection  of 
H,D  and  therefore  and  such  that  the  matrix  A^  has  a  satisfactory 
spectrum  where: 


N  =  [Z 


V] 


W  W 

r  P 

Y  U 


-1  &  [N  Nr],  N  6R(n-r)xP,Nr€R(n-r)xr 


=  A22  -  NrA12 


(2.3.5) 


Using  the  formula  for  the  inverse  of  a  partitioned  matrix  the  expression 
for  Ar  may  be  written  A^  +  BqEA^  where: 

AL  -  Aj2  -  ZY*1^ 

Bq  =  V-ZY_1U  6R(n‘r)xp  (2.3.6) 

P  =  L(Y-UL)"1  €Rpxr  ,  L  -  w“1Wr  €  RPXr 

Thus  the  satisfaction  of  the  condition  of  Theorem  2.2.1  that  A^  have  an 
acceptable  spectrum  is  reduced  to  the  solution  of  an  output  feedback  pole- 
placement  problem.  In  the  next  chapter  this  problem  will  be  solved  by  an 
algorithm  which  computes  P  as  the  sum  of  a  sequence  of  dyadic  feedbacks  and 
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determines  the  compensator  dimension  p.  It  is  noted  here  that  the  number 
of  columns  of  Bq  is  equal  to  the  compensator  dimension  and  that  rank  (Bq) 
is  always  maximal  since  the  columns  of  Bq  are  eigenvectors  of  a  particular 
matrix  arising  in  the  block  triangularization  of  F  [1], 

Assuming  that  this  pole-placement  problem  has  been  solved  and  that 
the  dimension  p  and  the  gain  P  are  known,  the  parameters  of  the  compensator 
are  given  by: 


h  *  vy'p1  •  Ho  ■ 

d*updo-  -  ovV-iy-^+up] 

L  =  (I  +  PU) ~  1PY 


N  =  N  W-1  ,  N  *  (V  -  ZY_1U) (I  +  PU) 
P  PO  P  Po  P 


Nr  =  ZY-1  -  (V  -  ZY_lU)P 


(2.3.7) 


A 


[K2  :  Ky  :  0] 


Thus  the  matrix  P  determines  the  compensator  up  to  a  similarity  transforma¬ 
tion  Wp  which  may  be  used  to  obtain  a  favorable  representation  of  the  pair 
(H,D).  The  final  closed  loop  spectrum  consists  of  the  r  +  p  retained  optimal 
eigenvalues  together  with  the  spectrum  of  A^. 

Before  considering  the  solution  of  the  pole-placement  problem  in 
the  next  chapter  some  remarks  regarding  the  resultant  compensator  dimension 


p  are  appropriate. 


By  Theorem  1.3,  if  psn-2r+l  then  the  spectrum  of  A^  may  almost 
always  be  assigned  arbitrarily.  Thus  a  bound  on  the  dimension  of  compensator 
required  to  satisfactorily  control  the  system  is  n-2r+l.  In  the  case  this 
bound  is  achieved,  the  resultant  closed  loop  system  will  retain  r  +  p=n-r  +  l 
optimal  eigenvectors  and  will  admit  an  arbitrary  complementary  spectrum 
through  Ar>  As  will  be  seen  in  numerical  examples  in  Chapter  4,  acceptable 
designs  may  be  obtained  with  compensators  of  dimension  well  below  this 
bound,  particularly  when  the  goal  is  to  place  all  the  eigenvalues  in  a  pre¬ 
scribed  region  of  the  complex  plane  rather  than  at  prescribed  locations. 

In  the  case  p=n-r  all  n  optimal  eigenvectors  will  be  retained 
and  the  resultant  compensator  may  be  identified  with  the  reduced  order 


Luenberger  observer.  To  make  the  correspondence  explicit  consider  an  observer 


given  by  [33] 


w  *  Ew  +  Gy  +  Ru  ,  E  €  R(n'r)x(n‘r) ,  G  €  R(n"r)xr,  R  6  R(n-r)xm 


x,  =  w  +  Sy  ,  S  6  R 


(n-r)xr 


(2.3.8) 


where  x^  is  an  estimate  of  the  unmeasured  state  variables  x^.  (Recall 
C»[I  0]).  Defining  an  error  e  =  x-x  and  introducing  the  partitions: 


An  A12 


Si  A2 


F11  F12 


F21  F22 


,  A.-JR™  B.JR™ 


(2.3.9) 


-1  T 

,  F  -  A  -  BR  B  M 

9  n 
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it  may  be  shown  that 

e  -  Ee  +  (A22  -  SA12  -  E)x2  +  (A^  -  SAU  -  G+ES)xL 

+  (B2-SB1-R)u  (2.3.10) 

Thus  the  error  will  converge  asymtopically  to  zero  independent  of  x  and  u 
provided 


E  -  A22  -  SA12 


G  “  A21  ‘  SA11  +A22S  "  SA12S 


R  «  B2  -  SB1 


(2.3.11) 


and  E  is  a  stable  matrix.  Observability  of  (A,C)  implies  that  the  pair 
(A22>A^2)  is  observable  and  so  the  pole-placement  problem  for  E  has  a 
solution. 


Implementing  the  control  law  as 
/x. 


u  *  r 


(K1  Kj) 


r  -  (I^  :  Kj+KjS) 


w 


(Kj  K2)  -  R-  13TMc 


(2.3.12) 


and  using  (2.3.11)  the  total  closed  loop  system  becomes 


d 

dt 

w 

X1 

sa 

F22‘SF12  F2l"SFll  +  f22S'SF12S  0 

“Bj^Kj  Au  -B1(Kl+K2S)  A12 

w 

X1 

_X2_ 

-B^  A21-B2(K1+K2S)  A22 

_X2_ 

39 


B2  “  SB1 


(2. 3-.  13) 


By  comparison,  the  compensator  of  this  section  is  of  the  form 
v  *  Hv  +  Dy  and  may  be  viewed  as  producing  an  estimate  =  N  v  +  N^y. 
Under  the  control  law 


u  =  r-(0  KL  Kj)  P  xx  *  r  -  (K^  l  I^+K^  :  0) 


(2.3.14) 


the  closed  loop  system  is 


dt  X1  *  “BlK2Np  All'Bl(Kl+  *2**^  A12  x!  +  B1  r 


W,  A2rVKl+K2Nr>  A22  1 1  X2  B2 


(2.3.15) 


Provided  N^1  exists,  then  under  the  transformation  (v.x^.x^  -  (n”1v,x1,x2) 

(2.3.15)  becomes 


N  HN_1 
P  P 


X1  “  ^i^  A11 'Bi(Kl+ W  A12  X1  +  B1  r 

*2  -B2K2  A2l  "  B2 (KL  +K2Nr)  A22  X2  B2 


Now  in  [1]  it  is  shown  that  the  transformation  T 


&■:] 


(2.3.16) 


upper 


block  triangularizes  the  optimal  closed  loop  state  feedback  matrix  provided 


N-ZY  exists  and  is  a  matrix  of  optimal  eigenvectors  of  the  closed 


l 


loop  system.  Applying  this  result  to  the  compensated  system  of  this  section 
gives 


H 

D 

0 

H 

D 

0 

T 

0 

F11 

F12 

t”1  = 

F10N 

12  p 

Fll+F12Nr 

F12 

0 

F21 

F22 

f (N  ,N  ) 
p  r 

g(N  ,N  ) 
p  r 

F22  '  NrF12 

(2.3.17a) 


where 


T 


0 


I 


r 


0 


-N 


r 


I 


n-r 


(2.3.17b) 


f<yv  -  <-y<+F22-NtF12)Np  ,0 

(2.3.17c) 

g(Np,Nr)  =  <-V+F21-NrFll  +  F22Nr-NrF12V  *° 

Identifying  S~Nr  and  comparing  (2.3.13)  and  (2.3.16)  it  follows 
that  the  compensator  of  this  section  is  a  Luenberger  reduced  order  observer 
in  an  appropriate  basis  provided  the  reference  input  is  zero.  Furthermore 
the  dynamics  of  the  error  equation  are  governed  by  the  matrix- 


Ar  -  Aj2  "  Nr^l2  *  +  B0FA12  (2.3.18) 

whose  eigenvalues  may  be  arbitrarily  assigned  by  state  feedback  pole-placement 
since  BQ  is  square  and  of  full  rank  n-r. 


ud 
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CHAPTER  3 

FURTHER  ASPECTS  OF  THE  DESIGN  METHODOLOGY 

3.1.  An  Algorithm  Solving  the  Associated  Pole-Placement  Problem 

This  section  presents  an  algorithm  which  solves  the  output  feed¬ 
back  pole-placement  problem  of  section  2.3,  associated  with  the  design  of 
compensators  for  suboptimal  regulators.  The  problem  is  to  determine  the 
compensator  dimension  p  and  a  matrix  P € R^Xr  such  that  the  spectrum  of  Ar * 

A^  +  BqPA^  i-s  satisfactory  (see  equation  (2.3.6).  A  solution  is  obtained 
here  by  iterating  the  proof  of  Theorem  1.3  in  the  manner  suggested  in  [2], 
Recalling  that  the  number  of  columns  of  Bq  is  the  dimension  of  the  compensa¬ 
tor,  the  idea  is  to  start  with  BQ  a  column  vector  and  successively  increase 
the  column  span  of  Bq,  solving  a  partial  pole -placement  at  each  stage,  until 
satisfactory  spectral  characteristics  are  obtained  for  A^,.  The  degree  of 
the  required  compensator  will  then  be  the  number  of  columns  of  Bq  and  the 

parameters  H,D,K  ,K  of  the  design  will  be  given  in  function  of  P  by  (2.3.7). 
z  y 

The  algorithm  consists  of  the  solution  of  a  sequence  of  pole- 
placement  problems  defined  as  follows.  Given  from  the  first  part  of  the 
design  are  the  observable  pair  (A22>A^2^  an<*  t*ie  °Ptimal  closed  loop  eigen¬ 
values  and  eigenvectors  where  the  are  distinct.  It 

is  assumed  that  the  X,  and  v.  have  been  ordered  such  that  (X  ,X  , ,)  is  not  a 

i  i  ' r’  r+1' 

complex  pair  and  such  that  the  matrix  Y  defined  in  (2.3.3)  is  invertible. 

The  ordering  may  aLso  have  taken  into  account  the  relative  importance  which 
the  designer  may  attach  to  the  optimal  eigenvalues.  Define: 


y  €r 


rxr 


Z  €  R(n_r)xr, 


T  defined  in 

(2.2.2) 


i 


n-r 


i  v 

I  n-r 


=  (vt+1  u(€Rrl11,  ..En'"-'*11 


n  r  i 


(3.1.1) 


bj.=vi-N0ui,  Nq=ZY  ,  Bi  “  (bi  ••••  bi)»  i  =  l,...,n-r 


Since  A^2  bas  tbe  role  of  an  output  matrix  and  since  it  will  be 

required  for  the  output  matrix  to  have  full  rank  in  the  pole- placement  pro- 

Xxr 

cedures  to  be  described,  denote  X  =  rank(A^2),  let  T€R  be  any  matrix  of 
full  rank,  and  define  a  new  output  matrix  C=TA^2*  Also  since  the  pole- 
placement  problem  for  the  triple  (Aj,Bq,Aj2)  may  be  solved  as  a  state  feed¬ 
back  pole-placement  problem  if  Xa n-r,  this  case  will  be  discussed  at  the 
end  of  the  section.  Presently  it  will  be  assumed  that  X  +  r<n. 

The  algorithm  may  be  stated  in  terms  of  two  pole-placement  proce¬ 
dures,  to  be  defined  below,  as  follows. 

Algorithm 

0.  Initialize  i = 0. 

1.  Let  i = i+1. 

2.  Using  either  procedure  1  or  2  below,  compute  a  dyadic  feedback 

fl8i,  f ^  €  Rixl,  gi€RlxX  (3.1.2) 


to  assign  i  +  X-1  desired  eigenvalues  to  the  matrix: 


Ai+l“Ai  +  BifigiC 


(3.1.3) 


i+1 


3.  If  the  resultant  spectrum  of  A 


is  unsatisfactory  or  if  the  algorithm 
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may  not  terminate  because  X^,  X^+^  a  comPlex  pair,  8°  to  1. 

4.  Let  the  compensator  dimension  be  p  * i  and  let 


P  = 


i 

£ 

j=l 


f . 

..1. 


gj  T  *  Ar=Ai+l 


(3.1.4)  a 


Since  the  spectrum  of  A  may  be  assigned  arbitrarily  if  p+2-l  =  n-r  this 
algorithm  will  terminate  after  at  most  n-r-i+1  stages. 

Procedures  1  and  2  apply  Theorem  1.3  to  the  systems  (A^,B  ,C)  de¬ 
fined  by  the  algorithm.  At  the  ith  stage  the  system  (A^jB^.C)  has  i  inputs 
and  L  outputs.  Thus,  with  possible  exceptions  due  to  the  system  lying  on  a 
hypersurface  on  which  the  results  of  the  theorem  fail  to  hold,  i  +  X-1 
eigenvalues  of  Ai+1  may  be  assigned  arbitrarily  closely  to  desired  values. 
The  procedure  then  is  to  choose  an  input  (or  output)  space  projection  which 
reduces  the  system  to  a  single-input  (or  single -output)  system  and  renders 
i  -  1  eigenvalues  uncontrollable  (or  i- 1  eigenvalues  unobservable).  These 
eigenvalues  correspond  to  eigenvalues  already  assigned  during  the  (i-l)st 
stage.  Next  feedback  gains  are  computed  to  place  at  desired  locations  as 
many  additional  eigenvalues  as  there  are  outputs  (or  inputs).  Because 
uncontrollable  and/or  unobservable  eigenvalues  are  invariant  under  output 
feedback  this  will  result  in  the  assignment  of  i  +  X-1  eigenvalues  of  the 
matrix  A^+^.  Both  state  space  and  frequency  domain  based  procedures  will  be 
given. 

It  will  be  assumed  that  (A^,b^)  is  a  controllable  pair,  and  there¬ 
fore  that  each  pair  (A^»B^)  *-s  controllable.  The  case  when  (A^,b^)  is 


uncontrollable  will  be  discussed  at  the  end  of  the  section. 


Procedure  1 


The  system  (A^,B^,C)  is  transformed  to  the  single-input  system 


,ixl 


(A^jB.f  ,C)  where  f^GR  is  chosen  to  render  i-1  eigenvalues  of  A^ 

,lxi 


uncontrollable.  A  vector  g^^  6 R  is  then  found  to  place  l  eigenvalues  of 
Ai+^ = A^ + B^f If  i  = 1,  then  f ^ =  1,  and  the  first  part  of  the  procedure 
is  omitted. 

State  Space  Solution 

Let  i  be  the  eigenvalues  of  A^,  the  first  i-1  of  which  are 

to  be  retained,  and  let  {  vkGRlx^n  r\  be  the  corresponding  left 

eigenvectors.  With  Tr  given  by  (2.2.2)  define 


V  =  T 

VI  lv 


V1 

6R(i-l)x(n-r)>  =t^T 

Vi 

» 

♦ 

— r  “ 

_vi-l_ 

V 

n-r 

€R(n-r-i+l)x(n-r) 

(3.1.5) 


Recalling  that  an  eigenvalue  is  uncontrollable  if  the  corresponding  left 
eigenvector  is  in  the  null  space  of  the  input  matrix,  select  a  vector 


f.^  €  Rix^  satisfying  V^Bj^  =  0.  Then  [sk}k_^  will  be  invariant  under  further 


feedback  and  for  each  k  such  that  the  kth  entry  of  nonzero>  the 

eigenvalue  will  be  controllable.  A  solution  of  V^B^f^  *  0  always 

exists  since  this  corresponds  to  finding  an  i-vector  orthogonal  to  (i-1) 
i-vectors.  However  conditions  under  which  the  n-r-i+1  remaining  eigenvalues 
will  be  controllable  are  not  available. 


Equations  (1.5)  and  (1.6)  determine  the  gain  g^  to  place  l  eigen¬ 


values  of  the  triple  (A^.B^f,  ,  c)  flt  desired  locations  (If  § 

entries  of  V2B^f^  are  nonzero  with  §<X  then  §  should  be  used  in  place  of 
l  in  the  following  formulas.)  Let  the  characteristic  polynomial  of  A,  be 
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n-r 


i 

P  (X)  * ,  a,  X  ,  a  =1 
Fov  '  k=0  k  *  n-r 


(3.1.6) 


and  define  the  matrices: 


Q  -  IB.fi 

A. B. f »  ••••• 

ill. 

A?~r 

1 

-1Bifi]  6R(n-r>x(n- 

— 

— 

1  a  . 

a0  an 

n-r-1. 

2  1 

1. 

•  • 

«  * 

a2 

• 

•  • 

•  • 

• 

• 

• 

• 

• 

•  • 

€R(n-r)x(n-r) 

0 

« 

•  • 

• 

o 

• 

1 

— 

1  _ 

n-r-1  n-r-1 

CT1  jfc 

n-r  n-r 

al  . ai 

•  • 

•  • 

•  • 

•  • 

Tr  5R<rl't)xi, 

A 

s  = 

•  • 

•  • 

a ,  a. 

°1  .  i 

1  l 

1  .  1 

1  .  1 

(3.1.7) 


Tr€R 


(n-r+l)xA 


A  A  IVX 

P  -  (P0(a1),...,p0(ajl))Tr=(l,an_r_1,...,a1>a0)S6R 
where  Tr  is  given  in  (2.2.2).  Then  the  solution  of 

g^(CQRS)  -  P  (3.1.8) 

will  assign  to  the  spectrum  of  *  A^  +  B.^g^C  the  ^  +  i"l  eigenvalues 

^ak^k*l  and  ^k^k=l  ’ 

Frequency  Domain  Solution 

The  transfer  function  of  the  system  (A^,B^,C)  may  be  written 
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C(XI  -Ai)'1Bi*Na)  /  d(X) 
n-r 

d(X)  =  £  a,  XK,  a  =1 
k=0  k  n-r 


(3.1.9) 


where  the 

roots  o 

f  d(^)  are 

^^k*!  tbe  first  i"l  are 

to  be 

retained. 

The  numerator  polynomial  matrix  may  be  written 

N(X) 

*n(X) 

• 

• 

• • •  (^) 

• 

• 

n-r-l  , 

=  k5>  Nk  X 

(3.1 

nXl(X) 

*  *  *  nXi 

n-r-1 


nij^)=k=0  nijkX  ’  \  = 


‘Ilk 


lik 


nXlk  ‘ ‘ *  niik 


(In  n, the  index  i  indicates  the  row  of  N(X),  j  the  column  of  N(X),  and  k 

tJK 

the  degree  of  the  term  of  n. .  (X)  in  which  n. appears.) 

Recall  that  the  goal  is  to  render  i-1  eigenvalues  uncontrollable  by 
a  choice  of  f^  and  to  choose  g^  to  place  &  eigenvalues.  Since  an  eigenvalue 
of  a  single-input  system  is  uncontrollable  if  the  numerator  of  the  transfer 
function  is  zero  evaluated  at  that  eigenvalue,  it  follows  that  for  k=^ 
to  be  uncontrollable  eigenvalues  of  the  single-input  system  (A^B^f^  }c  )  it 
must  hold  that  N(ek)f^  =  0  for  k  =  l,...,i-l.  By  Lemma  1.2  rank  N(efc)  =  1. 

Let  ljki£=^  be  any  sequence  such  that  the  j£h  row  of  N(ek)  is  nonzero,  and 
define 


M  = 


e  N(«.) 
J1 


*1 

J  (1-D 


d (i~l)xi 
R  ,  e, 


=  (0, 


0,1,0,. ..,0)  (3.1.11) 

t  kCtl  position 
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Then  may  be  computed  as  the  solution  of 

T  ^Mf.  =0 
r  1. 


(3.1.12) 


where  Tr  is  given  in  (2.2.2).  For  purposes  of  numerical  evaluation  the 
rows  of  the  matrix  Tr^M  may  be  obtained  as 


where 


(1 


e 

k 


e 


n-r- 

k 


if 


1  Re(ek)  , 
0  lm(ek),. 


_  .n-r-1 

Im(*  k) 


N . 


e,  is  real 
k 

(3.1.13) 

if  k,e  k+1^  iS  a  comPlex 

pair 


"j- 


n 


n 


V° 

i 

j l(n-r-l) 


jiO 


nji(n"r"1) 


€R(n'r)xi(J.l.l4) 


It  should  be  noted  that  the  same  difficulties  regarding  the  choice  of  f^ 
are  present  here  as  in  the  state  space  solution.  That  is,  M  may  be  rank 
deficient  in  which  case  there  will  be  a  multiplicity  of  choices  for  f^,  and 
not  all  the  remaining  eigenvalues  k=^  need  be  controllable  with  respect 
to  the  pair  )•  Their  controllability  may  be  verified  by  computing 

the  vectors  N(«k)f^. 

i 

To  assign  the  eigenvalues  {_  ak3k=i  the  8a^n  8^  Is  computed  from 
equation  (1.31)-(1.35).  Since  the  transfer  function  for  the  single  input 
system  (A^,B^f^,  C )  is 

C(U-Ai)"1Bifi  =  V  t  ^  /  d(X) 


(3.1.15) 


Y  —  (3g  a^  •••  “  (0  •••  0  dg  d^  ... 

The  remaining  spectrum  of  is  determined  as  the  roots  of  the  polynomial 

P2(^)- 

If  §  eigenvalues  of  (A^,B^f^  )  are  controllable  with  g<X  then  in 
general  (3.1.17)  will  be  inconsistent  and  must  be  modified  by  replacing  1 
with  g. 

To  illustrate  the  notation  consider  an  example. 

Example  3.a.l:  Suppose  the  system  (A^,B^,C)  at  the  ith  stage  of  the  algorithm 
is  given  by 
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T 

This  has  solution  f2*(3,l)  and  the  resulting  transfer  function  for  the 
single  input  system  Is 


(3.1.21) 
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If  the  remaining  poles  are  to  be  placed  at  -1  +  jl,  -3,  the  desired  partial 

3  2 

'characteristic  equation  is  p^(X)*X  +  5X  +  8x  +  6.  Equation  (3.1.15) 
becomes 


[g2  :  h] 


2-533 


0  2  10 


0  4  4  1 


6  8  5  1 


=  (0  0  -2  1)  -  (0  :  6  8  5) 


(3.1.22) 


and  has  solution 


=  (-6,  -50,  12),  hQ  =  2 .  The  eigenvalues  of  may  be 


verified  to  be  -1+lj,  -3  together  with  the  root  -2  of  X  +  hQ.  □ 

Procedure  2 

The  procedure  is  the  dual  of  Procedure  1.  The  idea  is  to  transform 

lxX 

the  system  (A^.B^C)  to  the  single-output  system  (A^B^g^c)  where  g^Ql 

XX 1 

is  chosen  to  render  X-l  eigenvalues  of  A^  unobservable.  A  vector  f ^  6  R 
is  then  found  to  place  i  eigenvalues  of  A^+^  =A^  +  B^f^gX).  If  X  =  1  then 
the  first  part  of  the  procedure  is  omitted. 

State  Space  Solution 

Let  be  th®  eigenvalues  of  A^,  the  first  X  -  1  of  which  are 

to  be  retained,  and  let  vk  ^  ^  be  corresponding  eigenvectors. 


Define 


vi  *  <vi  vI>Ires<n'r,x<1'1) 

V,  -  (v,  v  _  )I  SR<-r)x(«-r-l+l) 

4  f  •  I!1  IT  IT 


(3.1.23) 


To  render  unobservable  select  g^  such  that  g^CV^O.  The  observ¬ 


ability  of  the  remaining  eigenvalues  may  be  checked  by  computing  the  vector 
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gjCVa* 

Let  be  the  i  desired  eigenvalues  to  be  assigned  to  At+1 

A 

by  f,.  By  the  dual  of  equations  (1.5),  (1.6)  f^  is  given  by  (SRQB^f^P 
where : 

n-r 


Po(X)*lu-AiUkS0  akX  ,  an_r  “  1 


S  “  T 


n-r-1  , 

ln-r  _  , 

°i  . “i  ) 

- ixn-r  *  _T 

•  «  t  «  1 

• 

„  »  • 

€  R  ,  S  =  T 

•  a 

n-r-1  * 

r 

•  m 

n-r  _  , 

ifi  . °i  i] 

ai  ’  *  *  *  CTi  1 

6Rixn-r+l 

(3.1.24) 


«  T 

P0(<Jl) 

• 

A 

1 

ixl  _ 

8ic 

P  = 

• 

• 

• 

P0(Ci)_ 

=  S 

a  , 

n-r-1 

• 

• 

al 

ao 

€R  ,  Q  - 

fi^i _ 

6R 


n-rxn-r 


R 


a  .  1 

n-r-1. 


*1 


€  R 


n-rxn-r 


If  g  eigenvalues  of  are  observable  with  §<i  then  g  should  be  used  in 
place  of  i  in  (3.1.24). 

Frequency  Domain  Solution 

Let  the  transfer  function  of  (A^B^C)  be  given  by  (3.1.9), 
(3.1.10),  and  let  the  eigenvalues  of  A^  be  the  first  SL  -  1  of  which 
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If 

f. » 


are  to  be  rendered  unobservable  by  the  selection  of  g^.  Then  must 
satisfy 


giN(ek)=0,  k=l,...,X-l 


(3.1.25) 


and  may  be  obtained  as  the  solution  of  g  JIT^  *  0  where  the  kth  column  of  M 


r  1  £-1 

is  a  linearly  independent  column  of  N(ek),  k  =  l,...,A-l.  If  t  j^J  is  a 


.  th 


sequence  such  that  the  jk  column  of  N(e^)  is  nonzero  then  the  columns  of 


MTr  may  be  computed  as 


N. 


1 

e. 


.  n-r-1 


if  e,  is  real 
k 


U. 


R?<V 

Re(ek’rl) 


i»(  ek) 

T "  /.  n-r-1 
lm(e 


) 


(3.1.26) 


if  (ek>  ek+i)  are  a  complex  pair 


where 


N. 

J 


lljO 


nijO 


nlj  (n-r-1) 


lXj (n-r-1) 


(3.1.27) 


To  assign  the  eigenvalues  gain  f^  is  computed  from 

the  dual  of  equations  (1.31)— (1.35) .  The  transfer  function  of  the  single¬ 
output  system  (A^.B^.g^C)  is 


n-r-1 


gtC(Xl-Ai)Bi  -  klQ  8iNkX  /  d(X) 


(3.1.28) 


A 
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Define  k  -  0, . . . ,n-r-l  and  let  p£  (X)  *  P1(X)P2  (A.)  where 


Pltt)-t2ia-ak)-k50dkXk,  <1.-1 

p2  ^  *  k*0  ^k  ^  '  'V-r-i  * 


(3.1.29) 


Then  f^  is  the  solution  of 


-  Y  ,  h  =  (hQ  h1  -  hn-r-i-l)J 


xo 

d  0  ....  0 

u  •  • 

•  • 

1 

d.,  d„  ' .  ‘ 

1 

\  °.  *0 

1 

1 

1 

d1  dl  '• 

4“l  .  •.*  do 

1 

1 

1  :  *’di 

1 

•  d.  : 

:  i-i  : 

•  e  • 

1 

1 

:  1  ''A 

:  :  *.  di-i 

•  •  • 

1 

•  • 

•  • 

x  1 

0  0  ....  1 

n-r-1 

,  Y 


o 

“0  ~ 
• 

c 

• 

• 

0 

do 

1 

dl 

• 

• 

n-r-l_ 

• 

• 

_di-l_ 

(3.1.30) 


If  ^  eigenvalues  of  A^  are  observable  with  5  <i  then  §  should  be  used  in¬ 
stead  of  i  in  equation  (3.1.30). 

In  the  discussion  of  Procedures  1  and  2  it  has  been  assumed  that 
X<n-r.  In  the  case  that  X^n-r  the  spectrum  of  A^  may  be  assigned  by 
state  feedback  and  the  algorithm  will  end  at  the  first  stage.  Let  p'g  R 
be  a  vector  assigning  a  desired  spectrum  to  A^  ■  A^ +b^P' .  Then  P  may  be 
taken  to  be  any  solution  of 


lxn-r 


PTA12  =  P' 


(3.1.31) 
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and  Che  degree  of  Che  required  compensaCor  will  be  one.  EquaCion  (3.1.31) 
always  has  a  soluCion,  being  a  sec  of  n-r  equacions  in  &  unknowns. 

If  (A^,b^)  is  noC  conCrollable  Che  algoriChm  may  scill  be  applied 
only  wich  Che  disCincCion  ChaC  ac  Che  iCh  sCage  any  unconCro liable  eigen¬ 
values  of  (A^.B^)  will  be  conCained  in  Che  specCrum  of  A^+^  in  addicion  Co 
Chose  placed  by  Che  designer.  This  may  be  used  Co  advancage  if  Che  unconCrolla- 
ble  eigenvalues  are  '’accepCable".  If  Chey  are  noC  Chen  ic  may  be 
advanCageous  Co  choose  a  differenc  ordering  of  Che  eigenvecCors  iv^i.^  of 
Che  orginal  sysCem.  Ic  should  be  noCed  ChaC  Che  sysCem  C^n_r»®n_r>^)  i-s 
conCrollable  since  Che  columns  of  Bn_r  span  Rn  r  [1].  The  algoriChm  will 
cherefore  evenCually  reach  a  scage  where  (A^B^C)  is  conCrollable  regard¬ 
less  of  Che  ordering  of  Che  eigenvecCors. 

In  summary,  Chis  algoriChm,  when  incorporaced  inCo  Che  meChodology 
for  designing  low  order  dynamic  regulacors,  decermines  Che  degree  and 
impliciCly  Che  parameCers  of  Che  required  compensaCor  ChaC  shapes  Che  encire 
specCrum  of  Che  resulcing  closed  loop  sysCem.  This  is  achieved  by  compuCing 
a  feedback  maCrix  P  as  a  sum  of  dyadic  produces  such  ChaC  Che  specCrum  of 
Ar*A^  +  BgPAj2  i-s  saeisfacCory.  The  specCrum  of  Che  CoCal  closed  loop  sys¬ 
Cem  is  Chen  deCermined  as  Che  specCrum  of  Ar  CogeCher  wich  chose  eigenvalues 
corresponding  Co  Che  selecCed  Nq  and  and  Che  parameCers  of  Che 

compensaCor  are  fixed  as  funcCions  of  P. 

3.2.  Review  of  Che  Design  MeChodology 

Before  considering  several  numerical  examples  in  Che  nexC  chapCer, 
iC  may  be  useful  Co  sunmarize  Che  design  meChodology  ChaC  has  been  presenCed. 


Because  the  information  structure  of  most  linear  systems  prohibits 
the  implementation  of  optimal  control  laws  based  on  state  feedback,  a  design 
criterion  has  been  defined  to  be  the  construction  of  suboptimal  control  laws 
that  retain  as  large  an  invariant  subspace  of  an  optimal  state  regulator 
as  possible  in  the  resulting  closed  loop  system. 

A  preliminary  study  of  the  system  should,  in  addition  to  an 
identification  of  the  controllability  and  observability  structure,  include 
an  analysis  of  the  possibility  of  satisfactorily  shaping  the  closed  loop 
spectrum  by  static  output  feedback.  As  will  be  seen  in  the  examples,  in¬ 
sight  helpful  in  the  selection  of  A^  and  A^  may  be  gleamed  from  such  an 
analysis. 

Having  defined  a  state  feedback  regulator  problem  through  the 
selection  of  weighting  matrices  Q,R,  the  Ricatti  equation  must  be  solved 
and  the  resultant  closed  loop  eigenvectors  and  eigenvalues  computed.  The 
possible  choices  for  will  be  no  more  than  C(n,r),  and  by  inspection  of 
the  optimal  eigenvectors  may  be  easily  identified  in  accordance  with  the 
requirement  that  complex  pairs  not  be  split  and  that  the  matrix  Y  be  in¬ 
vertible.  On  the  basis  of  the  spectra  of  the  resultant  matrices  A^,  r 
eigenvalues  must  be  selected  for  retention.  If  none  of  the  spectra  are 
acceptable,  and  a  compensator  is  to  be  designed,  then  this  choice  may  be 
guided  by  an  "identification"  of  those  eigenvalues  which  have  contributed 
most  to  the  unacceptability  of  the  spectrum  of  A^,  in  their  departure  under 
output  feedback  from  their  optimal  locations.  The  selection  of  A^  may  also 
be  based  on  information  obtained  from  the  preliminary  pole-placement 
analysis,  or  on  the  retention  of  dominant  eigenvalues. 
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When  this  first  part  of  the  design  is  completed,  a  decision  must 
be  made  whether  to  improve  the  dynamics  of  the  system  by  retaining  p 
additional  optimal  eigenvectors  through  the  introduction  of  a  dynamic  com¬ 
pensator  or  by  reducing  the  number  of  retained  eigenvectors  from  r  to  £<r. 
Because  the  preliminary  static  output  feedback  pole-placement  problem 
corresponds  to  the  case  .4  =  0,  unless  this  problem  had  a  satisfactory  solution 
the  design  of  a  dynamic  compensator  should  be  undertaken. 

If  a  compensator  is  to  be  designed,  the  remaining  n-r  eigenvalues 
must  be  ordered,  the  vectors  b^  computed,  and  the  pole-placement  problem  for 
the  triple  (Aj^BqjA^)  solved.  This  ordering  may  again  be  based  on  the 
desire  to  retain  at  their  optimal  locations  those  eigenvalues  most  contribut¬ 
ing  to  the  unacceptability  of  the  spectrum  of  the  matrix  corresponding  to 

the  selected  A  .  It  would  be  desirable  to  also  take  into  consideration  the 
r 

controllability  properties  of  the  pairs  (A^jB^),  but  a  convenient  criterion 
for  ordering  the  vectors  b^  to  enhance  the  solvability  of  the  pole-placement 
problem  is  unfortunately  not  available.  The  pole-placement  problem  is 
solved  by  computing  a  sequence  of  dyadic  feedbacks,  at  each  stage  increasing 
by  one  the  number  of  optimal  eigenvectors  retained  (as  well  as  the  number  of 
assignable  eigenvalues  of  Af  and  the  number  of  columns  of  Bq),  until  a 
satisfactory  tradeoff  is  achieved  between  the  spectrum  of  A^  and  the  dimension 
p  of  the  compensator.  At  each  stage  of  the  algorithm  certain  previously 
assigned  eigenvalues  are  chosen  to  be  retained  and  a  number  of  additional 
eigenvalues  of  Af  are  specified.  This  provides  the  designer  with  consider¬ 
able  freedom  to  meet  design  specifications  for  the  n-(r+p)  remaining  eigen¬ 
values  of  the  eventual  closed  loop  system.  In  particular  though  arbitrary 
pole-placement  for  the  matrix  will  not  be  possible  if  p<n-2r+l,  this 


freedom  may  be  used  to  place  the  eigenvalues  of  Ar  in  desired  regions  of  the 
complex  plane,  as  may  be  determined  for  example  by  minimum  damping  ratio 
requirements,  or  other  considerations. 

Once  the  pole-placement  problem  has  been  solved  satisfactorily, 
the  dimension  of  the  desired  compensator  has  been  determined  and  the  param¬ 
eters  H,D,K  ,K  may  be  computed  in  any  convenient  basis,  completing  the 
z  y 

design  procedure. 

The  computational  aspects  of  the  design  procedure  are  straight 
forward,  involving  only  the  solution  of  eigenvector  equations,  the  solution 
of  linear  systems  of  equations,  and  associated  algebraic  manipulations. 
Following  a  transformation  of  basis  to  bring  the  output  matrix  to  the  form 
[Ir  0] ,  the  solution  of  a  Ricatti  equation,  and  the  determination  of  the 
optimal  eigenvectors,  at  most  C(n,r)  eigenvalue  calculations  are  required 
to  compute  the  spectra  of  the  matrices  A^.  The  solution  of  the  pole- 
placement  problem  by  state  space  procedures  requires  at  each  stage  the 
computation  of  the  left  (or  right)  eigenvectors  and  the  characteristic  equa¬ 
tion  of  A^,  the  solution  of  a  homogeneous  system  of  equations  of  order  r 
(or  i)  to  determine  f^  (or  g^)  and  finally  the  solution  of  an  inhomogenous 
system  of  equations  of  order  i  (or  r)  to  find  g^  (or  f^).  As  the  pole- 
placement  problem  is  solved  interactively,  allowing  for  the  repeated  execution 
of  each  stage  of  the  algorithm  until  the  designer  is  satisfied  with  the 
spectrum  of  A^_^,  the  solution  may  be  costly  if  the  dimension  n-r  of  A^  is 
large  and  many  repetitions  are  employed. 

The  final  computation  of  the  parameters  of  the  compensator  in¬ 
volves  only  the  algebraic  manipulation  of  matrices. 


The  design  methodology  will  be  illustrated  by  several  nontrivial 


examples  in  the  next  chapter. 


3.3.  Extension  to  Stabilizable  Systems 

Although  it  has  been  assumed  that  the  triple  (A,B,C)  is  controlla¬ 
ble  and  observable,  this  restriction  may  be  relaxed.  In  this  section  the 
application  of  the  design  methodology  to  the  class  of  observable,  stabiliz¬ 
able  systems  (A,B,C)  is  considered  and  it  is  shown  that  for  such  systems  the 
dynamic  compensator  of  section  2.3  possesses  a  separation  property. 

Let  the  system  be  represented  in  the  canonic  form 


A  n  x  1  A  ilxI 
xL  €R  ,  .€  R 

(3.3.1) 


A 

y  = 


where  A^  is  a  stable  matrix,  the  triple  (A22,B2,C2^  is  controllabl-e  and 
observable,  and  the  pair  *-s  observable. 

Consider  first  the  solution  of  the  state  feedback  regulator  problem 
for  this  system.  It  may  be  assumed  without  loss  of  generality  that  R  =  I 
(under  a  transformation  u.=Vru).  Compatibly  partitioning  the  solution  M 
of  the  Ricatti  equation  (2.1.2)  gives  the  three  equations: 

A22M22  +  **22^22  _M22^2^2  M22  +Q22  =  0 
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^11MU  +  M11^11  +  ^2lM2l+M2l^2l  _M2l®2®2M2l+^lP  *  0 


the  first  of  which  is  the  Ricatti  equation  for  the  subsystem  (A^.B^)  wit^ 
penalty  Q22-  The  optimal  control  is  then 

,  («j_  K2)  =  (B^M^  j  B2M22)  (3.3.3) 

and  the  closed  loop  matrix  is 


B2K1 


(3.3.4) 


Thus  the  optimal  control  has  the  property  that  the  controllable  subsystem 

^  A 

F22  iS  t*le  °Ptimal  closed  loop  system  matrix  for  the  pair  (A22»®2^  with 
penalty  Q22.  T^e  °Ptimal  control  shapes  the  dynamics  of  A22  exactly  the 
same  as  if  there  were  no  driving  uncontrollable  subsystem  (^^  =  0),  but  also 
expends  energy  in  shaping  the  eigenvectors  (and  therefore  the  response)  of 
the  uncontrollable  state  variables.  This  is  true  even  if  the  uncontrollable 
states  are  not  penalized  (Q  =  dg(0  Q22>>* 

Let  the  controllable  and  uncontrollable  eigenvalues  of  the  optimal 
closed  loop  system  be 


(3.3.5) 


The  optimal  controllable  eigenvectors  may  be  written  explicitly  as 
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(3.3.6) 


and  the  uncontrollable  eigenvectors  have  the  form 


c  ?  T  c 

ui  Anui  '  ui 


x  lj . . . ) n. 


(3.3.7) 


c  _  c  c  c 

vi  F12ui  +  F22vi  =  vi  ai 


In  order  to  apply  the  methodology  for  the  design  of  dynamic  compensators  it 
is  necessary  that  the  matrix  Y  in  equation  (2.3.3)  be  invertible.  In  an 
arbitrary  basis  this  requires  the  selection  of  r  eigenvectors  such  that  the 
matrix  CU  appearing  in  equation  (2.2.11)  is  invertible.  Assuming  that  a 
design  criterion  is  the  retention  of  as  many  optimal  controllable  eigen¬ 
vectors  as  possible,  the  maximum  number  of  controllable  eigenvectors  which 
may  be  retained  in  Ar  is  therefore  where 

r2  *  rankO^tVj0  vnC])  (3.3.8) 

A  A  C  n2 

Since  the  pair  (A^jC^)  is  observable  and  the  vectors  v^  span  R  ,(3.3.8) 
simplyfies  to  r^rank^).  Thus  r2  controllable  eigenvectors  may  be 
retained  in  and  the  remaining  r^ = r-r2  eigenvectors  must  be  selected 
from  the  uncontrollable  subsystem. 

To  obtain  the  separation  property  for  the  compensator,  the  output 
matrix  is  first  transformed  to  {I  0],  Let  rank  ((?2)=r2  and  let  S  be  any 
full  rank  output-space  transformation  such  that 


6] 


,  rank  (C22)  * (3.3.9) 

Then  under  the  transformation  y = Sy  the  output  matrix  takes  the  form 

r  xr  r  xr 

'  tfll€R  »  C22  €r  »  rl  +  r2  =  r 

(3.3.10) 

Now  let  and  T2  be  any  matrices  such  that  the  transformation  x’  =T£  is 
invertible  where 


T 


Under  this  transformation  the  system  (3.3.1)  becomes 


(3.3.11) 


(3.3.12) 


Introduce  the  partitions: 
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A1 

A11 


A11  A12 

A31  A32 

A33  A34 

A  '  = 

A2l 

A22  * 

_^2l  A22__ 

_\l  A42_ 

_^43  A44_ 

B2 


c 

B3 

x'  - 

X1 

B4_ 

1 

o’ 

X2_ 

-I  2 


x;  = 


(3.3.13) 


-  r.xl  -  n.-r,xl  „  r,xl  „  n,-r,xl 

1  ,  1  1  ,  Xic6r  2  ,  x/cr2  2 


Then  under  the  permutation  of  states: 


r 

* 


c 

X1 

|— 

—m 

I 

rl 

0 

0 

0 

c 

X1 

c 

= 

0 

0 

S 

0 

x2 

0 

Inrri 

0 

0 

X2C 

0 

0 

0 

ru  “*  rn 

2  2 

(3.3.14) 


the  system  is  represented  as: 
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a  c 
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a  C 

vi  »  i  =  1»  •  •  •  >n-i 


(3.3.16) 


At  the  first  stage  of  the  design  it  is  necessary  to  select  Y,Z  and 
to  compute  ct(A^).  From  (3.3.8)  and  the  structure  of  the  output  matrix  in 
(3.3.15)  it  follows  that  at  most  r^  controllable  eigenvectors  may  be  retained 
in  Y.  In  accordance  with  the  design  criterion  of  retaining  as  many  optimal 
controllable  eigenvectors  as  possible,  let: 


Z  = 


Y  C 

1 

O 

v- 

r  1  rl1 

•••  \  1  » 

,,  c  r  1 

Y2  -  [03  ... 

V 

r  1  'l, 

[«3  ^3  ] 

0 

*i5- 

1  rl 

[^2  • • •  ^  1  > 

_  c  r  1 

Z2  =  [04  ... 

h' ° 

Z2l 

’  hz- 

1  rl 

K  •••  i 

(3.3.17) 


u .  = 
X 


i+r. 


vi  = 


i+r„ 


,  r  —  1 ,  •  • • , n2  ~  r 2 


i+(n2-r2) 


i+r. 


i+r. 


i+(n2*r2^ 


i+r, 


i+rl 

|^_4  __! 


,  i-1... 


•Vri 


Then  from  equation  (3.1.1): 
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N  *  ZY 
o 


N 


H2l  No 


“o'  *  2lCfflC)'1 


H21  -  z2c(Ylc)'1-z2c(Y2c)‘1Y2c(Y1i:)"1 


80C  -  Z2CW2C):1 


/ 

0 

\ 

.  c 

i * l,...,n2-r2 

LO 

\ 

i 

b! 

i 

i  =n2~r2+l, . . .  ,n-r 

J 

.  c  „  i+t2  „  C  .  1+I2 

b .  ■  0,  - N  0. 

13  04  (3.3.18) 


y  hc- 


i+r-n»  —  i+r-n„ 
“*2  -8o  ”l  2 


b?  =  ID 
i  4 


i+r-n  i+r-n 

"N21  “l 

,  c  i+r-n2 
+  Nq  ^3 


A22  "No  A12 


A42  ‘  (N21A12  +No  A32)  A44  “  No  A34 


—  nl”rl  n2-r2 

Note  that  the  vectors  and  span  R  and  R  respectively  (see 

remarks  following  (2.3.6)). 

Assuming  that  only  controllable  eigenvalues  are  to  be  retained  in 
c  c 

A  ,  define  B.  = [b.  . ....  b.  ].  Then  the  pole-placement  problem  is  that  of 

P  *  I  •  »  1 

finding  a  feedback  P=[P^  P2J  to  satisfactorily  assign  the  spectrum  of 


0 

^12 

0 

A  «  A.  + 
r  1 

<P1  P2> 

_Bi_ 

A32 

A34_ 

(3.3.19) 


£ 

It  should  be  noted  that  Af  contains  an  uncontrollable  subsystem  A22  “N0  A^2 
whose  spectrum  is  the  n^  -  r^  uncontrollable  eigenvalues  not  retained  in  A^. 
The  remaining  spectrum  of  Af  may  be  shaped  by  solving  a  pole -placement 
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c 

problem  for  the  triple  (A^“Nq  B^,  A^4) . 

Assuming  the  solution  of  this  reduced  pole-placement  problem  is 
P2,  the  required  compensator  parameters  may  be  computed  from  equation  (2.3.7) 
with  P*[0  P2].  Let: 

1 +r2  r2+P  1+r  P+r 

U2  =  [03  ...  03  ],  v2“[04  •••  04  ] 

(3.3.20) 


where  A  ,  A  correspond  to  the  r.  uncontrollable  and  r.  controllable 
rl  r2 

eigenvalues  retained  in  the  first  stage.  Then  it  may  be  shown  that: 


H  =  H 
o  o 


HqC  =  IAp-LCAr2(Y2C)'1U2][I+P2U2] 


Do  =  [DoC  “o*1  *DoC  "  (L\  "  VC)(Y2C)“l(I+U2P2) 


»*  *  [LCA  -  LCA  (Y2C)_1Y2C]  (YlC)’L 


LC  *  [I  +  P2U2]_1P2Y2C  ,  LC  =  [I  +  P2U2]_1P2Y2C 


(3.3.21) 


N 


NpC  *  IV2  -N0CU2HI+P2U2]  -b;[I+P2U2] 


N 


N21  "r 


-8oC-<V2-“oCVP2-,,oC-i,;P2 


Partitioning  the  optimal  state  feedback  gains  in  this  basis  as  K  =  [K,  Kj,], 

* [K^C  K^C] ,  K2  “  K^] ,  the  suboptimal  feedback  gains  are  given  by: 
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"  h  +  V*  '  “=/  Ky°l 


Kz  -  V 


K 

P  z 


c  c  "S  c  c 
K  =  Kx  +K2  Nq  +  K,  Nn 


(3.3.22) 


c  c  c  c 
Ky  =  K,  +K,  Nr 


c  c  c 
K  -  K_  N 
z  z  p 


It  is  noted  that  Ky  is  independent  of  the  solution  P  of  the  pole-placement 
problem. 

By  the  superposition  principle  for  linear  systems,  the  controller 
may  be  represented  as  two  compensators  in  parallel: 


e  ,  c 

u  =  u  +  u 


(3.3.21a) 


uC*K  CvC  +  K  CyC,  vc  =  H  V+D  V  (3.3.21b) 
z  y  *  o  o 

7T  c  c  c  c  .c  c  3  c  c 
u  -K  v  +K  y  .  V  =H  v  +D  y  (3.3.21c) 

z  y  *  o  o 


The  compensator  given  in  (3.3.21b)  is  precisely  that  compensator  which 
would  have  been  designed  if  the  controllable  and  observable  subsystem 
^A22,B2’^Ir  ’  were  operating  in  isolation.  The  second  compensator 

(3.3.21c)  represents  a  modification  in  the  control  scheme  due  to  the  presence 
of  the  driving  uncontrollable  subsystem  Ajj.  A  block  diagram  of  this  con¬ 
troller  is  given  in  Figure  3.3.1. 

In  summary  it  has  been  shown  that  the  design  methodology  may  be 
applied  to  stabilizable  systems,  and  that  the  resultant  compensator  satisfies 
a  separation  property.  It  is  noted  that  the  restriction  that  the  system  be 
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observable  may  be  relaxed  provided  Che  designer  ensures  that  Che  macrix  Y 
in  (2.3.3)  is  invercible. 

3.4.  SofCware  Supporc 

The  examples  in  Che  nexC  chapcer  were  solved  using  an  inCeracCive 

Foreran  program  DSGN'FOR  wriCCen  Co  implemenC  Che  meChodology  described  in 

Chis  chesis.  The  algorichm  solving  Che  pole-placemenc  problem  was  implemenCed 

using  che  state  space  formulacions  of  Procedures  land  2.  A  flow  chare  for 

Che  program  is  given  in  Figure  3.4.1. The  user  is  assumed  Co  have  previously 

obCained  Che  solution  M£  to  the  Ricatti  equation  and  to  have  stored  on  disc 

in  column  major  order  the  matrices  A,  B,  R,  Mc  in  a  basis  representation  in 

which  C*=  [I  0] .  The  matrix  Q  is  not  needed, 
r 

All  eigenvalue  and  eigenvector  calculations  are  performed  using 
the  IMSL  subroutine  EIGRF,  which  returns  an  estimate  (on  the  user's  console) 
of  the  accuracy  of  the  computed  eigenvectors.  The  eigenvector  computation 
is  satisfactory  if  this  estimate  is  less  than  one,  fair  if  between  one  and 
100,  and  poor  if  over  100.  The  eigenvalues  are  then  ordered  by  increasing 
real  part,  and  their  corresponding  eigenvectors  are  normalized  to  unit 
length.  In  the  case  of  a  complex  eigenvalue,  the  corresponding  eigenvector 
is  represented  in  its  real  and  imaginary  parts  and  is  only  determined  to 
within  a  complex  multiplicative  constant  cos  9+j  sin  9.  All  matrix 
inversions  are  performed  using  the  IMSL  subroutine  LINV2F  which  also  returns 
an  error  code  which  should  be  zero. 

In  the  solution  of  the  pole-placement  problem  for  Ar=A^  +  BoPA^2 
the  vectors  f^  are  computed  as  zero  eigenvectors  of  the  matrices  obtained  by 
augmenting  the  homogeneous  systems  V^B^f^  =  0  with  a  zero  row.  The  vectors 
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are  computed  as  the  solutions  of  the  square  inhomogeneous  systems 
(CQRS)T  gT“PT  using  the  subroutine  LINEQ  available  from  [34].  The  routine 
returns  a  condition  number  for  the  coefficient  matrix  which  is  typed  on  the 
user's  console. 

All  output  is  to  the  lineprinter  only,  though  the  user  has  the 
possibility  of  suppressing  the  listing  of  any  undesired  data. 
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CHAPTER  4 
EXAMPLES 


4.1.  Saturn  V  Booster  Mudel 

In  this  example  a  second  order  compensator  is  designed  for  a 
seventh  order,  single-input,  two-output  model  of  a  Saturn  V  booster.  The 
model  has  appeared  in  three  articles  [23] , [24] , [35]  on  output  feedback  pole- 
placement  in  which  numerical  algorithms  were  employed  to  stabilize  the 
system.  In  [24]  and  [35]  the  real  part  of  the  least  stable  eigenvalue 
of  the  closed  loop  system  was  minimized  by  two  different  methods  and  in 
[23]  the  eigenvalues  of  the  closed  loop  system  were  constrained  to  lie  in 
a  prescribed  region  of  the  complex  plane. 

The  model  is  given  by  x^Ax+Bu,  yBCx  where  A,B,C  are  given  in 
Table  4.1.1. 

As  a  preliminary  analysis  the  possibility  of  stabilizing  the 
system  by  output  feedback  was  considered.  In  both  [23]  and  [24]  the  real 
part  of  the  least  stable  eigenvalue  was  required  to  be  less  than  -0.07,  and 
this  resulted  in  a  damping  ratio*  of  less  than  0.02  (p =  89°) .  Relaxing 
the  requirement  that  all  the  eigenvalues  of  the  closed  loop  system  lie  to 
the  left  of  a =-0.07,  and  attempting  instead  to  increase  the  damping  ratio, 
the  results  given  in  Table  4.1.2  were  obtained  using  the  pole-placement 
subroutine  of  the  compensator  design  software.  The  damping  ratio  has  been 
increased  to  0.1  (p=84°)  while  the  least  stable  eigenvalue  has  been  shifted 


The  damping  ratio  of  a  stable  matrix  A  is  here  taken  to  mean  the 
smallest  damping  ratio  of  all  the  complex  eigenvalues  of 
A:  5=  min  {-o//a^+ cu^}.  The  associated  angle  n  *  cos“^£  is  measured 

o+ju£a(A) 

from  the  negative  real  axis. 


J 


Table  4.1.1.  System  matrices  for  the  Saturn  V  booster  model 
The  model  is  x  =  Ax+Bu,  y  =  Cx 

0  1  0  0  0  0  0  1  'c 

0  0  0.2  -0.65  -0.002  2.6  0  C 

-0.014  1  -0.041  0.0002  -0.015  -0.033  0  C 

A  =  0  0  0  0  1  0  0B-C 

000  -45.0  -0.13  255.0  0  C 

0  0  0  0  0  0  1  C 

0000  0  -50.0  -10.0  ] 


0 

0.2 

0 

-0.65 

0 

-0.002 

0 

2.6 

0 

0 

-0.041 

0.0002 

-0.015 

-0.033 

0 

0 

0 

1 

0 

0 

0 

-45.0 

-0.13 

255.0 

0 

0 

0 

0 

0 

1 

0 

0 

0 

-50.0 

-10. 

)  0  0 

0  0 

0  0 

0  0 

Table  4.1.2.  Comparison  of  static  output  feedback  compensator 
designs  for  the  Saturn  V  boosters  Solution  using 


Here  F  »  A-BKC 


1. 

2. 

3. 

o(F)  4. 


Sirensa  and  Choi 

Miller,  et.  al. 

Solution  using 

PPL  subroutine 

(-20.31,  -16.56) 

(-26.68,-16.27) 

(-152.541,-42.623) 

(-4.841,  5.433) 

(-4.823,  5.401) 

(-4.340,  6.018) 

(-4.841,-5.433) 

(-4.823,-5.401) 

(-4.340,-6.018) 

(-0.125,  0.497) 

(-0.118,  0.642) 

(-0.471,  4.683) 

(-0.125,-0.497) 

(-0.118,-0.642) 

(-0.471,-4.683) 

(-0.098,  0  ) 

(-0.105,  6.204) 

(-0.250,  2.400) 

(-0.070,  6.204) 

(-0.105,-6.204) 

(-0.250,-2.400) 

(-0.070,-6.204) 

(-0.078,  0  ) 

(-0.050,  0  ) 

0.0113 

0.0169 

0.1001 

89° 

89° 

84° 

to  -0.05.  The  gains  required  to  achieve  this  solution  are  large  given 

the  slight  improvement  in  the  pattern  of  the  closed  loop  spectrum  over 

the  solutions  of  [23]  and  [24]. 

In  view  of  these  results  a  linear  quadratic  regulator  problem 

was  formulated  with  the  expectation  that  a  low  order  compensator  would  be 

required  to  satisfactorily  shape  the  dynamics  of  the  final  closed  loop 

T 

system.  Solving  the  state  feedback  regulator  problem  with  Q=aC  C  and  R=1 

3 

for  a  few  values  of  a  indicated  that  unless  a  were  of  the  order  of  10  or 
4 

10  ,  the  optimal  solution  would  possess  a  complex  pair  with  a  small  damping 
ratio,  and  this  led  to  the  selection  of  '  / 

J 500  0  * 

R  -  0.01  Q  =  C  C  .  (4.1.1) 

_  0  100. 

The  solution  of  the  state  feedback  regulator  problem  for  this  choice  of 
Q  and  R  is  given  in  Tables  4.1.3  and  4.1.4.  In  comparison  to  the  minimum 
energy  solution  (Q=0,  R=I)  which  simply  reflects  the  unstable  eigenvalues 
about  the  imaginary  axis,  this  optimal  regulator  has  perturbed  the  eigen¬ 
values  at  -5+j5  and  -0.014  only  slightly  to  X^,  X^,  and  X^,  whereas  the 
pair  -0.065 + j 6. 708  has  been  moved  to  the  pair  X  , X  and  the  two  real 

->  D 

eigenvalues  -0.475  and  -0.420  have  formed  a  complex  pair  and  moved  to 
X^,X^  (see  Figure  4.1.1).  It  is  noted  that  the  least  stable  eigenvalue  is 
very  near  an  invariant  zero  of  the  system.* 

Since  r=2  and  the  optimal  spectrum  of  F  contains  three  complex 
pairs  of  eigenvalues  there  are  only  three  choices  for  A  .  Computing  the 

* 

A  survey  of  the  literature  on  zeros  of  linear  time- invariant 
multivariable  systems  is  available  in  [36]. 


i J 


Table  4.1.3.  Open  loop  eigenvalues,  optimal  closed  loop  eigenvalues, 
and  invariant  zeros  of  the  Saturn  V  booster 


Closed  Loop  Spectrum  Invariant  Zeros 
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Table  4.1.4.  Optimal  state  feedback  regulator  solution 


The  Riccati  solution  is  M  = 


620.881 

326.193 

11.062 

-11.812 

-4.778 

-40.707 

-2.235 

326.188 

305.360 

23.048 

-9.677 

-4.775 

-48.934 

-2.826 

11.061 

23.047 

10.431 

-0.461 

-0.384 

-4.713 

-0.289 

-11.812 

-9.677 

-0.461 

0.510 

0.140 

0.336 

-0.013 

-4.778 

-4.775 

-0.384 

0.140 

0.082 

0.988 

0.054 

-40.705 

-48.933 

-4.713 

0.335 

0.988 

18.839 

1.158 

-2.235 

-2.826 

-0.289 

-0.013 

0.054 

1.158 

0.082_ 

K  =  R_1BTM 

c 

= 

223.486  - 

282.557 

-28.919 

-1.343 

5.370 

115.817 

8.211] 

closed  loop  matrix 

is  F 

=  A-BK  = 

0.000 

1.000 

0.000 

0.000 

0.000 

0.000 

o.ooo' 

0.000 

0.000 

0.200 

-0.650 

-0.002 

2.600 

0.000 

-0.014 

1.000 

-0.041 

0.0002 

-0.015 

-0.033 

0.000 

0.000 

0.000 

0.000 

0.000 

1.000 

0.000 

0.000 

0.000 

0.000 

0.000 

-45.000 

-0.130 

255.000 

0.000 

0.000 

0.000 

0.000 

0.000 

0.000 

0.000 

1.000 

223.486 

282.557 

28.919 

1.343 

-5.370 

-165.817 

-18.211 

The  optimal  eigenvectors  are: 


VX2 


'  o.ooi 

"  0.001 

'-0.002' 

0.000 

-0.021 

0.000 

-0.010 

-0.003 

0.004 

-0.012 

0.037 

-0.018 

0.000 

-0.001 

0.000 

0.000 

-0.024 

0.005 

0.026 

±3 

0.138 

VA  A  “ 

-0.122 

±j 

-0.014 

Vl5'\f 

0.204 

±j 

-0.401 

-0.755 

-0.587 

J  » 

0.388 

-0.903 

-0.030 

0.872 

0.030 

0.023 

0.002 

0.017 

0.033 

-0.076 

.-0.253 

_  0.018 

-0.134 

-0.022 

.  0.004 

0.162 

0.000; 
0.000 
0.685 
0.717 
0.033 
0.126 
L-0.006 


matrices  and  their  spectra  shows  that  in  each  case  is  unstable  and 

so  as  expected  a  compensator  will  be  required  to  stabilize  the  system. 

The  data  for  this  first  stage  of  the  design  are  given  in  Table  4.1.5,  and 

Figure  4.1.2.  In  all  three  cases  it  should  be  noted  that  there  is  an 

eigenvalue  very  near  the  invariant  zero  at  -0.046.  It  is  therefore  expected 

that  in  designing  a  dynamic  compensator  the  spectrum  of  A^  will  "likely" 

contain  an  eigenvalue  near  this  location.  Hence  A  should  not  contain 

P 

A^  = -0.046  but  rather  eigenvalues  which  depart  more  from  their  optimal 
locations  under  static  output  feedback  (see  Figure  4.1.2). 

This  was  confirmed  as  it  was  not  possible  to  design  a  first  order 
compensator  by  stabilizing  A  f°r  any  of  the  orderings  ^jc**k+l,'S^  ’ 
k*  1,3,5.  Using  two  columns  of  Bq  however,  A^_  was  easily  stabilized  for 
several  choices  of  and  A^.  Two  compensators  based  on  the  orderings 
^1^2^5^6  an<*  ^3^4^5^6  discussed  here.  In  both  designs  the  spectrum 

of  A^  was  shaped  in  accordance  with  two  criteria:  first  that  as  many  eigen¬ 
values  as  possible  be  placed  at  the  locations  of  the  n-r-p  unretained 
optimal  eigenvalues,  and  second  that  the  damping  ratio  of  Aq  (as  defined 
above)  be  no  less  than  that  of  the  optimal  solution.  In  view  of  the  latter 
requirement  it  is  noted  that  the  damping  ratio  of  the  optimal  closed  loop 
system  is  5*  0.289  (T|—  73°). 

Retaining  A^^  ant*  ordering  the  remaining  eigenvalues 
A^AgA^A^A^,  with  the  real  eigenvalue  \y  placed  third  in  anticipation  of 
designing  either  a  second  or  a  third  order  compensator,  the  pole-placement 
subroutine  was  used  to  place  the  spectrum  of  Since  the  first  row  of 

is  all  zeros  it  was  necessary  to  introduce  P*PT,  T*  (0,1)  and  to  solve 
the  pole-placement  problem 


Table  4.1.5.  Data  for  various  choices  of  Ar 


Based  on  retention  of  A^.A^: 
The  spectrum  of  is: 


1. 

( 

-0.194, 

7.095) 

-0.771 

-0.113 

2. 

( 

-0.194, 

-7.095) 

128.864 

8.119 

3. 

( 

-0.065, 

0.000) 

N  = 

-374.840 

45.949 

4. 

( 

0.247, 

0.729) 

o 

14.231 

-1.854 

5. 

( 

0.247, 

-0.729) 

85.580 

33.161 

-0.018 

-0.073 

-0.015 

0.260 

0.000 

-1.624 

5.278 

1.016 

-21.111 

0.000 

A, 

s 

-9.190 

-15.133 

-0.038 

135.534 

0.000 

1 

0.371 

-1.205 

-0.004 

4.820 

1.000 

-6.632 

21.555 

0.066 

-136.220 

-10.000 

Based  on  retention  of  AyA^: 
The  spectrum  of  is: 


1. 

(  -4.767, 

3.087) 

"  -0.168 

-0.021 

2. 

(  -4.767, 

-3.087) 

80.281 

1.251 

3. 

(  -0.047, 

0.000) 

N  • 

-79.796 

74.516 

4. 

(  2.010, 

2.973) 

0 

-4.519 

-1.401 

5. 

(  2.010, 

-2.973) 

_  89.397 

1.939 

‘  -0.037 

-0.014 

-0.015 

0.022 

0.000 

-0.250 

0.813 

1.003 

-3.252 

0.000 

A,  * 

-14.903 

3.436 

0.019 

61.257 

0.000 

1 

0.280 

-0.911 

-0.003 

3.643 

1.000 

.  -0.388 

1.260 

0.004 

-55.041 

-10.000 

Based  on  retention  of  A^jA^: 
The  spectrum  of  A^  is: 

1.  (  -5.565,  8.109) 

2.  (  -5.565,  -8.109) 

3.  (  -0.050,  0.000) 

4.  (  2.261,  4.187) 

5.  (  2.261,  -4.187) 


0.016  -0.184  -0.016 

-4.616  15.002  1.046 

10.076  -77.746  -0.231 

-0.880  2.861  0.009 

1.869  -6.075  -0.019 


0.630 

30.743 

-86.803 

6.124 

.-16.553 

0.704 

-60.009 

385.985 

-11.443 

-25.702 


-0.284 

23.080 

-50.379 

4.401 

-9.346 

0.000 

0.000 

0.000 

1.000 

-10.000 


J 


Based  on  retention  of  X, A,  Based  on  retention  of  X  ,X  Based  on  retention  of  X 


Figure  4.L.2.  Spectra  of  suboptimal  closed  loop  systems  under  static  compensation 
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A  =  A,  +  B  f(TA.,)6r  .  (4.1.2) 

r  1  o  12 

Thus  at  the  ith  stage  of  the  algorithm  only  i  poles  could  be  arbitrarily 

assigned.  Since  at  least  two  columns  of  B  were  to  be  used,  at  the  first 

o 

stage  a  real  pole  was  arbitrarily  placed  at  -2  using  Procedure  1.  Then 
using  Procedure  2  a  complex  pair  was  placed  at  -1+  j3.5  (5  =  0.275, 
n  -  74°)  resulting  in  another  complex  pair  at  -1.625  +  j5.41  (5=  0.288, 
n=73°)  and  a  real  pole  very  near  A^  (and  the  invariant  zero  at  -0.0462). 
Attempts  at  placing  poles  near  A^,A^  resulted  in  an  unstable  A^ .  In  view 
of  the  criteria  above,  this  solution  of  the  pole-placement  problem  was 
satisfactory  and  the  degree  of  the  required  compensator  taken  to  be  p=2. 
Data  for  the  pole-placement  problem  is  given  in  Table  4.1.6. 

The  parameters  of  the  compensator  are: 


1.959] 


-2:2m]  kz-  ‘-1-221  i-959) 

9o  *  [-1212 : 909  Ky-  [-36.437  -30.255] 


(4.1.3) 


and  the  total  closed  loop  matrix  is 


-6.258 

-24.226 

0.315 

-2.602 

-267.925 

-1212.909 

-31.216 

-59.339 

0.000 

0.000 

0.000 

0.000 

0.000 

0.000 

0.000 

1.000 

0.000 

o'.OOo 

0.000 

0.000 

0.000 

0.000 

0.200 

-0.650 

0.000 

0.000 

-0.014 

1.000 

-0.041 

0.000 

0.000 

0.000 

0.000 

0.000 

0.000 

0.000 

0.000 

0.000 

0.000 

0.000 

0.000 

-45.000 

0.000 

0.000 

0.000 

0.000 

0.000 

0.000 

1.221 

-1.959 

36.437 

30.255 

0.000 

0.000 

0.000 

0.000 

o.ooo' 

0.000 

0.000 

0.000 

m 

0.000 

0.000 

0.000 

-0.002 

2.600 

0.000 

-0.015 

-0.033 

0.000 

1.000 

0.000 

0.000 

-0.130 

255.000 

0.000 

0.000 

0.000 

1.000 

0.000 

-50.000 

-10.000. 

81 


Table  4.1.6.  Data  for  pole-placement  problem 
based  on  ordering 


fx  =  [1.000] 

gx-  [1.621] 


g2=  [1.000] 


p 


0.000  4.651 

0.000  21.055 


The  spectrum  of  A 2  is: 


1. 

(-2.000, 

0.000) 

2. 

(-0.039, 

0.000) 

3. 

(-0.038, 

7.130) 

4. 

(-0.038, 

-7.130) 

5. 

(  1.107, 

0.000) 

The 

spectrum 

of  =  A 

1. 

(-1.625, 

5.407) 

2. 

(-1.625, 

5.407) 

3. 

(-1.000, 

3.500) 

4. 

(-1.000, 

-3.500) 

5. 

(-0.055, 

0.000) 

-0.036 

0.003 

2.624 

-0.279 

■ 

-9.639 

1.742 

0.402 

-0.111 

_  0.585 

0.731. 

Wi 

81(TA12) 

= 

-0.030 

-0.035 

-0.015 

0.108 

0.000 

-0.773 

2.512 

1.008 

-10.049 

0.000 

-12.316 

-4.974 

-0.007 

94.898 

0.000 

0.501 

-1.629 

-0.005 

6.516 

1.000 

- 

-6.443 

20.939 

0.064 

-133.755 

-10.000 

- 

a3*  a2+  tbib2  f 2 

82^A12^ 

3 

-0.038 

-0.009 

-0.015 

0.003 

0.000 

-0.360 

1.170 

1.004 

-4.679 

0.000 

-10.818 

-9.841 

-0.022 

114.363 

0.000 

0.276 

-0.896 

-0.003 

3.586 

1.000 

-3.009 

9.781 

0.030 

-89.123 

-10.000 

is 
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Referring  to  Table  4.1.7  and  Figure  4.1.3  it  is  verified  that  the  compen¬ 
sated  system  has  retained  from  the  optimal  state  regulator  a  four 
dimensional  subspace  spanned  by  the  eigenvectors  corresponding  to 
Ag.  Also  the  compensator  is  open  loop  stable. 

Recalling  that  H ,  D,  and  are  determined  only  up  to  a  similarity 

transformation  under  W  ,  the  pair  (D,K  )  may  be  balanced  by  introducing 

p  z 


H  -  W  H  W_1,  D  =  W  D  , 
z  pop  p  o 


K 


K  W-1. 
z  p 


(4.1.5) 


For  example,  if  W  I.  _  then  H  and  K  are  unchanged  and 

p  ZP  Z*Z  o  y 


D  = 


-10.717  - 

-48.516 


-1.249  ] 
-2.374 J 


Kz  =  [-30.525  48.975] 


(4.1.6) 


For  purposes  of  comparison  a  design  is  given  based  on  the  ordering 
A^,A^, A^,A^,A^,A^,A2.  The  pole-placement  problem  was  solved  by  first 
placing  a  pole  arbitrarily  at  -1.0,  and  then  in  the  second  stage  of  the 
algorithm,  placing  a  complex  pair  at  -1.0+ j 3.0.  This  resulted  in  a 
spectrum  for  which  met  the  criteria  given  before,  and  again  placed  a  real 
pole  near  A^.  The  data  for  the  pole-placement  problem  are  given  in  Table 
4.1.8.  The  parameters  of  the  resultant  compensator  are: 


f  -30.530  4.925] 

o  L  -68.884  8.858 J 

f  -970.743  224.189] 

0*1-7204.778  582.655 J 


K  =  [  -15.881  3.919] 

z 


K  -  [-544.318  76.387] . 

y 


(4.1.7) 


Again  the  matrices  D  and  K 

o  z 


6 

o 


-80.895 

-183.732 


may  be  scaled  to  give  (W 


18 

48 


.682] 
.555  J 


[-190.572 


“  12  X2x2): 
47.028], 


(4.1  .8) 


Comparing  the  two  compensator  designs  it  is  seen  that  that  based  on 
retention  of  ^*2*5^6  resu^te<^  feedback  gains  more  than  an  order  of 
magnitude  smaller  than  those  obtained  in  the  design  based  on  retention  of 


Table  4.1.7.  Spectra  of  closed  loop  compensated  Saturn  V 
booster  and  of  open  loop  compensator 


a(A  )  5. 
c  , 


o(v  i: 


Design  based  on  retention 

of  A- A0 A c A- 
12  5  6 

Design  based 
of  XX 

on  retention 
4X5X6 

(-5.106, 

4.483) 

:  X 

(-18.750, 

0.000) 

(-5.106, 

-4.483) 

:  \\ 

(  -2.908, 

0.000) 

(-1.757, 

0.820) 

:  X25 

(  -2.305, 

7.648) 

x3 

(-1.757, 

-0.820) 

■  4 

(  -2.305, 

-7.648) 

X4 

(-1.625, 

5.407) 

(  -1.757, 

0.820) 

x5 

(-1.625, 

-5.407) 

(  -1.757, 

-0.820) 

X6 

(-1.000, 

3.500) 

(  -1.000, 

3.000) 

(-1.000, 

-3.500) 

(  -1.000, 

-3.000) 

(-0.055, 

0.000) 

(  -0.061, 

0.000) 

(-4.430, 

2.070) 

(-17.808, 

0.000) 

(-4.430, 

-2.070) 

(  -3.863, 

0.000) 
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[ 


* 


r 

i 


Q 


v= Spectrum  of  Ar 
a = Retained  Optimal  Eigenvalues 
a=Unretained  Optimal  Eigenvalues 
o=Open  Loop  Spectrum  of  H0 

fP-«  336 


Figure  4.1.3.  Spectrum  of  compensated  Saturn  V  booster  under  retention 
of 
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Table  4.1.8.  Data  for  pole-placement  problem  based  on  ordering 


f1  =  [  1.000] 


9.781' 

33.041. 


gx  =  [4.766] 

g2  =  [1.000] 


0.000  14.548 

0.000  33.041 


The 

spectrum  of 

is : 

The 

spectrum 

o 

> 

u> 

II 

> 

1. 

(  -9.739, 

0.000) 

1. 

(-18.750, 

0.000) 

* 

2. 

(  -1.000, 

0.000) 

2. 

(  -2.908, 

0.000) 

3. 

(  -0.305, 

3.496) 

3. 

(  -1.000, 

3.000) 

4. 

(  -0.305,  - 

-3.496) 

4. 

(  -1.000, 

-3.000) 

5. 

(  -0.060, 

0.000) 

5. 

(  -0.061, 

0.000) 

I 


B 

o 


-0.027  0.005 

1.852  -0.392 

-4.468  2.197 

-0.010  -0.100 

1.819  0.181 


-0.062 

0.069 

-0.015 

-0.307 

0.000 

1.515 

-4.925 

0.985 

19.701 

0.000 

= 

-19.162 

17.276 

0.062 

5.894 

0.000 

0.271 

-0.879 

-0.003 

3.517 

1.000 

1.346 

-4.375 

-0.013 

-32.499 

-10.000 

"  -0.083 

0.135 

-0.015 

-0.574 

0.000 

2.547 

-8.277 

0.975 

33.107 

0.000 

SB 

-13.386 

-1.494 

0.004 

80.976 

0.000 

-0.413 

1.341 

0.004 

-5.363 

1.000 

_  6.102 

-19.831 

-0.061 

29.324 

-10.000 

1 

I 


I 

I 
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A-A.A.A,.  Recall  that  in  comparing  the  optimal  state  feedback  regulator 
with  the  minimum  energy  solution,  it  was  noted  that  the  pair  A^^  waa 
slightly  perturbed  from  -5+ j5  whereas  the  pair  A^A^  arose  from  two  real 
poles  near  -0.45  forming  a  complex  pair  and  leaving  the  real  axis.  It 
may  be  inferred  that  this  has  been  reflected  in  the  second  design  by  a 
large  expenditure  of  energy  to  retain  the  pair  A,^.  Apparently  the  first 
design  is  preferable. 

It  should  also  be  noted  that  in  all  of  the  static  and  dynamic 
designs,  including  the  unstable  designs  based  on  a  zero  order  compensator, 
there  is  a  real  pole  near  the  invariant  zero  at  -0.0462.  This  is  con¬ 
sistent  with  the  fact  that  under  high  gain  output  feedback  a  number  of  the 
poles  of  the  closed  loop  system  tend  to  the  finite  invariant  zeros  [36]. 

It  is  interesting  to  note  that  the  system  appears  to  have  two  widely 
separated  time  scales  in  the  sense  that  the  gains  of  interest  are  "high 
gain"  relative  to  the  zero  at  -0.0462,  but  not  relative  to  the  other  two 
zeros,  which  do  not  have  any  easily  discemable  influence  on  the  system. 

In  summary  the  second  order  compensator  design  obtained  for  the 
ordering  ^^^5^6  comPares  favorably  in  dimension  with  the  order  of  the 
reduced  order  observer  Cn-r=5)  and  with  the  bound  on  the  dimension 
required  for  arbitrary  placement  of  the  spectrum  of  (n-r-Jl+1  *  5) . 

The  compensator  is  open  loop  stable,  retains  a  four  dimensional  invariant 
subspace  of  the  optimal  state  feedback  regulator  in  which  there  is  no  cost 
degradation,  requires  modest  gains,  and  achieves  a  damping  ratio  of 
0.275  (n  -74°). 


87 


4.2.  Fifth  Order  Example 

This  is  a  trivial  example  illustrating  the  application  of  the 
design  methodology  to  stabilize  systems.  Let  the  system  be  given  by 


* 

X  : 


0 

1  1  0 

0 

ol 

r°i 

-6 

-5  1  0 

0 

0 

0 

0 

0  1  0 

1 

0 

X 

+ 

0 

0 

0  1  0 

0 

1 

0 

L  1 

0  1  0 

0 

-1 

Li 

the  permutation  of  states 

r° 

1 

0 

0 

0 

0 

0 

1 

0 

0 

T  * 

1 

0 

0 

0 

0 

0 

0 

0 

1 

0 

_0 

0 

0 

0 

1 

fo  ljo  0  ol. 

Lo  oil  o  oj*  <4-2-1> 


(4.2.2) 


to  obtain  the  system  x=Ax+Bu,  y  =  Cx  where 


-5 

0  -6 

0 

0* 

‘o' 

0 

0  0 

1 

0 

0 

1 

0  0 

0 

0 

B  = 

0 

0 

0  0 

0 

1 

0 

_  0 

0  1 

0 

-1 

.1. 

c-[i  l  2  S  §]•  w-2-« 


In  this  basis  x^,x^  are  the  state  variables  of  the  uncontrollable  subsystem 
with  eigenvalues  at  -2,-3,  while  x2,x^,x^  are  the  state  variables  of  the 
controllable  subsystem  with  eigenvalues  at  -1,0,0. 

The  solution  of  the  state  feedback  regulator  problem  with  11=1  and 
Q»  dg(l,5,0,2,0)  is  given  in  Tables  4.2.1  and  4.2.2.  It  is  noted  that  the 
optimal  controllable  eigenvectors  have  the  form  given  in  (3.3.16).  To 
insure  the  invertibility  of  Y  at  most  one  controllable  eigenvector  may  be 
retained  in  Ar»  and  the  possible  choices  are  (o^,ct2) , (a^,X^) , (a^.A^) . 

The  spectra  of  the  resulting  matrices  are  given  in  Table  4.2.2.  It  is 
noted  that  the  controllable  subspectrum  of  is  the  same  in  both  the  cases 
“  dg(Oj^)  and  ■  dg(a2^)  as  predicted  by  the  form  of  the  block  diagonal 


entries  of  in  (3.3.18) 
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Table  4.2.1.  Solution  of  state  feedback  regulator  for  fifth 
order  example 


0.1146 

0.0518 

0.0817 

0.1826 

0.1309 

0.0518 

8.7407 

0.5519 

6.6400 

2.2361 

- 

0.0817 

0.5519 

1.0782 

1.3728 

0.8610 

c 

0.1826 

6.6400 

1.3728 

9.3716 

3.9090 

_  0.1309 

2.2361 

0.8610 

3.9090 

1.9695_ 

K  = 

[  0.1309 

2.2361 

0.8610 

3.9090 

1.9695] 

"-5. 0000 

0.0000 

-6.0000 

0.0000 

o.oooo' 

0.0000 

0.0000 

0.0000 

1.0000 

0.0000 

F  * 

1.0000 

0.0000 

0.0000 

0.0000 

0.0000 

0.0000 

0.0000 

0.0000 

0.0000 

1.0000 

-0.1309 

-2.2361 

0.1390 

-3.9090 

-2.9695 

The  optimal  closed  loop  eigenvectors  are  (see  Table  4.2.2) 


"  0.9361 

0.8057' 

0.0170 

0.0948 

-0.3120 

v  ■ 

-0.4028 

-0.0510 

a2 

-0.1896 

L  0.1530 J 

L  0.3791 

’  o.oooo' 

0.4332 

o.oooo' 

0.3711 

'  0.0000' 
0.1866 

0.0000 

V  ■ 

A  X 

0.0000 

±j 

0.0000 

-0.5547 

2  3 

-0.5031 

0.2196 

L  0.7104 

L  0.2017 J 

L-0.6968 

* 


Table  4.2.2.  Eigenvalues  of  various  matrices  of  the  fifth 
order  example 


Spectrum  of  A 


L.  (-3.000, 

0.000) 

>.  (-2.000, 

0.000) 

3.  (-1.000, 

0.000) 

(  0.000, 

0.000) 

>.  (  0.000, 

0.000) 

Spectrum  of  Ac 


Spectrum  of  F 


(-3.000,  i 
(-2.000,  i 
(-1.281,  i 
(-0.844,  : 

(-0.844,  -; 


Spectrum  of  H0 


.  (-3.000, 

0.000) 

1.  (-2.485,  2. 

.  (-2.000, 

0.000) 

2.  (-2.485,  -2. 

1.  (-1.500, 

1.500) 

.  (-1.500, 

-1.500) 

.  (-1.281, 

0.000) 

».  (-0.844, 

1.016) 

(-0.844, 

-1.016) 

Spectra  of  possible  matrices  A^ 


Eigenvalues  retained  in  A 


(Vi> 


<Vi> 


Spectrum  of  Aj 


(-0.702,  0.274) 

(-0.702,  -0.274) 
(  0.404,  0.000) 

(-2.000,  0.000) 

(  0.140,  0.583) 

(  0.140,  -0.583) 

(-3.000,  0.000) 

(  0.140,  0.583) 

(  0.140,  -0.583) 
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Since  all  three  static  compensators  are  unstable,  a  dynamic  com¬ 
pensator  is  designed.  Consistent  with  the  desire  to  retain  as  many  control¬ 
lable  optimal  eigenvectors  as  possible  the  design  is  based  on  the  ordering 

^°2^1^2^3al^ '  (^e  or^er^n8  (fT2^i^2^3a2^  wou^  serve  equally  well.)  A 

second  order  compensator  will  be  required  since  (A^.A^)  is  a  complex  pair. 

The  reduced  pole-placement  problem  of  (3.3.19)  is  solved  in  two 
stages.  Since  the  final  closed  loop  spectrum  will  retain  the  entire  optimal 
spectrum  of  F,  the  choice  of  poles  to  assign  is  arbitrary  and  taken  to  be 
-1.5  +  J1.5.  At  the  first  stage  a  pole  is  arbitrarily  placed  at  0.0  using 
Procedure  1,  and  at  the  second  stage  the  complex  pair  is  assigned  using 
Procedure  2.  The  data  for  this  pole-placement  problem  are  given  in  Table 
A. 2. 3.  In  the  notation  of  (3. 3.21) , (3. 3.22)  the  parameters  of  the  resultant 


compensator  are: 


-0.3698 

1.2551" 

cT-0. 23321 

c  r-1.279ll 

.-8.4828 

-4.5997. 

o  L  0.0830 J 

o  L 20. 1197 J 

'  0  * 

r  0.0000 

0.0000 ' 

- 

-1.2453 

-0.1537 

NC 

-0.7261 

-1.1634 

L  p. 

S21  Br 


-0.5000 

0.0000  ‘ 

-0.0846 

~2.0o00 

0.2777 

2.5000 

K^-  [-6.2978  -2.8921] 

Ky-  [-0.0835]  Ky  *  [14.9777] .  (4.2.4) 

In  the  original  basis  the  total  closed  loop  matrix  is 


-0.3698 

1.2551 

o  ; 

-0.2332 

-1.2791  ! 

0 

0 

-8.4828 

-4.5997 

0  ! 

0.0830 

20.1197  ! 

0 

0 

0 

0 

0 

1 

0 

0 

0 

0 

0 

-6 

-5 

0 

0 

0 

0 

0 

0 

0 

0 

1 

0 

0 

0 

0 

0 

1 

.  6.2978 

2.8921 

0.0835 

-14.9777 

0 

-1 

(4.2.5) 
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Table  4.2.3.  Data  for  pole-placement  problem  for  fifth  order  example 


Data  for  full  problem  (see  equation  (3.3.19)): 


[-0.5000 

-0.00001 

[  0.0000 

0.0000 

N  * 
o 

-0.0846 
_  0.2777 

-1.2805 

1.6398. 

B  -1-0.0279 

0  [-0.4069 

0.4586 

-1.0028. 

[-3.0000 

0.0000 

o.ooool 

f  0.0000 

L  0.0000 

17.1739 

-6.1098. 

V 

-0.5078 
.  2.6660 

1.2805 

-1.6398 

1.0000 

-1.0000 

Data  for  reduced  problem: 

£1  -  [  1.0000] 


18.0003 

-6.1098 


- 


gx  -  [-0.826] 

g2  -  [1.000] 

[  1.3036 

1.0000  ] 

-2.0000 

1.0000  ] 

A2  *  L-1.3036 

-1.0000  J 

V. 

-2.5000 

-1.0000J 

The  spectrum  of  A2  is:  The  spectrum  of  is: 

1.  (0.000,  0.000) 

2.  (0.304,  0.000) 


1.  (-1.500,  1.500) 

2.  (-1.500,  -1.500) 


This  compensator  retains  all  the  optimal  eigenvalues  of  the  state  regulator 
and  is  open  loop  stable. 

4.3.  Twelfth  Order  Nuclear  Reactor  Model 

In  this  example  a  single  input,  three  output,  twelfth  order 
model  of  a  nuclear  reactor  is  considered.  In  a  preliminary  analysis  the 
system  structure  of  the  model  will  be  discussed  and  the  output  feedback 
pole-placement  problem  solved,  and  then  a  linear  quadratic  regulator  problem 
will  be  defined  on  the  basis  of  which  both  a  static  and  a  first  order 
dynamic  compensator  design  will  be  given. 

The  model  is  taken  from  [24]  and  referring  to  Table  4.3.1  is 

given  by 


The  system  consists  of  a  seventh  order  uncontrollable  subsystem  driving 
a  fifth  order  controllable  subsystem  The  eigenvalues  of  these  two 

systems  will  be  denoted  {o  and  {*^1=1  respectively  and  are  given  in 

Table  4.3.2  and  Figure  4.3.1.  The  observability  structure  of  the  system 
may  be  obtained  by  inspection  of  the  eigenvectors  given  in  Table  4.3.3: 

a)  The  uncontrollable  eigenvalues  (0^,0^)  are  also  unobservable. 

b)  The  first  row  of  C  observes  a  fifth  order  uncontrollable  subsystem 

with  eigenvalues  (o^  5,a7^’ 

c)  The  second  row  of  C  observes  a  first  order  controllable  subsystem 
with  a  zero  eigenvalue. 


Table  4.3.1.  System  matrices  for  the  nuclear  reactor  model 


/*,,  0 

0 

c  0 

The  mode] 

is  x  * 

lx  +  [ 

u.  y 

“  1 

\  A_.  A-- 

\0  c„ 

\  21  22 

\2/ 

\  2 

"-0.4044 

0.000 

0.000 

0.4044 

0.000 

0.000 

0.000 

-0.4044 

0.000 

0.000 

0.4044 

0.000 

0.000 

0.000 

-0.4044 

0.000 

0.000 

0.4044 

0.01818 

0.000 

0.000 

-0.5363 

0.000 

0.000 

0.000 

0.0818 

0.000 

0.4545 

-0.5363 

0.000 

0.000 

0.000 

0.0818 

0.000 

0.4545 

-0.5363 

0.000 

0.000 

0.000 

0.000 

0.150 

0.000 

’  0.000 

0.000 

0.000 

0.000 

0.000 

0.000 

0.000 

-7.500 

0.000 

0.000 

75.000 

0.000 

0.000 

0.000 

0.000 

0.000 

0.000 

0.000 

0.000 

0.000 

0.000 

0.000 

0.000 

0.000 

.  0.000 

0.000 

0.000 

0.000 

0.000 

0.000 

”  0.000 

0.000 

0.000 

0.000 

o.ooo' 

600.000  - 

74.995 

0.033 

0.346 

0.621 

0.000 

2.475 

-0.033 

0.000 

0.000 

0.000 

25.950 

0.000 

-0.346 

0.000 

.  0.000 

46.570 

0.000 

0.000 

-0.621 

c^  [0000001] 


r 10000 

C2  *  LOIOOO  . 


where: 


0.000' 

0.000 

0.000 

0.4545 

0.000 

0.000 

-0.150 


0.000 

0.000 

0.000 

0.000 

0.000 
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Table  4.3.2.  Open  and  closed  loop  eigenvalues  of  the  nuclear  reactor  model 


The  uncontrollable  and  controllable  eigenvalues  of  A  are  the 
eigenvalues  of  A-q  and  A22  respectively.  The  eigenvalues  of 
the  optimal  state  feedback  regulator  are  those  of  A^  and  F22* 
The  eigenvalues  of  the  final  compensated  system  are  those  of 
An  together  with  the  controllable  eigenvalues  of  Ac. 


Spectrum  of  A^  (a^) 

Spectrum  of  (X^) 

Spectrum  of  F 

22  CV 

(-0.664,  0.000) 

BK 

(-75.502,  0.000) 

BS 

(-75.197, 

0.000) 

(-0.631,  0.195) 

MM 

(  -0.446,  0.000) 

21 

(-13.051, 

12.119) 

(-0.631,  -0.195) 

82 

(-0.0476,  0.000) 

Si 

(-13.051, 

-12.119) 

(-0.379,  0.034) 

MM 

(  0.000  ,  0.000) 

Si 

(  -0.399, 

0.000) 

(-0.379,  -0.034) 

mm 

(  0.000  ,  0.000) 

Si 

(-0.034, 

0.000) 

(-0.277,  0.000) 

(-0.0112,  0.000) 

1 

1 

Controllable  eigen- 

Open  loop  spectrum 

values  of  Ac 

of 

V 

1. 

(-75.197,  0.000) 

(-76.132,  0) 

2. 

(-13.051,  12.119) 

3. 

(-13.051,-12.119) 

4. 

(  -7.000,  0.000) 

5. 

(  -0.408,  0.000) 

■ 

6. 

(  -0.034,  0.000) 

| 
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o  Uncontrollable  Eigenvalues  of  A 
oControllable  Eigenvalues  of  A 

□  Controllable,  Optimal  Eigenvalues  Not 
Retained  in  o-(Ac) 

a  Controllable,  Optimal  Eigenvalues 
Retained  in  <r(Ac) 

<>  Remaining  Controllable  Eigenvalues  of  Ac 
(the  Controllable  Eigenvalues  of  Ar) 

v  Open  Loop  Eigenvalues  of  the  Compensator  H0 

*P-«537 


Figure  4.3.1.  Eigenvalues  of  the  nuclear  reactor  model 


0.0000 

-0.8146 

0.0267 

0.2813 

0.5066 


0.0000 

0.0027 

-0.0161 

-0.7003 

0.7137 


0.0000 

0.0048 

-0.8179 

0.4205 

0.3927 


0.0000 

-0.0077 

-0.5774 

-0.5774 

-0.5773 


The  eigenvectors  corresponding  to  uncontrollable  eigenvalues 


of  are: 


-0.0000" 

0.0000 

-0.8417 

0.0000 

-0.0000 

0.5399 

-0.0000 

0.0000 

-0.0000 

0.0000 

-0.0000 

o.oooo. 


0.0071' 

0.0029 

-0.0210 

-0.0049 

0.0043 

0.0137 

-0.0003 

0.0000 

0.0039 

-0.0139 

-0.2167 

0.1087 


0.0020 

-0.0123 

-0.0041 

0.0023 

0.0082 

-0.0078 

-0.0027 

0.0000 

0.0007 

-0.0072 

-0.2075 

.-0.9470- 


0.0550 

-0.0398 

-0.0001 

0.0072 

-0.0009 

-0.0023 

0.0001 

0.0000 

-0.0003 

0.0008 

-0.5548 

-0.0966 


"-0.0453 

-0.0189 

0.0272 

0.0018 

-0.0045 

0.0017 

0.0030 

0.0000 

-0.0017 

0.0122 

0.7594 

-0.3133. 


-0.0000 

-0.0000 

0.9537 

0.0000 

0.0000 

0.3007 

0.0000 

0.0000 

0.0000 

-0.0000 

0.0000 

.-0.0000. 


0.0003 

0.0003 

0.0003 

0.0003 

0.0003 

0.0003 

0.0003 

0.0000 

-0.0064 

-0.7223 

-0.4926 

-0.4853 
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d)  The  third  row  of  C  observes  a  tenth  order  system  consisting  of  two 
fifth  order  controllable  and  uncontrollable  subsystems  with  eigen¬ 
values  -j»  ^4**5)  and  (ct2  3»°4  5,a7^  respectively. 

Since  output  feedback  only  affects  the  eigenvalues  of  the  con¬ 
trollable  and  observable  subsystem,  a  preliminary  analysis  of  the  pole- 
placement  problem  under  static  output  feedback  compensation  may  be  carried 

out  on  the  controllable  and  observable  triple  under  feedback 

3x1 

u=-(k2k3)C2*  (Here  K=  (k^  ^2  ,k3)£R  .)  When  the  linear  quadratic 

regulator  problem  is  defined,  however,  it  will  not  be  possible  to  exploit 
the  uncontrollability  of  to  reduce  the  dimensionality  of  the  system 
since  uncontrollable  eigenvectors  are  not  invariant  under  feedback,  and  are 
shaped  in  accordance  with  the  selected  cost  structure  (Q,R).  It  will  then 
be  advantageous  to  retain  the  feedback  gain  k^  through  C^. 

The  matrices  of  the  subsystem  (A22»B2’^2^  are  g*ven  by: 


22 


0 

0 

0 

0 

0  " 

'l' 

Y 

-(81+e2+e3) 

°l 

“2 

a3 

0 

0 

6i 

- 

’“I 

0 

0 

V 

0 

0 

CM 

ca 

0 

"a2 

0 

0 

0 

S3 

0 

0 

“a3_ 

.0. 

al  * 

0.033 

61 

=  2. 

475 

a2 

0.346 

62 

-  25. 

95 

a3  “ 

0.621 

®3 

=  46. 

57 

a 

l+a2+a3  “ 

1 

Y 

=  600 

.0 

0  0  0  0 
1  0  0  0. 


(4.3.2) 


and  have  the  controllability  canonic  form  : 
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1 

0 

0 

0 

0 


0  a2  a3  a4  1 

5  e  y  y  0 


6  »  4.2543 

e  =  148.066 

a.  -  75.995 

4 

a3  =  37.2614 

a2  =  1.60241 


(4.3.3) 


from  which  the  open  loop  characteristic  equation  is 


/i\  , 5  ,  .4  ,  , 3  ,  .2 

Pq(A)  =  A  +  a^A  +  a3A  +  a^A  . 


(4.3.4) 


Under  output  feedback  u=-(k2  k3)C2  c^e  closed  loop  characteristic  equation 


is 


Pc(A)  =  A5  +  (k2+a4)A4  +  (k2a4+Yk3+a3) A3  +  (k2a3+k3Y+a2) A2 


+  (k2a2+k3e)A  +  k35. 


(4.3.5) 

Before  investigating  the  root  locus  of  this  equation  it  is  of 
interest  to  determine  if  the  system  has  any  invariant  zeros.  Since  this 
is  the  controllable  and  observable  subsystem,  the  invariant  zeros  may  be 
determined  as  the  roots  of  the  greatest  common  divisor  of  all  the  minors  of 
full  order  of  the  system  matrix 


'^22_-X_I5_ 

-v 

-  °2” 

0 

(4.3.6) 


7*6 


or  equivalently  in  this  case  as  the  intersection  of  the  invariant  zeros  of 
the  two  square  systems  [35] 

Ta^-AI  «  -B,l 

(4.3.7) 


r  a„„- ai,  ' 

-B_1 

A_„-aI  1 

-B,  1 

22  5  1 

2 

22  1 

2 

.  10000 

0  _ 

_01000  1 

0 
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The  zeros  of  Che  first  system  are  the  distinct  eigenvalues  of  ^2? 

( those  of  the  second  system  are  the  parameters:  (a^c^a^) . 
Since  these  two  sets  are  disjoint  it  follows  that  the  controllable  and 
observable  subsystem  has  no  invariant  zeros. 

As  the  closed  loop  characteristic  equation  is  a  function  of  the 
two  parameters  a  convent*ona^  root  locus  plot  may  be  replaced  with 

a  graph  of  those  values  of  gain  (l^k^)  for  which  the  eigenvalues  of 


^22~®2^k2k3^2  *n  Prescr*bed  regions  of  the  complex  plane.  This  will 
be  done  here  for  the  stability  region  (the  left  half  plane)  and  for  the 
region  to  the  left  of  the  hyperbola  a*  /ai^  +  having  as  assymptopes 
the  constant  damping  lines  a  =  +oi. 

As  the  number  of  assignable  eigenvalues  is  two,  these  may 
momentarily  be  denoted  -s  and  -t  and  the  closed  loop  characteristic  equation 
written 


Pc(A)  -  (A+s) (A+t) (A34h  AZ*Hi  A-Hi  ) . 

Equating  (4.3.5)  and  (4.3.8)  gives  the  linear  system  of  equations 


(4.3.8) 


£]■ 


-6 


-a2  -« 
-a3  -Y 

-4  -Y 

-1  0 


st  0  0 

s+t  st  0 

1  s+t  st 
0  1 

0  0  1 


V 

0 

k3 

0 

'h’o" 

- 

”a2 

t 

hl 

a3~st 

J 

CM 

X 

_ 1 

a^-(s+t) 

(4.3.9) 


which  has  a  solution  for  all  values  of  s,t  for  which  |x|  #0.  (Compare 
equation  (1.35).)  The  boundary  of  the  region  of  stabilizing  gains  may  be 


found  as  the  parameterized  solution  (k^jk^)  *  ^(Sjt)  ,k.j(s,t))  under  the 


constraints 
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a)  a  real  pole  crosses  the  imaginary  axis:  s=0,  t>0 

2 

b)  a  complex  pair  crosses  the  imaginary  axis:  s+t=0,  st=w  ,  w*:0. 
Similarly  the  boundary  of  the  region  to  the  left  of  the  hyperbola 

o*  / uj^+a^  (a=  .030,  X  = -a+joi)  is  the  solution  under  the  constraints 

a)  a  real  pole  crosses  the  line  a  =  0.03:  s  =  0.03,  t>0.03 

2  2 

b)  a  complex  pair  crosses  the  hyperbola:  s+t=2o,  st  =  a  +w  ,  a>.03. 
A  graph  of  these  two  regions  is  shown  in  Figure  4.3.2,  and  suggests  that 

a  static  design  for  the  output  feedback  regulator  problem  may  be  found 
satisfactory.  It  should  be  noted  that  placing  the  eigenvalues  of  A 22  to 
the  left  of  the  hyperbola  has  required  the  gain  to  be  roughly  an  order 
of  magnitude  larger  than  the  gain  k^  Recalling  that  the  third  row  of  C 
observed  the  entire  subsystem  A^2  whereas  the  second  row  only  observed  a 
one  dimensional  generalized  eigenvector  of  ,-=0,  this  suggests  that  in 
designing  a  static  compensator  for  the  linear  quadratic  regulator  the 
observability  properties  of  the  system  may  be  reflected  in  disparate  values 
for  the  gains  k^  and  k^. 

A  linear  quadratic  regulator  problem  is  now  defined  with 

'o  0 

Q  *  R-l  Q22  -  dg[l,l,o1,o2,a3J.  (4.3.10) 

°  q22 

Only  the  states  of  the  controllable  subsystem  have  been  penalized,  but 
since  A^^Oj  the  optimal  solution  shapes  the  eigenvectors  of  the  uncon¬ 
trollable  eigenvalues  and  it  is  not  possible  to  exploit  the  uncontrol¬ 
lability  of  A^  to  reduce  the  regulator  problem  to  one  for  pair  (A22,B2). 

* 

It  is  true  that  there  is  a  two  dimensional  subsystem  with  eigen¬ 
values  (a^»ag)  which  could  be  removed  to  reduce  the  dimensionality  of  the 
problem  from  twelve  to  ten,  but  this  does  not  seem  worthwhile. 
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A  Real  Pole  Crossing  the  Imaginary  Axis 


A  Complex  Pair  Crossing  the  Imaginary  Axis _ 

A  Real  Pole  Crossing  the  Hyperbola  a-y  w2+(.03)2 
A  Complex  Pair  Crossing  the  Hyperbola  cr=Vw2+(.03)2 


Suboptimal  Pair  Based  on  Retention  of  the  Complex 
Pair  X2,3  of  the  Optimal  Feedback  Regulator 

'P-65S9 


Figure  4.3.2.  Stability  regions  for  the  nuclear  reactor  model 
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This  uncontrollability  will  prove  useful  however,  in  analyzing  the  pole- 
placement  problem  for  the  triple  (A^.B^A^). 

The  solution  to  the  regulator  problem  is  given  in  Table  4.3.4, 
and  the  resultant  optimal  eigenvectors  are  given  in  Table  4.3.5. 
Qualitatively,  the  optimal  solution  has  moved  the  double  pole  at  the  origin 
to  the  complex  pair  -13+ jl2  while  only  slightly  perturbing  the  remaining 
three  real  poles,  suggesting  that  an  output  feedback  design  should  concen¬ 
trate  on  retaining  this  complex  pair.  Examination  of  the  eigenvectors  of 
F  in  Table  4.3.5  shows  that  there  are  severe  restrictions  on  the  permissible 
choices  of  A^.  Since  must  satisfy  rank(CV^  ) =  3  the  choices  may  include 
at  most  two  eigenvalues  of  F22  and  are:  (Xi02a3^  ’  ^ia4a5^  ’  ^2X3a7^ 

(X^a2a3> .(X^o^a^) ,(X^  A^  o^) .  (Those  involving  X^  are  omitted  as  the 
eigenvector  corresponding  to  X^  is  nearly  in  the  null  space  of  C.)  Since 
it  is  desired  to  retain  as  many  optimal  eigenvectors  as  possible, a  design 
is  based  on  retention  of  with  the  expectation  that  if  a  compen¬ 

sator  is  used,  additionally  one  or  more  of  A^.A^.A^  may  be  retained. 

Based  on  this  choice  of  Ar  and  under  the  permutation  of  states 


0 

I, 


0 

0 


3  J 


(4.3.11) 


bringing  C  to  the  form  [1  0],  the  resulting  matrices  A, ,N  ,B  ,A,„  are 

l  o  o  iz 

given  in  Table  4.3.6.  The  spectrum  of  A^  consists  of  the  uncontrollable 
eigenvalues  together  with  three  controllable  eigenvalues 

located  at  -6.066,  -0.407,  -0.034.  Thus  the  system  is  stable  as  may  have 
been  anticipated  from  the  preliminary  output  feedback  analysis  and  the 
corresponding  gain  K“  (-4.502,-43.385,6.249)  lies  within  the  region  of 
Figure  4.3.2. 
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Table  4.3.4.  Optimal  state  feedback  regulator  solution 


0.0000 

0.0002 

0  0.0008  0.0009 

0 

0 

.0024 

0.0002 

0.0047 

0  0.0034  -0.0348 

0 

0 

.0115 

0.0000 

-0.0000 

0  -0.0000  0.0000 

0 

0 

.0000 

0.0008 

0.0034 

0  0.0214  0.0281 

0 

0 

.0521 

0.0009 

-0.0348 

0  0.0281  0.4209 

0 

0 

.0552 

0.0000 

-0.0000 

0  -0.0000  0.0000 

0 

0 

.0000 

0.0024 

0.0115 

0  0.0521  0.0552 

0 

0 

.1580 

0.0000 

-0.3064 

0  0.0568  3.1448 

0 

0 

.0001 

-0.0000 

-0.0069 

0  0.0001  0.0690 

0 

-0 

.0000 

0.0000 

-0.0002 

0  0.0000  0.0024 

0 

0 

.0001 

-0.0000 

-0.0037 

0  0.0001  0.0373 

0 

-0 

.0000 

.-0.0000 

-0.0084 

0  0.0001  0.0846 

0 

-0 

.0001 

0.0000 

-0.0000  0.0000  -0 

.0000 

-o.oooo- 

-0 . 3064 

-0.0069  -0.0002  -0 

.0037 

-0.0084 

0.0000 

0.0000  0.0000  0 

.0000 

0.0000 

0.0568 

0.0001  0.0000  0 

.0001 

0.0001 

3.1448 

0.0690  0.0024  0 

.0373 

0.0846 

0.0000 

0.0000  0.0000  0 

.0000 

0.0000 

0.0001 

-0.0000  0.0001  -0 

.0000 

-0.0001 

25.7364 

0.5511  0.0191  0 

.2980 

0.6757 

0.5511 

0.0271  0.0008  0 

.0126 

0.0291 

0.0191 

0.0008  0.4953  -0 

.0132 

-0.0175 

0.2980 

0.0126  -0.0132  0 

.3843 

-0.1897 

0.6757 

0.0291  -0.0175  -0 

.1897 

0.1615 

-  [0  -0.3064  0  0.0568  3.1448  0  0.000l| 

25.7364  0.5511  0.0191  0.2980  0.6757] 
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Table  4.3.5.  Optimal  eigenvectors  of  state  feedback  regulator  solution 


The  eigenvectors  corresponding  to  controllable  eigenvalues  of  F ^  are: 


"  0.0000 

■  0.0000' 

0.0000 

■  0.0000 ' 

'  0.0000 

0.0000 

0.0000 

0.0000 

0.0000 

0.0000 

0.0000 

0.0000 

0.0000 

0.0000 

0.0000 

0.0000 

0.0000 

0.0000 

0.0000 

0.0000 

0.0000 

0.0000 

0.0000 

0.0000 

0.0000 

0.0000 

V  -  V  * 

\  9  X 

0.0000 

±j 

0.0000 

vx  * 

0.0000 

VX  = 

0.0000 

0.0000 

2  3 

0.0000 

0.0000 

4 

0.0000 

5 

0.0000 

0.0004 

0.0305 

-0.0141 

-0.0004 

-0.0001 

0.8135 

0.2477 

-0.1862 

-0.0019 

-0.0002 

-0.0268 

-0.0429 

-0.0045 

0.0126 

0.9997 

-0.2820 

-0.4549 

-0.0535 

0.9202 

-0.0186 

.-0.5080 

.-0.8246. 

-0.1062 

.-0.3912 

-0.0178 

The  eigenvectors  corresponding  to  uncontrollable  eigenvalues  of 


*11 


are: 


'  0.0000' 

-0.2172 

"  0.0707' 

-0.0000 

0.1347 

0.3648 

0.4817 

0.6117 

-0.2534 

-0.0000 

0.0876 

-0.1441 

0.0000 

-0.2509 

-0.1394 

-0.5399 

-0.2205 

0.4361 

0.0000 

0.0521 

±j 

0.0646 

-0.0000 

°2  °3 

0.0330 

0.0220 

0.0000 

0.0000 

-0.0000 

0.0000 

-0.0000 

0.0000 

0.0000 

-0.0000 

-0.0000 

_-0.0000. 

.-0.0001, 

-0.0001 

0.3313' 

"  0.7326' 

0.0000' 

'  0.3678' 

0.3228 

-0.3781 

0.0000 

0.3791 

-0.2965 

-0.0812 

0.9537 

0.3907 

-0.0405 

0.0731 

-0.0000 

0.3576 

0.0516 

0.0035 

0.0000 

0.3685 

-0.0116 

-0.0298 

0.3007 

0.3799 

-0.0327 

±j 

-0.0071 

-0.0000 

0.3984 

-0.0024 

-0.0052 

V  = 
O, 

-0.0000 

V  B 

-0.0413 

0.0000 

0.0000 

O 

-0.0000 

7 

0.0000 

0.0000 

-0.0000 

0.0000 

0.0000 

0.0000 

-0.0000 

0.0000 

0.0000 

.-0.0000. 

.  0.0000. 

--0.0000. 

.  0.0000. 

ffl 


Table  A. 3. 6.  Data  for  pole-placement  problem  for  the  nuclear  reactor 


model 

based  on 

retention 

of 

^2 » ^3  > 

VX1 

'o 

0  0 

0  0.15 

0  • 

0 

0 

0 

A,-  * 

0 

0  0 

0  0 

°J 

0 

0 

0 

_0 

'-7_.  5  0 

0“  75.0 

"0  i 

0.033 

0.346“ 

0.632 

"  0.9232  t 

0.0000 

0.0000' 

‘  0.0000' 

0.9516  i 

0.0000 

0.0000 

0.0000 

0.9808  , 

0.0000 

0.0000 

0.0000 

0.8976 

0.0000 

0.0000 

0.0000 

0.9252  ' 

0.0000 

0.0000 

B  = 

0.0000 

0.9536  l 

0.0000 

0.0000 

o 

0.0000 

-074328“, 

-4.1715“ 

0.3397 

“  0".  30f4 

-4.6563 

-44.8802 

3.6811 

3.2578 

.-8.5479 

-82.3894 

6.8002. 

.  6.0054. 

“-0.4044  0.0000 

0.0000 

0.4044 

-0.1385 

0.0000 

0.0000  -0.4044 

0.0000 

0.0000 

0.2617 

0.0000 

0.0000  0.0000 

-0.4044 

0.0000 

-0.1471 

0.4044 

0.0182  -0.0000 

0.0000 

-0.5363 

-0.1346 

0.0000 

0.0000  0.0818 

0.0000 

0.4545 

-0.6751 

0.0000 

0.0000  -0.0000 

0.0818 

0.0000 

0.3115 

-0.5363 

“0.0000  “  275478 

0.0000 

0  70000 

'  -25.4130 

0.0000 

0.0000  27.6085 

0.0000 

0.0000 

-275.3962 

0.0000 

.  0.0000  51.0014 

0.0000 

0.0000 

-508.7316 

0.0000 

;  o.oooo 

0.0000 

0.0000“ 

,  0.0000 

0.0000 

0.0000 

l  0.0000 

0.0000 

0.0000 

I  0.0000 

0.0000 

0.0000 

1  0.0000 

0.0000 

0.0000 

!  o.oooo 

0.0000 

0.0000 

,  -0.0442 

-07 1175 

'  -0 . 2lf0 

,  -0.1215 

-1.6197 

-2.2860 

»  -0.2244 

-2.3529 

-4.8439 
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Qualitatively  this  output  feedback  solution  retains  the  eigen¬ 
vectors  corresponding  to  the  optimal  complex  pair  A2,A3and  places  two 
eigenvalues  (though  not  their  eigenvectors)  very  near  their  optimal 
locations  A^.A^.  However  it  has  also  shifted  the  fast  eigenvalue  of  the 
system  from  -75.2  to  only  -6.066.  This  suggests  trying  to  design  a  first 
order  compensator  to  additionally  retain  A^  with  the  expectation  that  two 
eigenvalues  may  still  remain  near  A^A.^.  The  data  for  the  pole-placement 
problem  are  given  in  Table  4.3.6. 

Because  of  the  block  structure  of  (A^.B^.A^)  which  clearly 

satisfies  the  requirement  that  the  spectrum  of  A^  contain  the  unretained 

uncontrollable  eigenvalues  (o, ,cr„,a,,a. ,cr  ,a,) ,  and  the  fact  that  only  the 

third  row  of  A^2  observes  the  controllable  subsystem  of  A^,the  pole-placement 

problem  may  be  replaced  with  the  single-input  single-output  pole -placement 

22  2  22 

problem  for  the  triple  (A^  ,Bq>Ax2).  This  problem  may  then  be  solved 

analytically  by  transforming  to  the  controllability  canonic  form; 

0  1  O' 

0  0  1 

_-aQ  -ax  -a2 

aQ  *  0.083374  cQ  =  0.067455  (4.3.12) 

ax  -  2.6912  cx  *  2.1602 

a2  -  6.5080  c£  -  4.8663 


[c0  C1  C2] 


L 


•  i 


The  closed  loop  characteristic  equation  is 


Pc(A)  «  A3  +  (a2-c3p3)A2  +  (a1-c1p;J)A  +  (aQ-coP3)  , 


(4.3.13) 


This  gives  rise  to  the  root  locus 


i3  ,  .  , 2 


A  +  a2A  +  a^A  +  a^ 


(4.3.14) 


c2A  +c1A  +  c() 


2 
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where  the  roots  of  the  numerator  are  the  controllable  eigenvalues  of  A^  and 
the  roots  of  the  denominator  are  -0.0338,  -0.4101.  This  near  pole-zero 
cancellation  confirms  the  expectation  that  a  first  order  compensator  would 
allow  retention  of  and  yet  still  keep  two  eigenvalues  very  near  their 
optimal  locations  A^,A,.  without  requiring  large  gains.  Analysis  of  the 
root  locus  shows  that  the  system  is  stable  for  p^< 1.235,  with  two  poles  very 
nearly  at  -0.41,  -0.034  for  -®<p^<1.0.  Arbitrarily  placing  a  pole  at 
-7.0  gives  p^® -0.192  and  the  controllable  subspectrum  for  A^  becomes 
-7.0,  -0.408,  -0.034. 

The  parameters  of  the  compensator  are  then : 

Hq  -  [~76.132] 

T)c  -  [  -9.186]  *  [-88.540  -1.103] 

"  (4.3.15) 

K~  ®  [  5.820] 

Rc  -  [  -4.502]  KC  «  [-43.385  7.216] 

y  y 

and  in  the  original  basis  the  final  closed  loop  system  is  |xj*^c 

where  A  ■ 
c 


-76.132 

0.000 

0.000 

0.000 

0.000 

0.000 

0.000 

0.000 

-0.4044 

0.000 

0.000 

0.4044 

0.000 

0.000 

0.000 

0.000 

-0.4044 

0.000 

0.000 

0.4044 

0.000 

0.000 

0.000 

0.000 

-0.4044 

0.000 

0.000 

0.4044 

0.000 

0.01818 

0.000 

0.000 

-0.5363 

0.000 

0.000 

0.000 

0.000 

0.0818 

0.000 

0.4545 

-0.5363 

0.000 

0.000 

0.000 

0.000 

0.0818 

0.000 

0.4545 

-0.5363 

0.000 

0.000 

0.000 

0.000 

0.000 

0.150 

0.000 

-5.820 

0.000 

0.000 

0.000 

0.000 

0.000 

0.000 

0.000 

-7.500 

0.000 

0.000 

75.000 

0.000 

0.000 

0.000 

0.000 

0.000 

0.000 

0.000 

0.000 

0.000 

0.000 

0.000 

0.000 

0.000 

0.000 

0.000 

0.000 

0.000 

0.000 

0.000 

0.000 
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-9.186 

-88.540 

-1.103 

0.000 

0.000 

o.ooo' 

0.000 

0.000 

0.000 

0.000 

0.000 

0.000 

0.000 

0.000 

0.000 

0.000 

0.000  . 

0.000 

0.000 

0.000 

0.000 

0.000 

0.000 

0.000 

0.4545 

0.000 

0.000 

0.000 

0.000 

0.000 

0.000 

0.000 

0.000 

0.000 

0.000 

0.000 

0.000 

0.000 

0.000 

0.000 

0.000 

0.000 

-0.150 

0.000 

0.000 

0.000 

0.000 

0.000 

4.502 

43.385 

-7.216 

0.000 

0.000 

0.000 

0.000 

600.000 

-74.995 

0.033 

0.346 

0.621 

0.000 

0.000 

2.475 

-0.033 

0.000 

0.000 

0.000 

0.000 

25.950 

0.000 

-0 . 346 

0.000 

0.000 

0.000 

46.570 

0.000 

0.000 

-0.621. 

Thus  by  taking  into  account  the  eigenvector  structure  of  the 
optimal  state  feedback  regulator  solution,  the  methodology  for  designing 
regulators  has  been  applied  to  a  nontrivial  system  having  uncontrollable 
and  unobservable  modes.  A  first  order  compensator  has  been  constructed 
having  the  property  of  retaining  from  the  optimal  regulator  solution  a 
four  dimensional  '.nvariant  subspace  spanned  by  eigenvectors  corresponding 
to  three  controllable  and  one  uncontrollable  optimal  eigenvalues,  at  the 
same  time  placing  eigenvalues  near  all  the  remaining  optimal  locations. 


4.4.  Two  Interconnected  Power  System  Model 

This  example  considers  a  model  of  a  power  system  consisting  of 
two  interconnected  steam  generators.  The  model,  derived  in  [37],  is  given 
by  x*Ax+Bu,  y«Cx,  A,B,C  defined  in  Table  4.4.1,  and  represents  a 
linearization  of  the  system  about  an  operating  point,  describing  the  system 
behavior  under  real  power  and  frequency  variations.  The  state  vector  has 
the  physical  correspondence: 

x^,x^  -  valve  displacements  in  areas  one  and  two 
X2»xg  "  power  displacements  of  high  pressure  turbines  in  areas 
one  and  two 


16) 
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Table  4.4.1.  System  matrices  for  the  two  interconnected  power 
system  model 


The  system  is  represented  as: 
A  -a. 


0 

a. 


0 

-a„ 


s  J 


X  + 


0 

0 

b 


c  0 

0  1 

.0  0 


0 

0 

c  J 


where : 


-2.0 

4.75 

0 


0 

-5.0 

0.16667 


0 

0 

-0.16667 


£» 

0 

0 

2.0 

0 

0.025 

0.02333 

al  * 

t  0 

0 

0 

0 

0.08333] 

a2  * 

t  0 

0 

0 

0 

22.21439] 

b  - 

[4.0 

0 

0 

0 

o]T 

c  -  1 

[  0 

0 

0 

0 

1] 

0 

0 

0 

2.0 

0.035 


-4.0 

0 

0 

0 

-0.1125 


x^,Xg  -  power  displacements  of  intermediate  pressure  turbines  in 


areas  one  and  two 


X4,X10  "  P°wer  displacements  of  low  pressure  turbines  in  areas 
one  and  two 

X5,X11  ~  fre£luency  deviations  in  areas  one  the  two 

x,  -  tie-line  power  flow  deviation  from  area  one  into  area  two 
o 

and  the  controls  are  the  set  point  adjustments  in  the  two  areas.  It  is 
assumed  that  the  tie-line  power  flow  and  the  two  frequency  deviations  are 
available  for  measurement. 

Since  the  system  is  open  loop  stable  a  preliminary  analysis  of 

the  static  output  feedback  pole-placement  problem  may  be  replaced  by  an 

investigation  of  the  eigenstructure  of  both  the  open  and  closed  loop 

systems.  A  clear  understanding  of  this  structure  will  be  essential  in 

applying  the  design  methodology. 

The  eigenvalues  of  the  open  loop  matrix  may  be  classified 

according  to  whether  the  tie-line  power  flow  variable  x,  is  zero  in  the 

o 

corresponding  eigenspace.  The  eigenvector  equation  for  the  open  loop 


matrix  is: 


A  v  -  aa,  “  Xv 
s  1 

a2v  -  a2w  =  aX 

aa,  +  A  w  =  Xw 
1  s 


These  equations  may  be  rewritten: 


a(2a0  adj (XI-A  )a,  +  Xdet(XI-A  ))  =  0 
i  si  s 

(A  -XI) v  -  aa 
s  1 

(A  -XI)w  *  -aa,  . 
s'  1 


-  _.5*1  __1><1 

v,wt  R  ,  a£R 


(4.4.1) 


(4.4.2a) 

(4.4.2b) 

(4.4.2c) 


Ill 


If  a  *  0  then  there  are  five  solutions : 


{xM-r  <>i.r  *1 


■  d 
vi 


LviJ 


.  d  ,d  d 
A  v.  =  X  v. . 
si  i  1 


(4.4.3) 


If  a f  0  then  (4.4.2a)  is  a  sixth  degree  polynomial  in  X  and  the  solution  of 
(4.4.2)  is 


{Xc}6  { if\b  if 

1  i'i-1*  1  Vi-1*  i 


rvp 


L"Vi 


[X^!-As  :  aL] 


rvC'’ 

1 


a . 

u  x  -i 


(4.4,4) 


The  numerical  values  of  X?,X^,^,^,a^  are  given  in  Tables  4.4.2,  4.4.3. 

Thus .given  any  initial  condition,  the  transient  response  of  the 
system  may  be  represented  as  the  sum  of  a  coupled  and  a  decoupled  response 


x(t)  -  x  (t)  +  x,(t) 
c  d 


r  c  i 

r  «h 

6 

x  (t)  =  Z 
c  i-1 

vi 

a . 
i 
c 

c# 

V  ’ 

5 

x  (t)  =  I 
d  i=1 

vi 

0 

d 

-  ViJ 

-Vi- 

d  ^ 

V 


(4.4.5) 


xQ  -  x(0) 


i-1 


r  c  - 

\r 

r  di 

Vi 

5  d 

vi 

a. 

+  ly 

0 

i 

c 

i*l  1 

d 

LViJ 

-vi- 

The  first  term  x£(t)  represents  a  projection  of  the  system  response  into 

a  subspace  in  which  the  tie-line  power  flow  variable  is  nonzero  (a  ^0,  x,#0) 

i  o 

and  in  which  the  response  of  the  second  generator  is  in  exact  opposition 

to  that  of  the  first  generator.  The  second  term  x, (t)  represents  the  pro- 

a 

jection  of  the  response  into  a  subspace  in  which  the  tie-line  power  flow 

variable  is  zero  (x  =0)  and  the  two  generators  operate  synchronously.  For 
o 

example,  consider  the  case  that  an  initial  displacement  from  the  nominal 
operating  point  occurs  symmetrically  in  the  two  areas.  If  x. (0) =  x. , , (0) , 
i*l,...,5,  x^(0)“0,  then  the  two  coupled  steam  generators  evolve  with  the 
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Table  4.4.3.  Open  loop  eigenvectors  of  the  two  interconnected  power 
system  model 


c  d 

See  Table  4.4.2  for  the  definition  of  the  eigenvalues 


eigenvectors  have  the  form  V  = 

Kl 

ai 

d 

’  ^i  = 

-v4i 

1 

0 

c 

:viJ 

< 

W-  p. 

t  — 

0.0041 

0.2034" 

f  0.0463" 

'-0.0027*1 

-0.7062 

0.3222 

0.0620 

-0.0027 

0.0242 

-0.0292 

-0.0073 

-0.4507 

-0.0160 

-0.2381 

0.5329 

-0.4914 

0.0031 

0.0003 

-0.0054 

0.0012 

-0.0276 

*2*3" 

-0.0121 

±j 

0.1197 

-0.3326 

-0.0041 

-0.2034 

-0.0463 

0.0027 

0.7062 

-0.3222 

-0.0620 

0.0027 

-0.0242 

0.0292 

0.0073 

0.4507 

0.0160 

0.2380 

-0.5329 

0.4914 

L-0.003iJ 

L-0.0003. 

.  0.0054 J 

.-0.0012 

The 


c 

6 


*i*i 


-0.0389' 

0.0497 

*-0.0047* 

0.0975* 

0.2414* 

-0.0162 

0.0541 

0.7065 

0.1702 

0.3703 

0.0046 

0.0016 

-0.0242 

-0.0129 

-0.0353 

0.0032 

-0.0015 

0.0160 

-0.5093 

0.0730 

0.0429 

-0.0052 

-0.0036 

A  A 

0.0079 

-0.0053 

-0.1502 

±3 

-0.9794 

0.0000 

II 

3  0 

3  C* 

0.0000 

±1 

0.0000 

0.0389 

-0.0497 

-0.0047 

tm  -J 

0.0975 

0.2414 

0.0162 

-0.0541 

0.7065 

0.1702 

0.3703 

-0.0046 

-0.0016 

-0.0242 

-0.0129 

-0.0353 

-0.0032 

0.0015 

0.0160 

-0.5092 

0.0730 

.-0.0429. 

.  0.0052. 

.-0.0036. 

.  0.0079. 

.-0.0053. 

0.2434' 

'-0.2146  ' 

0.2319 

-0.2175 

-0.2193 

-0.2775 

-0.2603 

-0.2796 

-0.1203 

0.0900 

0.0000 

1  + 
U 

0.0000 

0.2434 

-0.2146 

-0.2319 

-0.2175 

-0.2193 

-0.2775 

-0.2603 

-0.2796 

. -0.1203 . 

-  0.0900. 

dynamics  of  a  single  fifth  order  system.  Note  that  the  total  system  is  not 

observable  from  the  tie-line  power  flow  variable.  If  x^O)  *  -xi+g(0)  * 

i=l,...,5,  then  the  system  evolves  in  a  six  dimensional  invariant  subspace, 

the  two  steam  generators  operating  in  exact  opposition  to  each  other. 

d  c 

Inspection  of  the  numerical  values  of  A_^,A^  in  Table  4.4.2  shows 
that  the  eigenvalues  of  A  may  also  be  identified  in  pairs: 
c  cl 

{X  }  -  {Xx}  two  real  eigenvalues  at  ~-5 

{ two  complex  pairs  at  ~  -2  +  jO.l 
c  c  d  d 

{A  , A&}  -  {A4, A5}  two  different  dominant  complex  pairs. 

The  eigenstructure  of  the  system  indicates  that  rather  than  representing  a 
pairing  of  distinct  modes,  one  from  each  of  the  two  steam  generators,  this 
symmetry  represents  the  pairing  of  coupled  and  decoupled  modes.  For  example  , 
the  frequency  pair  -0.15+j0.15  may  be  associated  with  decoupled  synchronous 
mechanical  rotation  of  the  two  generator  shafts,  while  the  pair  -0.06+j2.0 
may  be  associated  with  inversely  coupled  mechanical  rotation  of  the  two 
shafts.  If  the  initial  frequency  deviations  in  the  two  areas  were  equal 
then  the  response  would  be  well  damped,  while  if  the  initial  frequency 
deviation  in  one  area  were  the  negative  of  that  in  the  other  area,  the 
response  would  be  slow  and  highly  oscillatory. 

It  may  be  shown  that  the  closed  loop  system  possesses  the  same 
eigenstructure  symmetry  as  the  open  loop  system  provided  the  linear  quadratic 
regulator  problem  is  defined  such  that  the  states  of  the  generators  are 
weighted  equally:  Q  ■  dg(Q' ,q,Q' ) .  The  resultant  closed  loop  matrix  has  the 


symmetric  form 
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-fl 

0 

f. 


-a„ 


(4.4.6) 


and  the  closed  loop  eigenvectors  and  eigenvalues  are  given  by: 


(A?)5  , 

1  i=l 


U?}6  , 

i=l 


{^>5  , 

i*l 


fr'}6  , 

1  i=l 


V 


r  d 
vi 
°d 

vil 


rv 


i 

c 

i  J 


(F1+Vvi 


.d  d 

Vi 


rvC-’ 

1 


(4.4.7) 


Therefore  all  the  remarks  made  regarding  the  open  loop  transient  response 
apply  also  to  the  closed  loop  response. 

A  linear  quadratic  regulator  problem  may  now  be  defined  by 
letting  Q=  dg(5, 0,0, 0,30, 10, 5, 0,0, 0,30) ,  R«  ^2*2'  The  so^ution  to  tlie 
Riccati  equation  is  given  in  Table  4.4.4,  and  the  optimal  closed  loop  eigen¬ 
vectors  are  given  in  Table  4.4.5.  Referring  to  Figure  4.4.1,  the  optimal 
solution  has  increased  the  damping  of  the  two  complex  pairs  associated  with 
coupled  and  decoupled  mechanical  rotation  of  the  generator  shafts.  The  real 
eigenvalue  associated  with  the  tie-line  interaction  has  increased  in  magni¬ 
tude:  A  » -0.220.  The  remaining  spectrum  consists  of  three  pairs  of  real 
o 

eigenvalues  at  -2, -5, -9  corresponding  to  coupled  and  uncoupled  modes.  It  is 

c  d 

noted  that  the  closed  loop  eigenvectors  4.4.5)  corresponding 

to  the  real  eigenvalues  at  -5  are  nearly  equal  to  the  eigenvectors 
(Table  4.4.3)  corresponding  to  the  two  open  loop  eigenvalues  near  -5. 

Since  the  optimal  control  has  expanded  no  energy  shaping  these  two  eigen¬ 
vectors  it  is  concluded  that  these  two  modes  might  be  neglected  in  designing 
an  output  feedback  compensator. 
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Table  4.4.4.  Optimal  state  feedback  regulator  solution 


With  Q*  dg(5, 0,0, 0,30, 10, 5, 0,0, 0,30)  and  R=I?x2  the  solution  of 
the  Riccati  equation  and  the  associated  feedback  matrix  are: 


'M1 

~ml 

V 

~K 

-0.8761  k. 

16.9456 

-m 

K  « 

i 

z 

2. 

Z 

k- 

0.8761  k. 

_M2 

”1 

Ml_ 

L  2 

1 J 

where: 


"  0.4627 

0.0292 

0.1811 

0.0719 

4.3802' 

0.0292 

0.0600 

0.3147 

0.1426 

10.3402 

0.1811 

0.3147 

6.2726 

0.5895 

6.1776 

0.0719 

0.1426 

0.5895 

0.3795 

23.7609 

_  4.3802 

10.3403 

6.1773 

23.7608 

2478.8810, 

'-0.0112 

-0.0222 

-0.0439 

-0.0626 

-4.0168' 

-0.0222 

-0.0458 

-0.0492 

-0.1233 

-9.2460 

-0.0439 

-0.0492 

-0.9310 

-0.2711 

11.1443 

-0.0626 

-0.1233 

-0.2711 

-0.3484 

-21.9648 

„-4.0170 

-9.2462 

11.1440 

-21.9648 

-2367.8620, 

[-0.2190 

-0.2945 

-3.9340 

-1.2252 

25.3958] 

"  0.2190' 

0.2945 
3.9342 
1.2252 
.-25.3955. 


kx-  [  1.8507  0.1169  0.7244  0.2875 
k2=  [-0.0449  -0.0889  -0.1755  -0.2502 


17.5206] 

-16.0673] 
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Open  loop  eigenvalues 


A  2 


1 


— 6 


-5 


-3 


-1 


---1 

A  -2 


Optimal  closed  loop  eigenvalues 


-B — $ - p - - 1 - § 

-9  -5  -3  -2  -1 


a  Coupled  Modes 
o  Decoupled  Modes 


-l 


A± 


-2 


PP-6538 


Figure  4.4.1.  Open  and  closed  loop  spectra  of  the  power  system  model 
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To  apply  the  design  methodology  a  permutation  of  states  x”Tx  is 
introduced  where 

'0  0001000000" 

00000100000 

00000000001 

10000000000 

01000000000 

T  *  00100000000  .  (4.4.8) 

00010000000 
00000010000 
00000001000 
00000000100  . 

.0  0  0  0  0  0  0  0  0  1  0. 

In  this  basis  the  observation  matrix  has  the  form  (I_  _!0)  and  the  state 
vector  is  (x5 ,*fi ,xu .x^.x^x^.Xg.x^) . 

At  the  first  stage  of  the  design  it  is  necessary  to  select  three 

eigenvectors  such  that  Y  1  exists,  where  the  ith  column  of  Y  consists  of 

the  5th,  6th,  and  11th  entries  of  the  corresponding  selected  eigenvector. 

In  view  of  the  above  discussion  of  the  closed  loop  eigenstructure,  it 

T 

follows  that  the  columns  of  Y  have  the  form  (6,0,8)  if  a  decoupled  eigen- 

T 

value  is  selected,  and  (6,Y,~S)  if  a  coupled  eigenvalue  is  selected.  Thus, 

for  Y  to  be  invertible,  A  must  contain  exactly  one  decoupled  eigenvalue  and 

two  coupled  eigenvalues.  In  particular,  the  selection  of  the  three  dominant 

c  c  c 

coupled  eigenvalues  A^,A^,Ag  is  not  possible. 

Subject  to  complex  pairing  there  are  21  possible  choices  for  A^, 

of  which  four  give  a  stable  spectrum  for  A^:  (A^.A^.A^)  .(A^A^,!^) , 

c  c  d  c  c  d 

(A^, Aj, A2) , ( A^ , Aj, A^) .  For  each  of  these  choices  however,  the  spectrum  of 
A^  contains  at  least  one  eigenvalue  with  real  part  to  the  right  of  a* -0.060. 
Since  the  spectrum  of  the  open  loop  system  lies  to  the  lect  of  the  line 
o* -0.060,  a  static  output  feedback  design  is  not  acceptable  and  a  dynamic 
compensator  will  be  designed. 


To  solve  the  pole-placement  problem  of  the  second  stage  of  the 


design,  A^  must  be  selected  and  the  remaining  optimal  eigenvalues  ordered. 
Because  of  the  nature  of  the  tie-line  equation  the  second  row  of  is  all 
zeros  and  it  is  necessary  to  introduce 


~c  O' 

fi  0 

°] 

TA12  " 

1 

Lo  0 

ij  ’ 

-°  Cl- 

c±  -  [0.035  0.02333  0.025  0] 


(4.49) 


Then  l *  2  and  at  the  ith  stage  of  the  algorithm  only  i+Z-1  ■  i+1  eigenvalues 
may  be  assigned  to  A^+^.  It  is  assumed  that  the  dimension  of  the  desired 
compensator  may  not  exceed  four.  At  first  designs  based  on  stable  matrices 
A^  are  discussed  and  then,  as  these  will  prove  unsatisfactory,  designs  based 
on  unstable  matrices  A^  will  be  investigated. 

Consider  first  the  four  choices  of  which  give  rise  to  a  stable 

matrix  A^.  As  the  optimal  control  expended  little  energy  shaping  the 

*)  O  jfl  o  c  d  c  c  <1 

eigenvectors  »  ^ t*le  two  choices  (l2»^3»^2^  ’  ^2’^3’^3^  not  Pursued 

c  c  d  ,  c  c  d 

here.  The  other  two  choices  (A^,  A^  ,  Ag)  » (A^  ,A^ » retain  the  dominant 

coupled  frequency  pair  and  since  compensator  designs  based  on  the  first 

choice  were  unsuccessful,  only  the  latter  will  be  considered  here. 

dec 

Let  A^  ■  dg(X2 ,X^,X^) .  Based  on  a  desire  to  retain  the  dominant 

eigenvalues,  the  pole-placement  problem  is  solved  for  the  ordering 
c  c  d  d 

*6*^3*  4**5'  cumulations  are  performed  using  Procedure  1,  and  the 

data  are  given  in  Table  4.4.6.  In  the  following  discussion  eigenvalues 
without  superscripts  refer  to  spectra  in  Table  4.4.6.  At  the  first  stage 
of  the  algorithm  a  complex  pair  is  placed  at  -0.1+j0.1  in  the  expectation 
of  removing  the  two  small  eigenvalues  A^,Ag.  Unfortunately  the  resultant 
spectrum  of  A2  still  contains  a  small  eigenvalue  Ag  * -0.060.  The  compen¬ 
sator  based  on  this  solution  is  (D  ,K  have  been  scaled) 

o  z 


i 


i 


4 


% 


* 


f 
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Table  4.4.6. 


Data  for  pole-placement  problem  based  on  the  ordering: 
\C  AC  Ac  Ac 

a2»a4>a5»a5>a3 


'-2.0000  0 

.1961 

0.1830 

0.2746' 

4’ 

4.7500  -2 

.3665 

2.4576 

3.6869 

0.0000  0 

.0833 

-0.2445 

-0.1168 

* 

A  4" 

_  0.0000  0 

.0473 

2.0441 

-1.9338. 

1 

A  A 

"  0.0000  -0 

.1838 

-0.1715 

-0.2573* 

A2  - 

0.0000  2 

.2413 

2.0916 

3.1378 

A1 

0.0000  -0 

.0851 

-0.0794 

-0.1192 

_  0.0000  0 

.0655 

0.0612 

0.0918_ 

1- 

-7.8452 

0.2596 

7.3487" 

'-0.2336 

0.1484* 

- 

105.3412 

-0.0249 

-89.6134 

-0.1529 

0.1192 

3.3363 

-0.0595 

3.4031 

0.4484 

-0.0101 

-1.8903 

-0.0280 

-2.6205 

B  - 
o 

0.4941 

-0.7080 

7.3517 

-0.2596 

-7.8421 

0.2336 

-0.1484 

-89.6515 

0.0249 

-105.3800 

0.1540 

-0.1192 

3.4045 

0.0595 

3.3377 

-0.4484 

0.0101 

- 

-2.6216 

0.0280 

-1. 8914_ 

.-0.4940 

0.7080. 

-  [1.0000] 

f  .  p 
Z2  [o 

. 7072  | 

.  7070 _ J  8l" 

[3914.940  3918.477] 

g2  “ 

[1014. 

327  1017.972] 

The 

i  spectrum  of  is: 

The 

spectrum 

of  is: 

The 

spectrum  of  A^  is 

1. 

(-5.174, 

0.000) 

1. 

(-5.193, 

0.000) 

1. 

(-5.197, 

0.000) 

2. 

(-2.101, 

0.000) 

2. 

(-2.102, 

0.000) 

2. 

(-2.102, 

0.000) 

3. 

(-1.922, 

0.000) 

3. 

(-1.921, 

0.000) 

3. 

(-1.921, 

0.000) 

4. 

(-1.799, 

0.306) 

4. 

(-1.780, 

0.315) 

4. 

(-1.741, 

0.297) 

5. 

(-1.799, 

-0.306) 

5. 

(-1.780, 

-0.315) 

5. 

(-1.741, 

-0.297) 

6. 

(-0.230, 

0.000) 

6. 

(-0.112, 

0.104) 

6. 

(-0.152, 

0.153) 

7. 

(-0.039, 

0.000) 

7. 

(-0.112, 

-0.104) 

7. 

(-0.152, 

-0.153) 

8. 

(-0.024, 

0.000) 

8. 

(-0.060, 

0.000) 

8. 

(-0.037, 

0.000) 
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[-0.1904] 

D0 

-  [-182.623, 

0.00297,  182.634] 

f -762.3500 

-0.4577 

-764.9227] 

f  20.0423] 

(4.4.10) 

L  760.9804 

0.4577 

763.5529 J 

*  “  L-20.0423 j  * 

Continuing  the  solution  of  the  pole-placement  problem  because  of  the 

large  gains  in  K  and  the  small  eigenvalue  A  ,  f  is  chosen  at  the  second 
y  o  L 

stage  to  render  A^  uncontrollable  and  a  complex  pair  is  placed  at  0.13  +  j0.13. 

The  compensator  based  on  this  solution  is  (Dn,K  have  been  scaled) 

U  z 


no* 

-0.2322  0.0632] 

[“-214.281 

0.005176 

-214.298 

1  | 

.-0.0567  -1.9364J 

0*|_  -18.645 

0.002220 

-18.605 

J 

-841.9190  -0.4577 

-844.7776] 

r  20.131 

-9.015  ] 

(4.4.11)  j 

K  * 

y 

840.5499  0.4577 

843.4084  J  z 

*1.-20.131 

9.015 J 

a(HQ)  *  (-1.934,  -0.234). 

The  degree  of  stability  of  A^  has  decreased: 

Xg  * -0.037 , 

and  the 

gains  K  j 

y  !  1 

have  increased. 

Attempts  at  continuing  the  solution  of  the  pole-placement  problem 
to  increase  the  degree  of  stability  of  the  closed  loop  system  and  reduce 
the  feedback  gains  were  unsuccessful.  Therefore,  a  different  ordering  for 
Ar»Ap  was  considered.  Since  the  major  objection  to  the  design  above  is  the 
poor  stability  margin,  and  based  on  the  suspicion  that  the  small  real  pole 
in  cr(A^)  is  due  to  the  tie-line  interaction  mode  A^  departing  from  its 

c 

optimal  location,  A  is  selected  to  retain  A,.  The  remaining  two  eigenvalues 
r  o 

are  chosen  based  on  the  expectation  that  if  the  spectrum  of  contains  a 

complex  pair  near  + j2  then  this  will  correspond  to  a  departure  of  the 
c  c 

frequency  pair  A^.A^  from  their  optimal  locations.  Based  on  physical 
considerations  it  should  then  be  relatively  easy  to  shape  the  real  part  of 
this  complex  pair,  though  the  complex  part  should  be  insensitive  to  feedback. 
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From  an  Inspection  of  the  spectra  of  the  21  possible  matrices  A^ 

dec 

(this  data  is  not  reported  here),  the  choice  A  *  dg(A., A_, A,)  was  made. 

r  3  J  b 

The  ordering  for  the  pole-placement  problem  is  taken  to  be  A^,A^,A^,A^. 

The  pole-placement  problem  is  solved  using  Procedure  1  at  each  stage  and 

the  data  are  given  in  Table  4.4.7.  In  the  following  discussion  eigenvalues 

without  superscripts  refer  to  spectra  in  Table  4.4.7.  At  the  first  stage  a 

complex  pair  is  placed  at  the  location  of  A^,\^.  At  the  second  stage  f^  is 

chosen  to  render  A^  uncontrollable  and  a  complex  pair  is  again  placed  at 
c  c 

^4*^5"  At  t*ie  third  stage  f^  is  chosen  to  render  uncontrollable  the  complex 
pair  just  assigned  and  then  eigenvalues  are  placed  at  -2,-5.  At  the  last 
stage  f^  is  chosen  to  render  A^,A^,Ag  uncontrollable  and  a  complex  pair  is 
placed  at  -3  + j2. 

Three  compensators  are  defined  by  this  solution  to  the  pole- 

placement  problem.  Based  on  the  first  stage  of  the  solution  (Drt,k  have 

u  z 

been  scaled) 


Hq  -  [-6.960] 

Do  - 

(50.941  -0.239 

-63.913) 

f -608. 8379 

C. 2478 

755.8289] 

r  56.571 

y  l -638. 6288 

-0.2479 

791.7451J  * 

L  58.91 J 

Based  on  two  stages  of  the  solution  (D.  K  have  been  scaled) 

O’  z 

r-6.960  16.1161  [50.93;3  -0.239  -63.9051 

0  L  »  °1  0  0  0  J 


f-608.7381  0.2478  755.7395' 

y  L-638.5250  -0.2479  791.6521. 


|~56.57  -273.271 
z  158.91  -115. 18J 


(4.4.12) 


(4.4.13) 


Based  on  the  complete  solution  (D  ,K  have 

w  2 


been  scaled) 


J 
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Table  4.4.7.  Data  for  pole-placement  problem  based  on  the  ordering: 

'\C  AC  A^  AC  A1^  A*^ 

A3  9  A3  9  Ag  9  A^  ,  A^  9  A^  i  A3 


-A1  A2  1 

-2.0000 

0.1279  0.1194 

0.1791 

V 

A1  A1 

.1 

4.7500 

-4.8888  0.1037 

0.1556 

2  1 

Al  = 

0.0000 

-0.3099  -0.6114 

-0.6672 

-A1  Al. 

_  0.0000 

1.6633  3.5522 

0.3287. 

'  0.0000 

-0.1215  -0.1134 

-0.170l’ 

II 

0.0000 

-0.1010  -0.0943 

-0.1415 

0.0000 

0.4756  0.4439 

0.6659 

.  0.0000 

-0.3211  -0.2997 

-0.4496. 

'  -5.1168 

-0.5450 

4.8607" 

'  0.4571  -0. 

4775  -0 

.1525 

0.0870' 

-4.4470 

-0.5345 

4.0414 

-0.5199  0. 

5240  -0 

.1311 

0.0583 

19.0616 

1.7236 

-19.0248 

0.0104  0.0315  0 

.2584 

0.3075 

No“ 

-66.5339 

1.3320 

12.8315 

Bo  = 

0.0728  -0. 

1166  5 

.9765 

-10.3782 

4.8610 

0.5450 

-5.1170 

0.4765  0. 

4759  -0 

.1525 

0.0870 

4.0418 

0.5345 

-4.4473 

-0.5412  -0.5222  -0 

.1311 

0.0583 

-19.0253 

-1.7236 

19.0622 

0.0091  -0. 

0315  0 

.2584 

0.3075 

fl  = 

'  12.8455 

[1.0000] 

-1.3308 

-66.5481 _ 

f  J-1 

2  [O 

0  -a- 
0  0 

0  0 

0  0 
•  • 

_  0.0776  0. 

V 

1164  5.9765 

-0.9954 

0.0628 
.  0.0728- 

-10.3783  . 

hm 

[745.263 

-925.893] 

g2=  [0.122  - 

0.109]  g3  = 

[748.068 

-929 

.488] 

0.2210 
-0.1526 
0.7809 
-  0.5640 


P  = 


'  8.9730  0.0000 

41.1179  0.0000 

84.3020  0.0000 

L 21. 5255  0.0000 


-11.2749 

-50.9564 

-105.4287 

-27.2383. 


g4=  [38.169  -48.298] 
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Table  4.4.7.  continued 


The 

spectrum 

of  is: 

The 

spectrum 

of  is 

1. 

(-5.090, 

0.000) 

1. 

(-5.079, 

0.000) 

2. 

(-5.000, 

0.000) 

2. 

(-1.979, 

0.000) 

3. 

(-1.999, 

0.000) 

3. 

(-1.203, 

2.501) 

4. 

(-1.978, 

0.000) 

4. 

(-1.203, 

-2.501) 

5. 

(-0.167, 

0.000) 

5. 

(-1.093, 

0.196) 

6. 

(-0.112, 

0.000) 

6. 

(-1.093, 

-0.196) 

7. 

(  0.002, 

2.026) 

7. 

(-0.241, 

1.943) 

8. 

(  0.002, 

-2.026) 

8. 

(-0.241, 

-1.943) 

The 

spectrum 

of  A^  is: 

The 

spectrum 

of  A,  is 
4 

1. 

(-5.079, 

0.000) 

1. 

(-5.000, 

0.000) 

2. 

(-1.979, 

0.000) 

2. 

(-3.224, 

2.336) 

3. 

(-1.203, 

2.501) 

3. 

(-3.224, 

-2.336) 

4. 

(-1.203, 

-2.501) 

4. 

(-2.001, 

0.000) 

5. 

(-1.093, 

0.196) 

5. 

(-1.884, 

0.000) 

6. 

(-1.093, 

-0.196) 

-  6. 

(-0.304, 

0.000) 

7. 

(-0.241, 

1.943) 

7. 

(-0.241, 

1.943) 

8. 

(-0.241, 

-1.943) 

8. 

(-0.241, 

-1.943) 

The 

spectrum 

of  Aj  is: 

1. 

(-4.999, 

0.000) 

2. 

(-3.000, 

2.000) 

3. 

(-3.000, 

-2.000) 

4. 

(-2.001, 

0.000) 

5. 

(-1.842, 

0.000) 

6. 

(-0.546, 

0.000) 

7. 

(-0.214, 

1.943) 

8. 

(-0.241, 

-1.943) 
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'-9.1432 

0.2209 

1.7208  -3.2424' 

7.3655 

0.0949 

-8.994 " 

V 

0.1214 

-8.1721 

7.3545  -13.8580 

n  a 

33.439 

0.4032 

-40.400 

-0.0662 

-0.3129 

-4.5239  8.2957 

Do 

-7.968 

0.8692 

11.217 

.-0.0142 

-0.0517 

-1.0476  1.6279. 

.-0.837 

0.236 

1.741. 

II 

22.6798 

.-50.4229 

0.2478 

-0.2479 

-28.51041 

61. 0083  J  z" 

'8.276 

.8.334 

-7.668 

6.804 

5.957 

-11.387 

-10.206] 
22.475 J 

(4.4.14) 

The  first  order  compensator  gives  a  nice  total  closed  loop  spectrum 

but  has  the  same  disadvantage  encountered  before  of  large  feedback  gains. 

The  second  order  compensator  has  no  apparent  advantage  over  the  first  beyond 

c  c 

the  retention  of  the  optimal  pair  The  fourth  order  compensator 

however  has  nice  properties.  It  requires  gains  no  larger  than  60  as  compared 

to  20  for  the  full  state  feedback  solution.  Furthermore  it  retains  a  seven 

c  c  c  d  d  d  d 

dimensional  optimal  invariant  subspace  corresponding  to  (X  ,X  ,X  ,X  ,X  ,X  ,X  ) 

1  3  o  1  J  4  j 

c  c  c 

and  the  spectrum  of  contains  eigenvalues  near  X^jX^.X,,.  The  spectrum  of 
the  resultant  closed  loop  system  is  given  in  Table  4.4.8  and  Figure  4.4.2. 

The  compensator  is  also  open  loop  stable. 

In  summary  if  large  feedback  gains  are  tolerable  then  the  first 
order  compensator  defined  in  (4.4.12)  is  satisfactory,  however,  if  it  is 
desired  to  reduce  the  magnitude  of  the  feedback  gains  then  the  dimension  of 
the  compensator  must  be  increased.  The  fourth  order  controller  defined  by 
(4.4.14)  is  one  possibility  and  its  degree  compares  favorably  with  that  of 
the  Luenberger  reduced  order  observer  (n-r  *  8) .  It  is  noted  that  there  are 
many  other  possibilities  for  designing  satisfactory  controllers.  The 
analysis  presented  here  has  only  served  to  illustrate  the  design  methodology. 
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Table  4.4.8.  Spectra  of  closed  loop  system  under  dynamic  compensation 


For 

the  compensator 

|  For 

the  compensator 

given  in  (4 

.4.12) 

given 

in  (4. 

4.14) 

The 

spectrum  < 

af  A  is: 
c 

The 

spectrum 

of  A  is 
c 

1. 

(-9.171, 

0.000) 

1. 

( 

-9.171, 

0.000) 

2. 

(-5.079, 

0.000) 

2. 

( 

-9.171, 
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Figure  4.4.2.  Closed  loop  spectrum  using  the  compensator  given  in  (4.4.14) 
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CONCLUSIONS 

A  design  oriented  methodology  for  the  construction  of  suboptimal 
linear  quadratic  regulators  has  been  presented.  A  design  criterion  has 
been  taken  to  be  the  retention  of  as  many  optimal  eigenvectors  as  possible 
from  a  reference  state  feedback  regulator.  This  gives  rise  to  an  associated 
output  feedback  pole-placement  problem  both  in  designs  using  static  and 
dynamic  compensation.  For  the  latter  case  an  algorithm  has  been  given  which 
solves  this  pole-placement  problem,  implicitly  fixing  the  parameters  of  the 
controller,  and  determining  its  dimension  without  a  priori  assumptions.  It 
has  also  been  shown  that  the  methodology  may  be  extended  to  the  class  of 
stabilizable  systems.  Finally  the  design  methodology  has  been  illustrated 
with  three  nontrivial  examples. 
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